Intermolecular Proton Transfer in Anionic Complexes of Uracil with Alcohols

Proton transfer from water to aromatic N-heterocyclic anions from DFT-MD simulations

The phenomenon of proton transfer from water to six N-heterocyclic anions and free energy landscapes of this process are studied using both electronic structure calculations and first principles molecular metadynamics simulations. Our investigation involves microhydrated and aqueous phase interaction of water with six aromatic heterocyclic anions relevant to chemistry and biology: imidazolide, pyrrolide, benzimidazolide, 2-cyanopyrrolide, indolide, and indazolide. The basic structures of all these heterocyclic anions differ by substituted functional groups as well as fused rings. We study the proton transfer reaction and the minimum number of required water molecules for the reaction in hydrated microclusters. We find out that at least four water molecules are necessary for hydrated clusters to facilitate the intracluster proton transfer reaction from water to anions except for pyrrolide, for which this magic number is 3. To obtain the reaction free energy and activation barrier of the proton transfer process in an aqueous solution, the metadynamics method based first principles molecular dynamics simulations were performed. The complete proton transfer was observed in aqueous solutions for all the anions. The water molecule directly involved in proton transfer becomes acidic due to the cooperative effect of neighboring water molecules. From the metadynamics simulation, we obtain the values of activation barrier for the proton transfer processes from neutral water to anions, and the highest activation barrier is obtained for benzimidazolide, whereas the lowest activation barrier is obtained for pyrrolide. The structures and free energy profiles of the process for all the anions are discussed, and a comparative outlook of the study is presented here.

The Role of Proton Transfer on Mutations

Hydrogen bonds play a critical role in nucleobase studies as they encode genes, map protein structures, provide stability to the base pairs, and are involved in spontaneous and induced mutations. Proton transfer mechanism is a critical phenomenon that is related to the acid-base characteristics of the nucleobases in Watson-Crick base pairs. The energetic and dynamical behavior of the proton can be depicted from these characteristics and their adjustment to the water molecules or the surrounding ions. Further, new pathways open up in which protonated nucleobases are generated by proton transfer from the ionized water molecules and elimination of a hydroxyl radical in this review, the analysis will be focused on understanding the mechanism of untargeted mutations in canonical, wobble, Hoogsteen pairs, and mutagenic tautomers through the non-covalent interactions. Further, rare tautomer formation through the single proton transfer (SPT) and the double proton transfer (DPT), quantum tunneling in nucleobases, radiation-induced bystander effects, role of water in proton transfer (PT) reactions, PT in anticancer drugs-DNA interaction, displacement and oriental polarization, possible models for mutations in DNA, genome instability, and role of proton transfer using kinetic parameters for RNA will be discussed.

Intramolecular electron-induced proton transfer and its correlation with excited-state intramolecular proton transfer

Electron-induced proton transfer depicts the proton motion coupled with the attachment of a low-energy electron to a molecule, which helps to understand copious fundamental chemical processes. Intramolecular electron-induced proton transfer is a similar process that occurs within a single molecule. To date, there is only one known intramolecular example, to the best of our knowledge. By studying the 10-hydroxybenzo[h]quinoline and 8-hydroxyquinoline molecules using anion photoelectron spectroscopy and density functional theory, and by theoretical screening of six other molecules, here we show the intramolecular electron-induced proton transfer capability of a long list of molecules that meanwhile have the excited-state intramolecular proton transfer property. Careful examination of the intrinsic electronic signatures of these molecules reveals that these two distinct processes should occur to the same category of molecules. Intramolecular electron-induced proton transfer could have potential applications such as molecular devices that are responsive to electrons or current.

A delicate case of unidirectional proton transfer from water to an aromatic heterocyclic anion

One of the hydroxyl modes of the surrounding water molecules with the lowest stretching frequency facilitates the proton transfer from water to an anion.

The onset of electron-induced proton-transfer in hydrated azabenzene cluster anions

The prospect that protons from water may be transferred to N-heterocyclic molecules due to the presence of an excess electron is studied in hydrated azabenzene cluster anions using spectroscopy and computational chemistry.

Molecular recognition of nucleotides in water by scorpiand-type receptors based on nucleobase discrimination

