Ab initio molecular orbital study on the pairing and stacking interactions between nucleic acid bases in relation to the biological activities

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.

Study of the hydrogen bond in different orientations of adenine-thymine base pairs: An ab initio 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 evidence of the existence of a diazafulvene intermediate in the reaction pathway of imidazoleglycerol phosphate dehydratase: design of a novel and potent heterocycle structure for the inhibitor on the basis of the electronic structure-activity relationship study

The reaction mechanism of imidazoleglycerol phosphate dehydratase has not yet been clearly revealed. Structural comparison between inhibitors and the substrate IGP implicates that the reaction involves a diazafulvene intermediate. Here, we present evidence to support this hypothesis by investigating the electronic structure-enzyme inhibitory activity relationship on inhibitors with different heterocycles using 6-31G** level theory of the ab initio molecular orbital method. The calculation results showed that potent inhibitors can be distinguished from weak ones by the atomic charge density and by the energy levels of the highest occupied lone-pair orbital on the nitrogen atoms in the heterocycles. Furthermore, very good correlations (r2=0.8-0.9) were found between the charge density on the nitrogen atom and the inhibitory activity. It was also revealed that the diazafulvene is electronically similar to the potent inhibitors. Thus, these results strongly suggest the existence of the diazafulvene as an intermediate possessing tight-binding affinity to the enzyme. Based on the electronic structural similarity between the potent inhibitors and the proposed intermediate, a novel heterocycle was designed and predicted its inhibitory activity prior to the synthesis. Then, activity of synthesized inhibitors showed excellent agreement with this prediction. Hence, from the theoretical studies and experimental results, we conclude to obtain evidence of the hypothesis that the enzyme reaction proceeds via the diazafulvene intermediate.

Hydrogen bonding and stacking of DNA bases: a review of quantum-chemical ab initio studies

Ab initio quantum-chemical calulations with inclusion of electron correlation made since 1994 (such reliable calculations were not feasible before) significantly modified our view on interactions of nucleic acid bases. These calculations allowed to perform the first reliable comparison of the strength of stacked and hydrogen bonded pairs of nucleic acid bases, and to characterize the nature of the base-base interactions. Although hydrogen-bonded complexes of nucleobases are primarily stabilized by the electrostatic interaction, the dispersion attraction is also important. The stacked pairs are stabilized by dispersion attraction, however, the mutual orientation of stacked bases is determined rather by the electrostatic energy. Some popular theories of stacking were ruled out: The theory based on attractive interactions of polar exocyclic groups of bases with delocalized electrons of the aromatic rings (Bugg et al., Biopolymers 10, 175 (1971), and the pi-pi interactions model (C.A. Hunter, J. Mol. Biol. 230, 1025 (1993)). The calculations demonstrated that amino groups of nucleobases are very flexible and intrinsically nonplanar, allowing hydrogen-bond-like interactions which are oriented out of the plane of the nucleobase. Many H-bonded DNA base pairs are intrinsically nonplanar. Higher-level ab initio calculations provide a unique set of reliable and consistent data for parametrization and verification of empirical potentials. In this article, we present a short survey of the recent calculations, and discuss their significance and limitations. This summary is written for readers which are not experts in computational quantum chemistry.

Base stacking and hydrogen bonding in protonated cytosine dimer: the role of molecular ion-dipole and induction interactions

An ab initio quantum-chemical study of stacked and hydrogen-bonded protonated cytosine dimer has been carried out. The calculations were made using the second-order Moller-Plesset perturbational method (MP2) with a medium-sized polarized set of atomic orbitals. H-bonded as well as stacked protonated base pairs are more stable than the neutral base pairs. Two energy contributions not present in the neutral base pairs stabilize the protonated base pairs: the molecular ion - dipole interaction, and the induction interaction. The molecular ion - dipole stabilization dominates in base pairs with highly polar neutral monomers, such as the C...CH+ base pair. The induction interaction is not included in the commonly used empirical potentials, which do not reproduce the changes in intermolecular stabilization due to protonation. We demonstrate that the base stacking of several consecutive C...CH + pairs, as proposed for polycytidylic acid and i-DNA, is strongly repulsive. We also show that the intermolecular interactions strongly prefer protonation of adenine in protonated adenine-cytosine pairs.

Elucidation of primary radiation damage in DNA through application of ab initio molecular orbital theory

This review summarizes the results of recent ab initio molecular orbital calculations performed on DNA constituents that attempt to further our understanding of the radiation-induced damage to DNA. The results reviewed include calculations performed on the four individual DNA bases, the base pairs in gas phase and modelled aqueous phase, the deoxyribose moiety, and a portion of the sugar-phosphate backbone. The emphasis is on the electron affinities and ionization potentials of the radical species calculated under various conditions (i.e. gas phase, aqueous phase, proton transfer, base stacking), as it has been shown that the initial ion radical distribution is largely a function of these two properties. Theoretical studies of the electronic excited states of the individual bases and radioprotection of the biomolecule by various thiol compounds are also reviewed. Finally, a summary is provided to allow for further elaboration of the current model for radiation damage to DNA and to show the present advantages and limitations of ab initio theory in the investigation of such processes.