Structure and Stability of Chemically Modified DNA Bases: Quantum Chemical Calculations on 16 Isomers of Diphosphocytosine
We studied for the first time 16 tautomers/rotamers of diphosphocytosine by four computational methods. Some of these tautomers/rotamers are isoenergetic although they have different structures. High-level electron correlation MP2 and MP4(SDQ) ab initio methods and density functional methods employing a B3LYP and the new M06-2X functional were used to study the structure and relative stability of 16 tautomers/rotamers of diphosphocytosine. The dienol tautomers of diphosphocytosine are shown to be much more stable than the keto-enol and diketo forms. The tautomers/rotamers stability could be ranked as PC3 = PC12 < PC2 = PC11 < PC1 < PC10 < PC8 < PC9 < PC15 < PC16 < PC6 ~ PC7 < PC13 < PC4 ~ PC14 < PC5. This stability order was discussed in the light of stereo and electronic factors. Solvation effect has been modeled in a high dielectric solvent, water using the polarized continuum model (PCM). Consideration of the solvent causes some reordering of the relative stability of diphosphocytosine tautomers: PC3 ~ PC12 ~ PC2 ~ PC11 < PC1 < PC10 < PC8 < PC9 < PC15 ~ PC16 < PC13 < PC6 ~ PC7 ~ PC14 < PC4 ~ PC5.
DNA is a naturally occurring biological macromolecule, containing thousands of nucleic acid bases, and it is of prime importance in genetic determination . The five nucleic acid bases, cytosine, thymine, uracil, adenine, and guanine, found in DNA and RNA control the replication of DNA, store information required to synthesize proteins, and translate this information to the protein. Tautomerism is a well-known phenomenon occurring in nucleic acid bases [2–16], in which proton transfer from the heterocyclic ring center to an exocyclic oxo- or imino- group leads to the formation of either an –OH or an –NH2 groups. These processes are known as keto-enol or imino-amino tautomerism, respectively. Tautomerism, a sort of isomerism, plays an important role in organic chemistry, medicinal chemistry, pharmacology, and molecular biology. Tautomerism partially explains the structure of nucleic acids and their mutations . In DNA bases, tautomerization results in altered base pairing configurations or mispairing due to changes in hydrogen-bonding capabilities. DNA mutations are likely to be caused by such alterations. Cytosine is one of the building pyrimidine nucleobases of RNA. A large amount of experimental [18–34] and theoretical work [35–58] has been carried out in order to elucidate the structure of cytosine and its tautomers. Chemically modified bases have attracted extensive interest due to their numerous pharmacological, biochemical, and biological capabilities [59–62]. Replacing nitrogen in pyrimidine and purine DNA bases by the biological element of interest, phosphorus, may produce new synthesizable bases of special interest. Modification of cytosine base to phosphorus-containing cytosine analogues, diphosphocytosine, is a good start point in this topic. Also, their existence in many tautomeric forms, like other nucleoside bases, seems to be crucial in order to explain the mutation occurring during DNA duplication [63–66]. Searching the literatures confirms that no previous studies include all numbers of cytosine tautomers. This is the first time to include 16 isomers of cytosine. Similar to cytosine, the diphosphocytosine should also have 16 tautomers. There are many interesting questions on the diphosphocytosine tautomers concerning, for example, their relative stability order, energy criteria, proton transfer processes, electronic structure, and chemical properties. In an attempt to answer some of these questions using computational studies, some DFT methods (B3LYP, M06-2X) and ab initio MP2 and MP4(SDQ) methods are chosen to study the stability order and tautomerization processes of all the possible 16 diphosphocytosine tautomers in the gas phase as well as in water. Compared to the popular B3LYP method, studies of tautomerism using the recent method M06-2X are very rare and it is worthwhile to consider them here in this work to assess their accuracy compared to MP4(SDQ) data.
Among cytosine tautomers, experimental investigations performed in gas phase, solid state, or in solution indicate that the canonic amino-oxo tautomer of cytosine is the most stable one . One of our interesting computational results [67, 68] of diphosphocytosine tautomers is the reverse stability order as compared to the parent analogue, cytosine. The origin of this stability order needs to be interpreted as it is not fully discussed. One of the main objectives of this work is to investigate this behavior on the basis of the unique bonding properties of the phosphorous atom.
2. Computational Methods
All quantum chemical calculations were performed using the Gaussian 09 suite of programs . Geometry optimizations for all tautomers have been performed using ab initio density functional theory (DFT) at the B3LYP functional [70–72] and the new M06-2X functional in conjunction with the 6-31G(d) basis set. For each stationary point, we carried out a force constant harmonic frequency calculation at the same levels to characterize their nature as minima or transition states and to correct energies for zero-point energy and thermal contribution. The frequencies are scaled by a factor of 0.98. The vibrational modes were animated using the ChemCraft program . Natural charges were computed within full Natural Bond Orbital Analysis, using NBO implemented in Gaussian 09 . The vibrational modes were examined by using the ChemCraft program . Partial charge distributions were calculated using the natural population analysis (NPA) method . Solvation effect has been modeled in water using the polarized continuum model (PCM) of Tomasi et al. [75–78] at the B3LYP/6-31+G(d,p) level for all tautomers/rotamers.
3. Results and Discussion
3.1. Geometrical Aspects of Diphosphocytosine Tautomeric Forms
Geometry of the 16 tautomers and rotamers was optimized at four different levels. Optimized structures of the 16 diphosphocytosine tautomers (at MP4(SDQ)/6-31+G(d)) are presented in Figure 1. For the notations of the tautomers and rotamers of diphosphocytosine, we use PC1, PC2, PC3, 16 for the sixteen forms. All structures with their notations and atomic numbering are shown in Figure 1. Structure geometrical parameters of all the 16 tautomers are listed in Table 1. Bond length changes among tautomers show the same trend as expected from formal bond orders. In tautomers PC2, PC3, PC11, and PC12, all CC and CP bond lengths are in the range characteristic for aromatic rings. C1-C2 bond length changes among tautomers depending on the amount of π-bond character the tautomer has. In tautomers PC2, PC3, PC11, and PC12 this bond distance is about 1.40 Å while in the rest of the tautomers this distance increased as shown in Table 1. In all the dienol forms, the P3-C2 and P3-C4 and also P5-C6 and P5-C4 distances are nearly similar while those for the other keto-enol and diketo forms show considerable differences of about 0.147 Å. This reflects the extent of conjugation and the increasing aromatic character of the six-membered ring. An inspection of the dihedral angles in Table 1 proves that all the tautomers/rotamers structures are nonplanar. It is reported that  results of electron correlation calculations (MP2/6-31G*) give nonplanar geometries for several nucleotide bases specifically, for cytosine, which showed a significant pyramidalization of the amino group from MP2/6-31G* calculations. Using enormous computer resources, a complete MP2/6-31G** frequency calculation was carried out  for cytosine, confirming that the planar structure is a saddle point. Bond angles in the ring show quite dramatic changes between tautomers and rotamers, by up to more than 10°. Owing to the lack of polarization functions on hydrogens, geometry parameters of the latter are expected to show larger differences, but these less important parameters were not listed in.
In the dienol form of diphosphocytosine, there are three double bonds giving it the six π-electrons needed to gain its aromatic character as in case of benzene. The phosphorous atom has an additional lone pair of electrons but not involved in the conjugated system. Once the keto form is formed in the diphosphocytosine, one resonance form is lost as the phosphorus atom is strongly pyramidal in the keto forms and its lone pairs are no longer available for conjugation. The extent of pyramidalization emerges from the sum of bond angles around the phosphorus atom. An extra destabilization occurs upon forming the second keto tautomer. Accordingly, diketo form is accompanied by destruction of conjugated π-electrons system in the ring of the diphosphocytosine. In case of cytosine, on the other hand, formation of keto or diketo forms still keeps planarity at the nitrogen atoms where they still able to resonate with the adjacent double bonds. Aromaticity is lost in both cytosine and diphosphocytosine, but the number of resonance structures decreased in the case of the latter system.
3.2. Relative Energies and Thermodynamic Parameters of Tautomers
Since our major objective of this study was to obtain information on the relative stabilities of diphosphocytosine tautomers and rotamers, we tried to go as high as possible, with the level of theory and use both DFT and high correlation ab initio MP2 and MP4(SDQ) methods. The results on the relative stabilities of diphosphocytosine tautomers and rotamers using different computational levels are listed in Table 2. For a better view, this is also shown in Figure 2. The results of the four computational methods are consistent with slight differences. The 16 diphosphocytosine tautomers and rotamers could be classified, on the basis of total relative energy, into three groups: (1) low-energy tautomers group containing four tautomers namely, PC2, PC3, PC11, and PC12. These tautomers: are those satisfying the aromaticity of the pyrimidine ring (the benzene-like structure), (2) moderate-energy tautomers group containing also four tautomers, PC1, PC8, PC9, and PC10. (3) high energy tautomers group containing the rest eight tautomers, PC4, PC5, PC6, PC7, PC13, PC14, PC15, and PC16. DFT methods predict PC2 to be the most stable lowest-energy tautomer. High correlation ab initio methods, on the other hand, predict the most stable is the hydroxy-amino form PC3. Energy difference is less than 0.50 kcal/mol. The structural difference between the two tautomers/rotamers is the proximity of the positively charged H to the positively charged phosphorus atom. The same argument is applied to PC11 and PC12, where the positively charged H is repulsively interacting with P3 (0.503, Table 3) in PC11 and P5 (0.531, Table 3) in PC12. Previous studies suggest that the canonical amino-oxo structure is the most stable one . PC5 is the most unstable highest-energy tautomer as predicted by all methods. This could be due to the destruction of the aromaticity and also attributed to the repulsive interaction of the positively charged H and the great positive charge accommodated on P3 (0.657, Table 3). As could be seen from Table 2 that including zero-point energy and thermal energy correction has a considerable effect on the molecular energy. The corrected energy is lower than the uncorrected by about 4.0 kcal/mol. Also, considering these corrections causes reordering of the tautomers stability. Figure 3 shows the MP2-corrected for zero-point and thermal energy versus uncorrected one for all the 16 diphosphocytosine tautomers. Based on the corrected energy, both B3LYP and M06-2X have the same order as PC2 < PC3 = PC12 < PC11 < PC1 < PC10 < PC8 < PC9 < PC15 < PC16 < PC7 < PC6 < PC13 < PC5 ~ PC4 ~ PC14. On the other hand, the relative stability order as predicted by ab initio Moller Plesset MP2 and MP4(SDQ) is PC3 = PC12 < PC2 = PC11 < PC1 < PC10 < PC8 < PC9 < PC15 < PC16 < PC6 ~ PC7 < PC13 < PC4 ~ PC14 < PC5. However, MP4(SDQ) method shows slight differences like PC15 ~ PC16 and PC4 < PC14 ~ PC5. Figure 4 presents tautomers’ order by DFT and MP2 methods. Enthalpy and free energy are consistent with the total energy as can be detected from Table 2. One of the electronic factors that affect the stability of the tautomers is the repulsion between the neighboring –OH and –PH or –CH groups, while the other factor is the attraction between the O–H or P–H bonds and the lone pair of the neighboring phosphorous or oxygen atoms. Another factor is the attractive interaction between the O–H or P–H bonds, on one side, and the lone pair(s) of the adjacent P or O atoms, on the other one .
The amino N accommodates the highest negative charge (up to −0.904); then the O atom comes in the second place (up to −0.778). Positive charges are accommodated on the P and H atoms. Very little negative charge on H bonded to P atoms as appeared in PC4, PC6, and PC8. The values of the NPA charges at hydrogen atoms of –OH groups are larger than the values at the hydrogen atom of the –CH group (see Table 3). As a result, the repulsion between the neighboring –OH and –CH groups is stronger than the repulsion of the neighboring –PH and –CH groups.
3.3. Solvation Study
One of the most crucial factors determining the tautomer distribution in the biological media is the environment. Solvent preferentially stabilizes structures with localized charges. Solvation effect has been explored in water (ε = 78.54) using the polarized continuum model (PCM) at the B3LYP/6-31+G(d) level for all the 16 tautomers. Table 4 collects the relative energies for different tautomers/rotamers (kcal/mol) at B3LYP/6-31+G(d) in water. In water, the most stable diphosphocytosine tautomers are PC3 and PC12 and the most unstable tautomers/rotamers are PC4 and PC5. As seen from Table 4, the relative stabilities in water go in the order PC3 ~ PC12 ~ PC2 ~ PC11 < PC1 < PC10 < PC8 < PC9 < PC15 ~ PC16 < PC13 < PC6 ~ PC7 ~ PC14 < PC4 ~ PC5 at B3LYP/6-31+G(d) (PCM/water). This order is slightly different from that predicted at the same level in gas phase (ε = 1). For example, PC3 ~ PC12 ~ PC2 ~ PC11. This result is easily understandable on the basis of the calculated dipole moment of these tautomers. Polar tautomers/rotamers are favored in polar solvents.
The NPA charges computed for tautomers/rotamers geometry optimized in water are given in Table 5.
The NPA charges computed for tautomers/rotamers geometry optimized in water are generally smaller than those computed for tautomers/rotamers optimized in gas phase. This is understandable as a role of the high dielectric constant solvent, water. The amino N accommodates the highest negative charge (up to −0.846), then the O atom comes in the second place (up to −0.707). Positive charges are accommodated on the P and H atoms. Very little negative charge on H bonded to P atoms as appeared in PC4, PC6, and PC8. The values of the NPA charges at hydrogen atoms of –OH groups are larger than the values at the hydrogen atom of the –CH group (see Table 5). As a result, the repulsion between the neighboring –OH and –CH groups is stronger than the repulsion of the neighboring –PH and –CH groups.
On the basis of the present theoretical study, we may confirm the following conclusions.(i)We studied for the first time 16 tautomers/rotamers of diphosphocytosine by four computational methods. Some of these tautomers/rotamers are isoenergetic although they have different structure.(ii)From the total energy point of view, the 16 tautomers/rotamers are categorized into three main groups. These are low-energy difference group, moderate-energy difference group, and high-energy difference group.(iii)PC3 and PC12 are the most stable tautomers/rotamers.(iv)Our systematic calculations show that the energy difference between components of a rotamer pair is very stable, independent of the level of theory. This leads to the important conclusion that besides the dominant PC3 form, its pair PC2 may lie not higher than ~0.4 kcal mol−1; therefore, any experimental study should consider this possibility when evaluating spectroscopic experimental studies.(v)Benzene-like structures are the most stable tautomers/rotamers. (vi)Charge is distributed on the ring atoms; −ve charges are on N and O, while +ve charges on P and H.(vii)Solvation does not affect the order significantly; PC3 ~ PC12 ~ PC2 ~ PC11 < PC1 < PC10 < PC8 < PC9 < PC15 ~ PC16 < PC13 < PC6 ~ PC7 ~ PC14 < PC4 ~ PC5.
This work was done at the Vienna Supercomputer Centre of the University of Vienna, Austria.
J. March, Advanced Organic Chemistry, Reaction, Mechanism and Structure, John Wiley & Sons, New York, NY, USA, 4th edition, 1993.
J. Elguero, C. Marzin, A. R. Katritzky, and P. Linda, The Tautomerism of Heterocycles, Advances in heterocyclic chemistry: Supplement I, Academic Press, 1976.
M. Dreyfus, O. Bensaude, G. Dodin, and J. E. Dubois, “Tautomerism in cytosine and 3-methylcytosine. A thermodynamic and kinetic study,” Journal of the American Chemical Society, vol. 98, no. 20, pp. 6338–6349, 1976.View at: Google Scholar
R. Stolarski, M. Remin, and D. Shugar, “Studies on prototropic tautomerism in neutral and monoanionic forms of pyrimidines by nuclear magnetic resonance spectroscopy,” Zeitschrift Für Naturforschung C, vol. 32, no. 5, pp. 894–900, 1977.View at: Google Scholar
J. Lin, C. Yu, S. Peng et al., “Ultraviolet photoelectron studies of the ground-state electronic structure and gas-phase tautomerism of hypoxanthine and guanine,” The Journal of Physical Chemistry, vol. 84, no. 9, pp. 1006–1012, 1980.View at: Google Scholar
M. J. Nowak, K. Szczepaniak, A. Barski, and D. Shugar, “Spectroscopic studies on vapour phase tautomerism of natural bases found in nucleic acids,” Zeitschrift fur Naturforschung, vol. 33, no. 11-12, pp. 876–883, 1978.View at: Google Scholar
P. O. LÖwdin, Quantum Genetics and Aperiodic Solids, Preprint No 32. Quantum Theory Project for Research in Atomic, Molecular and Solid state, Chemistry and Physics, University of Florida, Gainesville, Fla, USA, 1962.
J. Leszczynski, Ed., Computational Molecular Biology, vol. 8 of Theoretical and Computational Chemistry Book Series, Elsevier, 1999.
S. Kwiatkowski and J. Leszczynski, “An ab Initio quantum-mechanical study of tautomerism of purine, adenine and guanine,” Journal of Molecular Structure, vol. 208, no. 1-2, pp. 35–44, 1990.View at: Google Scholar
J. Leszczynski, “The potential energy surface of guanine is not flat: an ab Initio study with large basis sets and higher order electron correlation contributions,” Journal of Physical Chemistry A, vol. 102, no. 13, pp. 2357–2362, 1998.View at: Google Scholar
E. D. Radchenko, A. M. Plohotnichenko, G. G. Sheina, Blagoi, and P. Yu, “ab Initio study of the prototropic tautomerism of cytosine and guanine and their contribution to spontaneous point mutations,” Biophysics, vol. 28, p. 559, 1983.View at: Google Scholar
Szczepaniak and K. Szczepaniak M, “Matrix isolation infrared studies of nucleic acid constituents: Part 4. Guanine and 9-methylguanine monomers and their keto-enol tautomerism,” Journal of Molecular Structure, vol. 156, no. 1-2, pp. 29–42, 1987.View at: Google Scholar
M. J. Nowak, L. Lapinski, and J. Fulara, “Matrix isolation studies of cytosine: the separation of the infrared spectra of cytosine tautomers,” Spectrochimica Acta A, vol. 45, no. 2, pp. 229–242, 1989.View at: Google Scholar
W. Saenger, Principles of Nucleic Acid Structure, Springer, New York, NY, USA, 1994.
T. Shimanouchi, M. Tsuboi, and Y. Kyogoku, “The structure and properties of biomolecules and biological systems,” Advances in Chemical Physics, vol. 7, p. 435, 1964.View at: Google Scholar
R. C. Lord and G. J. Thomas Jr., “Raman spectral studies of nucleic acids and related molecules-I Ribonucleic acid derivatives,” Spectrochimica Acta Part A, vol. 23, no. 9, pp. 2551–2591, 1967.View at: Google Scholar
G. Lauer, W. Schäfer, and A. Schweig, “Functional subunits in the nucleic acid bases uracil and thymine,” Tetrahedron Letters, vol. 16, no. 45, pp. 3939–3942, 1975.View at: Google Scholar
D. Dougherty, K. Wittel, J. Meeks, and S. P. McGlynn, “Photoelectron spectroscopy of carbonyls. Ureas, uracils, and thymine,” Journal of the American Chemical Society, vol. 98, no. 13, pp. 3815–3820, 1976.View at: Google Scholar
A. Padva, T. G. O’Donnell, and P. R. LeBreton, “UV photoelectron studies of biological pyrimidines: the valence electronic structure of methyl substituted uracils,” Chemical Physics Letters, vol. 41, no. 2, pp. 278–282, 1976.View at: Google Scholar
M. Fujii, T. Tamura, N. Mikami, and M. Ito, “Electronic spectra of uracil in a supersonic jet,” Chemical Physics Letters, vol. 126, no. 6, pp. 583–587, 1986.View at: Google Scholar
Y. Tsuchiya, T. Tamura, M. Fujii, and M. Ito, “Keto-enol tautomer of uracil and thymine,” Journal of Physical Chemistry, vol. 92, no. 7, pp. 1760–1765, 1988.View at: Google Scholar
B. B. Brady, L. A. Peteanu, and D. H. Leavy, “The electronic spectra of the pyrimidine bases uracil and thymine in a supersonic molecular beam,” Chemical Physics Letters, vol. 147, no. 6, pp. 538–543, 1988.View at: Google Scholar
M. Kubota and T. J. Kobayashi, “Electronic structure of uracil and uridine derivatives studied by photoelectron spectroscopy,” Journal of Electron Spectroscopy and Related Phenomena, vol. 82, no. 1-2, pp. 61–70, 1996.View at: Google Scholar
M. Szczesniak, M. J. Nowak, H. Rostkowska, K. Szczepaniak, W. B. Person, and D. Shugar, “Matrix isolation studies of nucleic acid constituents,” Journal of the American Chemical Society, vol. 105, no. 19, pp. 5969–5976, 1983.View at: Google Scholar
Y. D. Radchenko, G. G. Scheina, N. A. Smorygo, Blagoi, and P. Yu, “Experimental and theoretical studies of molecular structure features of cytosine,” Journal of Molecular Structure, vol. 116, no. 3-4, pp. 387–396, 1984.View at: Google Scholar
R. D. Brown, P. D. Godfrey, D. McNaughton, and A. P. Pierlot, “Microwave spectrum of uracil,” Journal of the American Chemical Society, vol. 110, no. 7, pp. 2329–2330, 1988.View at: Google Scholar
I. Kulakowska, M. Geller, B. Lesyng, K. L. Wierzchowski, and K. Bolewska, “Barrier to rotation and conformation of the -NR2 group in cytosine and its derivatives. Part II. Experimental and theoretical dipole moments of methylated cytosines,” Biochimica Et Biophysica Acta, vol. 407, no. 4, pp. 420–429, 1975.View at: Google Scholar
P. Beak and J. M. White, “Relative enthalpies of 1, 3-dimethyl-2, 4-pyrimidinedione, 2, 4-dimethoxypyrimidine, and 4-methoxy-1-methyl-1-2-pyrimidinone: estimation of the relative stabilities of two protomers of uracil,” Journal of the American Chemical Society, vol. 104, no. 25, pp. 7073–7077, 1982.View at: Publisher Site | Google Scholar
J. S. Kwiatkowski, R. J. Bartlett, and W. B. Person, “Contributions from electron correlation to the relative stabilities of the tautomers of nucleic acid bases,” Journal of the American Chemical Society, vol. 110, no. 8, pp. 2353–2358, 1988.View at: Google Scholar
H. Basch, D. R. Garmer, P. G. Jasien, M. Krauss, and W. J. Stevens, “Electrical properties of nucleic acid bases,” Chemical Physics Letters, vol. 163, no. 6, pp. 514–522, 1989.View at: Google Scholar
A. Les, L. Adamowicz, and R. J. Bartlett, “Relative stability of cytosine tautomers with the coupled cluster method and first-order correlation orbitals,” Journal of Physical Chemistry, vol. 93, pp. 4001–4005, 1989.View at: Google Scholar
P. G. Jasien and G. Fitzgerald, “Molecular dipole moments and polarizabilities from local density functional calculations: Application to DNA base pairs,” The Journal of Chemical Physics, vol. 93, no. 4, pp. 2554–2560, 1990.View at: Google Scholar
J. Leszczynśki, “Structure and properties of uracil and its sulfur analogs: a systematic study of basis set effects in Ab Initio SCF calculations,” Physical Chemistry, vol. 40, supplement 18, pp. 9–21, 1991.View at: Google Scholar
A. R. Katritzky and M. Karelson, “Am1 calculations of reaction field effects on the tautomeric equilibria of nucleic-acid pyrimidine and purine-bases and their 1-methyl analogs,” Journal of the American Chemical Society, vol. 113, no. 5, pp. 1561–1566, 1991.View at: Google Scholar
D. A. Estrin, L. Paglieri, and G. Corongiu, “A density functional study of tautomerism of uracil and cytosine,” Journal of Physical Chemistry, vol. 98, no. 22, pp. 5653–5660, 1994.View at: Google Scholar
I. R. Gould, N. A. Burton, R. J. Hall, and I. H. Hillier, “Tautomerism in uracil, cytosine and guanine: a comparison of electron correlation predicted by ab Initio and density functional theory methods,” Journal of Molecular Structure, vol. 331, no. 1-2, pp. 147–154, 1995.View at: Google Scholar
L. Paglieri, G. Corongiu, and D. A. Estrin, “Solvent effects in density functional calculations of uracil and cytosine tautomerism,” International Journal of Quantum Chemistry, vol. 56, no. 5, pp. 615–625, 1995.View at: Google Scholar
M. Monshi, K. Al-Farhan, S. Al-Resayes, A. Ghaith, and A. A. Hasanein, “Excited states dipole moments and polarizabilities of uracil and cytosine 5-halo derivatives,” Spectrochimica Acta A, vol. 53, no. 14, pp. 2669–2677, 1997.View at: Google Scholar
S. I. Kawaharaa, T. Uchimarua, and M. Sekineb, “The hydrogen bond energy on mismatched base pair formation between uracil derivatives and guanine in the gas phase and in the aqueous phase,” Journal of Molecular Structure, vol. 530, pp. 1109–2117, 2000.View at: Google Scholar
Iwona Dabkowska, Janusz Rak, Maciej Gutowskia, J. Michael Nilles, S. T. Stokes, and K. H. Bowen, “Barrier-free intermolecular proton transfer induced by excess electron attachment to the complex of alanine with uracil,” Journal of Chemical Physics, vol. 120, no. 13, 2004.View at: Publisher Site | Google Scholar
R. Zhang, A. Ceulemans, and Nguyen, “A theoretical study of Uracil and its tautomers in their lowest-lying triplet state,” Molecular Physics, vol. 103, no. 6-8, pp. 983–994, 2005.View at: Google Scholar
J. K. Wolken and F. Turecek, “Proton affinity of uracil: a computational study of protonation sites,” Journal of the American Society For Mass Spectrometry, vol. 11, no. 12, pp. 1065–1071, 2000.View at: Google Scholar
S. D. Wetmore, R. J. Boyd, and L. A. Eriksson, “A theoretical study of 5-halouracils: electron affinities, ionization potentials and dissociation of the related anions,” Chemical Physics Letters, vol. 343, no. 1-2, pp. 151–158, 2001.View at: Google Scholar
P. U. Civcir, “AM1 and PM3 studies of some thio analogues of pyrimidine bases in the gas and aqueous phases,” Journal of Physical Organic Chemistry, vol. 14, no. 3, pp. 171–179, 2001.View at: Google Scholar
W. Saenger and D. Suck, “The Relationship between hydrogen bonding and base stacking in crystalline 4-thiouridine derivatives,” European Journal of Biochemistry, vol. 32, no. 3, pp. 473–478, 1973.View at: Google Scholar
A. G. Lezius and K. H. Scheit, “Enzymatic synthesis of DNA with 4-thio-thymidine triphosphate as substitute for dTTP,” European Journal of Biochemistry, vol. 3, no. 1, pp. 85–94, 1967.View at: Google Scholar
K. H. Scheit and E. Gartner, “Über die Eigenschaften von polynucleotiden, welche uridin und 4-thiouridin enthalten,” Biochimica Et Biophysica Acta, vol. 182, no. 1, pp. 10–16, 1969.View at: Google Scholar
E. M. Gottschalk, E. Kopp, and A. G. Lezius, “A synthetic DNA with unusual base-pairing,” European Journal of Biochemistry, vol. 24, no. 1, pp. 168–182, 1971.View at: Google Scholar
M. J. Frisch, G. W. Trucks, H. B. Schlegel et al., Gaussian 09, Revision A. 1, Gaussian Inc., Wallingford, Conn, USA, 2009.
A. D. J. Becke, “Density-functional thermochemistry. III. The role of exact exchange,” Chemical Physics, vol. 98, no. 7, pp. 5648–5652, 1993.View at: Google Scholar
C. Lee, W. Yang, and R. G. Parr, “Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density,” Physical Review B 37, pp. 785–789, 1988.View at: Google Scholar
P. J. Stephens, F. J. Devlin, C. F. Chabalowski, and M. J. Frisch, “ab Initio calculation of vibrational absorption and circular dichroism spectra using density functional force fields,” Journal of Physical Chemistry, vol. 98, no. 45, pp. 11623–11627, 1994.View at: Google Scholar
G. A. Zhurko and D. A. Zhurko, ChemCraft version 1. 5, 2005.
A. E. Reed, R. B. Weinstock, and F. Weinhold, “Natural population analysis,” The Journal of Chemical Physics, vol. 83, no. 2, pp. 735–746, 1985.View at: Google Scholar
S. Miertus, E. Scrocco, and J. Tomasi, “Electrostatic interaction of a solute with a continuum. A direct utilizaion of ab Initio molecular potentials for the prevision of solvent effects,” Journal of Chemical Physics, vol. 55, no. 1, pp. 117–129, 1981.View at: Google Scholar
J. Tomasi and M. Persico, “Molecular interactions in solution: An overview of methods based on continuous distributions of the solvent,” Chemical Reviews, vol. 94, no. 7, pp. 2027–2094, 1994.View at: Google Scholar
R. Cammi and J. Tomasi, “Remarks on the use of the apparent surface charges/ASC methods in solvation problems: iterative versus matrix-inversion procedures and the renormalization of the apparent charges,” Journal of Computational Chemistry, vol. 16, pp. 1449–1458, 1995.View at: Google Scholar
M. Cossi, V. Barone, R. Cammi, and J. Tomasi, “ab Initio study of solvated molecules: A new implementation of the polarizable continuum model,” Chemical Physics Letters, vol. 255, no. 4–6, pp. 327–335, 1996.View at: Google Scholar
J. Sponer and P. Hobza, “Nonplanar Geometries of DNA Bases. ab Initio Second-Order Moeller-Plesset Study,” Journal of Physical Chemistry, vol. 98, no. 12, pp. 3161–3164, 1994.View at: Google Scholar
E. L. Stewart, C. K. Foley, N. L. Allinger, and J. P. Bowen, “ab Initio calculations with electronic correlation (MP2) on the nucleic acid bases and their methyl derivatives,” Journal of the American Chemical Society, vol. 116, no. 16, pp. 7282–7286, 1994.View at: Google Scholar
V. B. Delchev and M. V. Nenkova, “Theoretical Modeling of the Ground State Intramolecular Proton Transfer in Cytosine: DFT Level Study,” Acta Chimica Slovenica, pp. 132–137, 2008.View at: Google Scholar
T. K. Ha, H. J. Keller, R. Gunde, and H. H. Gunthard, “Energy increment method based on quantum chemical results: a general recipe for approximative prediction of isomerization and tautomerization energies of pyrimidine and purine Nucleic acid bases and related compounds,” Journal of Physical Chemistry A, vol. 103, no. 33, pp. 6612–6623, 1999.View at: Google Scholar