The detection of nucleotides is of crucial importance because they are the basic building blocks of nucleic acids. Scorpiand-based polyamine receptors functionalized with pyridine or anthracene units are able to form stable complexes with nucleotides in water, based on coulombic, π-π stacking, and hydrogen-bonding interactions. This behavior has been rationalized by means of an exploration with NMR spectroscopy and DFT calculations. Binding constants were determined by potentiometry. Fluorescence spectroscopy studies have revealed the potential of these receptors as sensors to effectively and selectively distinguish guanosine-5'-triphosphate (GTP) from adenosine-5'-triphosphate (ATP).

Photoelectron spectroscopy and density functional theory studies on the uridine homodimer radical anions

We report the photoelectron spectrum (PES) of the homogeneous dimer anion radical of uridine, (rU)(2)(●-). It features a broad band consisting of an onset of ∼1.2 eV and a maximum at the electron binding energy (EBE) ranging from 2.0 to 2.5 eV. Calculations performed at the B3LYP∕6-31++G∗∗ level of theory suggest that the PES is dominated by dimeric radical anions in which one uridine nucleoside, hosting the excess charge on the base moiety, forms hydrogen bonds via its O8 atom with hydroxyl of the other neutral nucleoside's ribose. The calculated adiabatic electron affinities (AEAGs) and vertical detachment energies (VDEs) of the most stable homodimers show an excellent agreement with the experimental values. The anionic complexes consisting of two intermolecular uracil-uracil hydrogen bonds appeared to be substantially less stable than the uracil-ribose dimers. Despite the fact that uracil-uracil anionic homodimers are additionally stabilized by barrier-free electron-induced proton transfer, their relative thermodynamic stabilities and the calculated VDEs suggest that they do not contribute to the experimental PES spectrum of (rU)(2)(●-).

Effect of the methylation of uracil and/or glycine on their mutual interaction

In order to simulate the hydrogen bonding and proton transfer (PT) in protein-DNA/RNA interactions, a series of simplified models were employed and investigated in the gas phase. These models included various neutral, anionic and cationic glycine-uracil dimers, and their methylated derivatives generated by the mono- or dimethylation of glycine and/or uracil moieties of the dimer. The results reveal that the only process that can occur in the neutral complexes is a double-PT process leading to proton exchange between the two moieties (i.e., point mutation). The first methyl substitute can reduce the activation energy of the PT process and thus promote the isomerization of the two moieties; further methylation can reduce the isomerization in only some of the cases. In the anionic complexes, only the one-way PT (i.e., amino acid → nucleic acid base) process is energetically favorable, and this PT process is an interesting barrier-free one (BFPT), with the attached electron locating itself at the base moiety. Methylation will disfavor BFPT, but it cannot alter the nature of BFPT. In the cationic complexes, three different PT processes can occur. These processes can transform mutually by adjusting either or both of the methylated sites and methyl number, indicating that the methylation can regulate the dynamics of these PT processes.

Photoelectron spectroscopic studies of 5-halouracil anions

The parent negative ions of 5-chlorouracil, UCl(-) and 5-fluorouracil, UF(-) have been studied using anion photoelectron spectroscopy in order to investigate the electrophilic properties of their corresponding neutral halouracils. The vertical detachment energies (VDE) of these anions and the adiabatic electron affinities (EA) of their neutral molecular counterparts are reported. These results are in good agreement with the results of previously published theoretical calculations. The VDE values for both UCl(-) and UF(-) and the EA values for their neutral molecular counterparts are much greater than the corresponding values for both anionic and neutral forms of canonical uracil and thymine. These results are consistent with the observation that DNA is more sensitive to radiation damage when thymine is replaced by halouracil. While we also attempted to prepare the parent anion of 5-bromouracil, UBr(-), we did not observe it, the mass spectrum exhibiting only Br(-) fragments, i.e., 5-bromouracil apparently underwent dissociative electron attachment. This observation is consistent with a previous assessment, suggesting that 5-bromouracil is the best radio-sensitizer among these three halo-nucleobases.

Structures of alkali metal ion-adenine complexes and hydrated complexes by IRMPD spectroscopy and electronic structure calculations

Complexes between adenine and the alkali metal ions Li(+), Na(+), K(+), and Cs(+) have been investigated by infrared multiple photon dissociation (IRMPD) spectroscopy between 2800 and 3900 cm(-1), as have some singly hydrated complexes. The IRMPD spectra clearly show the N-H stretching and the NH(2) symmetric and asymmetric stretching vibrations of adenine; and for the solvated ions, the O-H stretching vibrations are observed. These experimental spectra were compared with those for a variety of possible structures, including both A9 (A9 refers to the tautomer where hydrogen is on the nitrogen in position 9 of adenine, see Scheme 1) and A7 adenine tautomers, computed using B3-LYP/6-31+G(d,p). By comparing the experimental and the simulated spectra it is possible to rule out various structures and to further assign structures to the species probed in these experiments. Single-point calculations on the B3-LYP/6-31+G(d,p) geometries have been performed at MP2/6-311++G(2d, p) to obtain good estimates of the relative thermochemistries for the different structures. In all cases the computed IR spectrum for the lowest energy structure is consistent with the experimental IRMPD spectrum, but in some cases structural assignment cannot be confirmed based solely upon comparison with the experimental spectra so computed thermochemistries can be used to rule out high-energy structures. On the basis of the IRMPD spectra and the energy calculations, all adenine-M(+) and adenine-M(+)-H(2)O are concluded to be composed of the A7 tautomer of adenine, which is bound to the cations in a bidentate fashion through N3 and N9 (see Scheme 1 for numbering convention). For the hydrated ions water binds directly to the metal ion through oxygen, as would be expected since the metal contains most positive charge density. For the hydrated lithium cation-bound adenine dimer, the water molecule is concluded to be hydrogen bonded to a free basic site of one of the adenine monomers, which is also bound to the lithium cation. Experimental and theoretical results on adenine-Li(+)-H(2)O suggest that the electrosprayed adenine-Li(+) resembles the lowest-energy solution phase ion rather than the lowest-energy gas-phase ion, which is the imine form.

Combinatorial-computational-chemoinformatics (C3) approach to finding and analyzing low-energy tautomers

Finding the most stable tautomer or a set of low-energy tautomers of molecules is critical in many aspects of molecular modelling or virtual screening experiments. Enumeration of low-energy tautomers of neutral molecules in the gas-phase or typical solvents can be performed by applying available organic chemistry knowledge. This kind of enumeration is implemented in a number of software packages and it is relatively reliable. However, in esoteric cases such as charged molecules in uncommon, non-aqueous solvents there is simply not enough available knowledge to make reliable predictions of low energy tautomers. Over the last few years we have been developing an approach to address the latter problem and we successfully applied it to discover the most stable anionic tautomers of nucleic acid bases that might be involved in the process of DNA damage by low-energy electrons and in charge transfer through DNA. The approach involves three steps: (1) combinatorial generation of a library of tautomers, (2) energy-based screening of the library using electronic structure methods, and (3) analysis of the information generated in step (2). In steps 1-3 we employ combinatorial, computational and chemoinformatics techniques, respectively. Therefore, this hybrid approach is named "Combinatorial*Computational*Chemoinformatics", or just abbreviated as C(3) (or C-cube) approach. This article summarizes our developments and most interesting methodological aspects of the C(3) approach. It can serve as an example how to identify the most stable tautomers of molecular systems for which common chemical knowledge had not been sufficient to make definite predictions.

Valence anion of thymine in the DNA pi-stack

Most of theoretical data on the stability of radical anions supported by nucleic acid bases have been obtained for anions of isolated nucleobases, their nucleosides, or nucleotides. This approach ignores the hallmark forces of DNA, namely, hydrogen bonding and pi-stacking interactions. Since these interactions might be crucial for the electron affinities of nucleobases bound in DNA, we report for the first time on the stability of the thymine valence anion in trimers of complementary bases possessing the regular B-DNA geometry but differing in base sequence. In order to estimate the energetics of electron attachment to a trimer, we developed a thermodynamic cycle employing all possible two-body interaction energies in the neutral and anionic duplex as well as the adiabatic electron affinity of isolated thymine. All calculations were carried out at the MP2 level of theory with the aug-cc-pVDZ basis set. The two-body interaction energies were corrected for the basis set superposition error, and in benchmark systems, they were extrapolated to the basis set limit and supplemented with correction for higher order correlation terms calculated at the CCSD(T) level. We have demonstrated that the sequence of nucleic bases has a profound effect on the stability of the thymine valence anion: the anionic 5'-CTC-3' (6.0 kcal/mol) sequence is the most stable configuration, and the 5'-GTG-3' (-8.0 kcal/mol) trimer anion is the most unstable species. On the basis of obtained results, one can propose DNA sequences that are different in their vulnerability to damage by low energy electron.

Barrier-free proton transfer in the valence anion of 2'-deoxyadenosine-5'-monophosphate. II. A computational study

The propensity of four representative conformations of 2(')-deoxyadenosine-5(')-monophosphate (5(')-dAMPH) to bind an excess electron has been studied at the B3LYP6-31++G(d,p) level. While isolated canonical adenine does not support stable valence anions in the gas phase, all considered neutral conformations of 5(')-dAMPH form adiabatically stable anions. The type of an anionic 5(')-dAMPH state, i.e., the valence, dipole bound, or mixed (valence/dipole bound), depends on the internal hydrogen bond(s) pattern exhibited by a particular tautomer. The most stable anion results from an electron attachment to the neutral syn-south conformer. The formation of this anion is associated with a barrier-free proton transfer triggered by electron attachment and the internal rotation around the C4(')-C5(') bond. The adiabatic electron affinity of the a_south-syn anion is 1.19 eV, while its vertical detachment energy is 1.89 eV. Our results are compared with the photoelectron spectrum (PES) of 5(')-dAMPH(-) measured recently by Stokes et al., [J. Chem. Phys. 128, 044314 (2008)]. The computational VDE obtained for the most stable anionic structure matches well with the experimental electron binding energy region of maximum intensity. A further understanding of DNA damage might require experimental and computational studies on the systems in which purine nucleotides are engaged in hydrogen bonding.

Effects of Microsolvation on the Adenine−Uracil Base Pair and Its Radical Anion: Adenine−Uracil Mono- and Dihydrates

Microhydration effects upon the adenine-uracil (AU) base pair and its radical anion have been investigated by explicitly considering various structures of their mono- and dihydrates at the B3LYP/DZP++ level of theory. For the neutral AU base pair, 5 structures were found for the monohydrate and 14 structures for the dihydrate. In the lowest-energy structures of the neutral mono- and dihydrates, one and two water molecules bind to the AU base pair through a cyclic hydrogen bond via the N(9)-H and N(3) atoms of the adenine moiety, while the lowest-lying anionic mono- and dihydrates have a water molecule which is involved in noncyclic hydrogen bonding via the O4 atom of the uracil unit. Both the vertical detachment energy (VDE) and adiabatic electron affinity (AEA) of the AU base pair are predicted to increase upon hydration. While the VDE and AEA of the unhydrated AU pair are 0.96 and 0.40 eV, respectively, the corresponding predictions for the lowest-lying anionic dihydrates are 1.36 and 0.75 eV, respectively. Because uracil has a greater electron affinity than adenine, an excess electron attached to the AU base pair occupies the pi* orbital of the uracil moiety. When the uracil moiety participates in hydrogen bonding as a hydrogen bond acceptor (e.g., the N(6)-H(6a)...O(4) hydrogen bond between the adenine and uracil bases and the O(w)-H(w)...N and O(w)-H(w)...O hydrogen bonds between the AU pair and the water molecules), the transfer of the negative charge density from the uracil moiety to either the adenine or water molecules efficiently stabilizes the system. In addition, anionic structures which have C-H...O(w) contacts are energetically more favorable than those with N-H...O(w) hydrogen bonds, because the C-H...O(w) contacts do not allow the unfavorable electron density donation from the water to the uracil moiety. This delocalization effect makes the energetic ordering for the anionic hydrates very different from that for the corresponding neutrals.

Investigation of proton transport tautomerism in clusters of protonated nucleic acid bases (cytosine, uracil, thymine, and adenine) and ammonia by high-pressure mass spectrometry and ab initio calculations

The energetics of the ion-molecule interactions and structures of the clusters formed between protonated nucleic acid bases (cytosine, uracil, thymine, and adenine) and ammonia have been studied by pulsed ionization high-pressure mass spectrometry (HPMS) and ab initio calculations. For protonated cytosine, uracil, thymine, and adenine with ammonia, the measured enthalpies of association with ammonia are -21.7, -27.9, -22.1, and -17.5 kcal mol-1, respectively. Different isomers of the neutral and protonated nucleic acid bases as well as their clusters with ammonia have been investigated at the B3LYP/6-31+G(d,p) level of theory, and the corresponding binding energetics have also been obtained. The potential energy surfaces for proton transfer and interconversion of the clusters of protonated thymine and uracil with ammonia have been constructed. For cytosine, the experimental binding energy is in agreement with the computed binding energy for the most stable isomer, CN01-01, which is derived from the enol form of protonated cytosine, CH01, and ammonia. Although adenine has a proton affinity similar to that of cytosine, the binding energy of protonated adenine to ammonia is much lower than that for protonated cytosine. This is shown to be due to the differing types of hydrogen bonds being formed. Similarly, although uracil and thymine have similar structures and proton affinities, the binding energies between the protonated species and ammonia are different. Strikingly, the addition of a single methyl group, in going from uracil to thymine, results in a significant structural change for the most stable isomers, UN01-01 and TN03-01, respectively. This then leads to the difference in their measured binding energies with ammonia. Because thymine is found only in DNA while uracil is found in RNA, this provides some potential insight into the difference between uracil and thymine, especially their interactions with other molecules.

Hydrogen-Atom Abstraction from the Adenine−Uracil Base Pair†

The hydrogen-abstracted radicals from the adenine-uracil (AU) base pair have been studied at the B3LYP/DZP++ level of theory. The A(N9)-U and A-U(N1) radicals, which correspond to hydrogen-atom abstraction at the adenine N9 and uracil N1 atoms, respectively, were predicted to be the two lowest-lying among the nine (AU-H) radicals studied in this study. The removal of the amino hydrogen of the adenine moiety that forms a hydrogen bond with the uracil O4 atom in the AU pair resulted in radical A(N6a)-U, which has the smallest base-pair dissociation energy, 5.9 kcal mol(-1). This radical is more likely to dissociate into the two isolated bases than to recover the hydrogen bond with the O4 atom through N6-H bond rotation along the C6-N6 bond. In general, the radicals generated by C-H bond breaking were higher in energy than those arising from N-H bond cleavage, because the unpaired electrons in the carbon-centered radicals were mainly localized on the carbon atom from which the hydrogen atom was removed. However, the highest-lying radical was found to arise from removal of the N3 hydrogen of uracil. The most remarkable structural feature of this radical is a very short C-H...O distance of 2.094 A, consistent with a substantial hydrogen bond. Although this radical lost the N1...H-N3 hydrogen bond between the two bases, its dissociation energy was predicted to be 12.9 kcal mol(-1), similar to that of the intact AU base pair. This is due to the transfer of electron density from the adenine N1 atom to the uracil N3 atom.

Proton-transfer and H2-elimination reactions of main-group hydrides EH4- (E = B, Al, Ga) with alcohols

The reaction of the isostructural anions of group 13 hydrides EH4- (E = B, Al, Ga) with proton donors of different strength (CH3OH, CF3CH2OH, and CF3OH) was studied with different theoretical methods [DFT/B3LYP and second-order Møller-Plesset (MP2) using the 6-311++G(d,p) basis set]. The results show the general mechanism of the reaction: the dihydrogen-bonded (DHB) adduct (EH...HO) formation leads through the activation barrier to the next concerted step of H2 elimination and alkoxo product formation. The structures, interaction energies (calculated by different approaches including the energy decomposition analysis), vibrational E-H modes, and electron-density distributions were analyzed for all of the DHB adducts. The transition state (TS) is the dihydrogen complex stabilized by a hydrogen bond with the anion [EH3(eta2-H2)...OR-]. The single exception is the reaction of BH4- with CF3OH exhibiting two TSs separated by a shallow minimum of the BH3(eta2-H2)...OR- intermediate. The structures and energies of all of the species were calculated, leading to the establishment of the potential energy profiles for the reaction. A comparison is made with the mechanism of the proton-transfer reaction to transition-metal hydrides. The solvent influence on the stability of all of the species along the reaction pathway was accounted for by means of polarizable conductor calculation model calculations in tetrahydrofuran (THF). Although in THF the DHB intermediates, the TSs, and the products are destabilized with respect to the separated reactants, the energy barriers for the proton transfer are only slightly affected by the solvent. The dependence of the energies of the DHB complexes, TSs, and products as well as the energy barriers for the H2 release on the central atom and the proton donor strength is also discussed.

Valence Anions in Complexes of Adenine and 9-Methyladenine with Formic Acid: Stabilization by Intermolecular Proton Transfer

Photoelectron spectra of adenine-formic acid (AFA(-)) and 9-methyladenine-formic acid (MAFA(-)) anionic complexes have been recorded with 2.540 eV photons. These spectra reveal broad features with maxima at 1.5-1.4 eV that indicate formation of stable valence anions in the gas phase. The neutral and anionic complexes of adenine/9-methyladenine and formic acid were also studied computationally at the B3LYP, second-order Møller-Plesset, and coupled-cluster levels of theory with the 6-31++G** and aug-cc-pVDZ basis sets. The neutral complexes form cyclic hydrogen bonds, and the most stable dimers are bound by 17.7 and 16.0 kcal/mol for AFA and MAFA, respectively. The theoretical results indicate that the excess electron in both AFA(-) and MAFA(-) occupies a pi* orbital localized on adenine/9-methyladenine, and the adiabatic stability of the most stable anions amounts to 0.67 and 0.54 eV for AFA(-) and MAFA(-), respectively. The attachment of the excess electron to the complexes induces a barrier-free proton transfer (BFPT) from the carboxylic group of formic acid to a N atom of adenine or 9-methyladenine. As a result, the most stable structures of the anionic complexes can be characterized as neutral radicals of hydrogenated adenine (9-methyladenine) solvated by a deprotonated formic acid. The BFPT to the N atoms of adenine may be biologically relevant because some of these sites are not involved in the Watson-Crick pairing scheme and are easily accessible in the cellular environment. We suggest that valence anions of purines might be as important as those of pyrimidines in the process of DNA damage by low-energy electrons.

On the Unusual Stability of Valence Anions of Thymine Based on Very Rare Tautomers: A Computational Study

We characterized anionic states of thymine using various electronic structure methods, with the most accurate results obtained at the CCSD(T)/aug-cc-pVDZ level of theory followed by extrapolations to complete basis set limits. We found that the most stable anion in the gas phase is related to an imino-oxo tautomer, in which the N1H proton is transferred to the C5 atom. This valence anion, aT(c5)(nl), is characterized by an electron vertical detachment energy (VDE) of 1251 meV and it is adiabatically stable with respect to the canonical neutral nT(can) by 2.4 kcal/mol. It is also more stable than the dipole-bound (aT(dbs)(can)), and valence anion aT(val)(can) of the canonical tautomer. The VDE values for aT(dbs)(can)and T(val)(can) are 55 and 457 meV, respectively. Another, anionic, low-lying imino-oxo tautomer with a VDE of 2458 meV has a proton transferred from N3H to C5 aT(c5)(n3). It is less stable than aT(val)(can) by 3.3 kcal/mol. The mechanism of formation of anionic tautomers with the carbons C5 or C6 protonated may involve intermolecular proton transfer or dissociative electron attachment to the canonical neutral tautomer followed by a barrier-free attachment of a hydrogen atom to C5. The six-member ring structure of the anionic tautomers with carbon atoms protonated is unstable upon an excess electron detachment. Within the PCM hydration model, the low-lying valence anions become adiabatically bound with respect to the canonical neutral; becomes the most stable, being followed by aT(c5)(nl), aT(c5)(n3), aT(can), and aT(c5)(nl).