Prebiotic chemistry presumably took place before formation of an oxygen-rich atmosphere and thus under conditions of intense short wavelength UV irradiation. Therefore, the UV photochemical stability of the molecular building blocks of life may have been an important selective factor in determining the eventual chemical makeup of critical biomolecules. To investigate the role of UV irradiation in base-pairing we have studied guanine (G) and cytosine (C) base pairs in the absence of the RNA backbone. We distinguished base-pair structures by IR-UV hole-burning spectroscopy as well as by high-level correlated ab initio calculations. The Watson-Crick structure exhibits broad UV absorption, in stark contrast to other GC structures and other base-pair structures. This broad absorption may be explained by a rapid internal conversion that makes this specific base pair arrangement uniquely photochemically stable.ab initio computation ͉ DNA base pairs ͉ IR-UV spectroscopy ͉ jet cooling ͉ photochemistry T he DNA bases involved in reproduction have short S 1 excited state lifetimes, of the order of 1 ps or less (1-7). It has been argued that this phenomenon serves to protect these bases against photochemical damage, because after excitation they do not cross to the reactive triplet state; instead, they rapidly internally convert to the electronic ground state (8). This mechanism may have been particularly significant under the conditions of the early earth, when purines and pyrimidines presumably were assembled into the first macromolecular structures, producing RNA. At that time, the earth was exposed to shorter wavelength UV radiation than it is today. For an analysis of possible prebiotic chemistry, it is necessary not only to consider the individual bases but also to study them as they interact and to do so without the RNA backbone (9). We achieve such an analysis by studying clusters of guanine (G) and cytosine (C) in the gas phase by using double resonant laser spectroscopy techniques. The experimental studies are accompanied by stateof-the-art quantum chemical and molecular dynamics calculations of the pairing of G and C. These two bases can form many different hydrogen-bonded structures, of which at least 20 are within 12 kcal͞mol of the global minimum. Thus, the familiar Watson-Crick (WC) structure may not be unique, based on energies alone. Here, we report a remarkable difference in excited state properties among the WC structure and other structures. The former exhibits broad UV absorption features; this is in contrast with the sharp UV spectra exhibited by non-WC structures. The broad absorption can be explained by recent theoretical results, predicting a pathway for rapid internal conversion (10). If the difference is solely due to lifetime broadening, these results correspond to a lifetime for the WC structures that is at least 2 orders of magnitude shorter than those of observed non-WC structures. Thus, it appears possible that the WC recognition mechanism involves a structure that was significantly more stab...

Altogether eight keto and enol tautomers of guanine were studied theoretically in the gas phase, in a microhydrated environment (1 and 2 water molecules) and in bulk water. The structures of isolated, as well as mono- and dihydrated tautomers were determined by means of the RI-MP2 method using the extended TZVPP (5s3p2d1f/3s2p1d) basis set. The relative energies of isolated tautomers included the correction to higher correlation energy terms evaluated at the CCSD(T)/aug-cc-pVDZ level. The relative enthalpies at 0 K and relative free energies at 298 K were based on the above-mentioned relative energies and zero-point vibration energies, temperature-dependent enthalpy terms and entropies evaluated at the MP2/6-31G level. The keto form having hydrogen atom at N7 is the global minimum while the canonical form having hydrogen atom at N9 represents the first local minimum at all theoretical levels in vacuo and in the presence of 1 and 2 water molecules. All three unusual rare tautomers having hydrogens at N3 and N7, at N3 and N9, and also at N9 and N7 are systematically considerably less stable and can be hardly detected in the gas phase. The theoretical predictions fully agree with existing theoretical as well as experimental results. The effect of bulk solvent on the relative stability of guanine tautomers was studied by self-consistent reaction field and molecular dynamics free energy calculations using the thermodynamic integration method. Bulk solvent, surprisingly, strongly favored these three rare tautomers over all remaining low-energy tautomers and probably only these forms can exist in water phase. The global minimum (tautomer with hydrogens at N3 and N7) is by 13 kcal/mol more stable than the canonical form (3rd local minimum). Addition of one or two water molecules does not change the relative stability order of isolated guanine tautomers but the respective trend clearly supports the surprising stabilization of three rare forms.

Altogether, 14 amino and imino tautomers of adenine were studied theoretically in the gas phase, in a microhydrated environment (one and two water molecules), and in bulk water environment using the thermodynamic integration method (MD-TI), conductor-like polarizable continuum model (C-PCM, COSMO), and a hybrid model (C-PCM + one to three explicit water molecules). The structures and relative energies of various tautomers were determined at the RI-MP2 level using the TZVPP basis set. The relative enthalpies at 0 K and relative free energies at 298 K were based on relative energies and zero-point vibration energies, temperature-dependent enthalpy terms, and entropies evaluated at the MP2/6-31G** level. The effect of bulk solvent on the relative stability of adenine tautomers was studied by molecular dynamics free energy calculations using the thermodynamic integration method and self-consistent reaction field. The dipole moment of the canonical form is rather small (2.8 D) but three rare imino tautomers have very large dipole moments (more than 10 D). The canonical form is the global minimum at all theoretical levels in the gas phase, in a microhydrated environment, and in the bulk water. Two unusual rare amino tautomers having hydrogens at N3 and N7, respectively, are less stable in the gas phase by more than 7 kcal/mol and represent the first and the second local minimum. Microhydration, as well as bulk water, stabilizes these unusual tautomers, and the energy gap between them and the canonical form is reduced, but the canonical tautomer remains the global minimum in all three phases. Relative free energies (T = 298 K) of these two unusual tautomers in the bulk water evaluated by molecular dynamics free energy calculations are 2.5 and 2.8 kcal/mol, which supports their coexistence in this phase. The C-PCM results agree well with the MD-TI data, and the agreement became close when considering not only the bare tautomers but their complexes with several water molecules representing first solvation shell. Other tautomers are considerably less stable (by 12−45 kcal/mol), and neither a microhydrated environment nor bulk water can change this unfavorable tautomeric equilibrium. The theoretical data predicting the coexistence of the canonical form and the N3 and the N7 tautomers in bulk water nicely agreed with experimental data obtained from NMR measurements of the adenine tautomers in DMSO (Laxer, A.; Major, D. T.; Gottlieb, H. E.; Fischer, B. J. Org. Chem. 2001, 66, 5463.)

Properties of isolated intercalators (ethidium (E), daunomycin (D), ellipticine (EL), and 4,6'-diaminide-2-phenylindole (DAPI)) and their stacking interactions with adenine...thymine (AT) and guanine...cytosine (GC) nucleic acid base pairs were investigated by means of a nonempirical correlated ab initio method. All intercalators exhibit large charge delocalization, and none of them (including the DAPI dication) exhibits a site with dominant charge. All intercalators have large polarizability and are good electron acceptors, while base pairs are good electron donors. MP2/6-31G*(0.25) stabilization energies of intercalator...base pair complexes are large (E...AT, 22.4 kcal/mol; D...GC, 17.8 kcal/mol; EL...GC, 18.2 kcal/mol; DAPI...GC, 21.1 kcal/mol) and are well reproduced by modified AMBER potential (van der Waals radii of intercalator atoms are enlarged and their energy depths are increased). Standard AMBER potential underestimates binding, especially for DAPI-containing complexes. Because the DAPI dication is the best electron acceptor (among all intercalators studied), this difference is explained by the importance of the charge-transfer term, which is not included in the AMBER potential. For the neutral EL molecule, the standard AMBER force field provides correct results. The Hartree-Fock and DFT/B3LYP methods, not covering the dispersion energy, fail completely to reveal any energy minimum at the potential energy curve of the E...AT complex, and these methods thus cannot be recommended for a study of intercalation process. On the other hand, an approximate version of the DFT method, which was extended to cover London dispersion energy, yields for all complexes very good stabilization energies that are well comparable with referenced ab initio data. Besides the vertical dependence of the interaction, an energy twist dependence of the interaction energy was also investigated by a reference correlated ab initio method and empirical potentials. It is concluded that, despite the cationic (E +1, D +1, DAPI +2) or polar (EL) character of the intercalators investigated, it is the dispersion energy which predominantly contributes to the stability of intercalator...base pair complexes. Any procedure which does not cover dispersion energy is thus not suitable for studying the process of intercalation.

Altogether 13 keto and enol tautomers of uracil and 13 keto and enol tautomers of thymine were studied theoretically in the gas phase, in a microhydrated environment (1 and 2 water molecules) and in a water environment. Bulk water was described using the thermodynamic integration method, Conductor-like polarizable continuum model (C-PCM, COSMO) and hybrid model (C-PCM + 1-2 explicit water molecules). The structures of various tautomers were determined at the RI-MP2 level using the TZVPP basis set while relative energies were determined at the CCSD(T) level. The relative free energies at 298 K were based on the relative energies mentioned above and zero-point vibration energies, and temperature dependent enthalpy terms and entropies evaluated at the MP2/6-31G** level. The effect of bulk solvent on the relative stability of uracil and thymine tautomers was studied using molecular dynamics free energy calculations by means of the thermodynamic integration method and self-consistent reaction field. Despite the completely different nature of these methods they provide comparable solvation free energies. Besides theoretical investigation, experimental detection of uracil and thymine tautomers was performed by means of steady-state fluorescence. We conclude that it is impossible to utilize the method used by Suwaiyan and Morsy (M. A. Morsy, A. M. Al-Somali and A. Suwaiyan, J. Phys. Chem. B, 1999, 103(50), 11205) for tautomer detection, even if a very sensitive fluorimeter is used. Theoretical relative energies and free energies for isolated uracil and thymine tautomers support the existence of the canonical form only. The microhydrated environment and bulk solvent stabilize enol forms more than the canonical keto one, but gas phase destabilization of these enol forms is too high. Population of rare enol forms of uracil and thymine in bulk water will thus be very low and canonical structure will also be dominant in this phase.

No abstract

No abstract

The results of harmonic and anharmonic frequency calculations on a guanine-cytosine complex with an enolic structure (a tautomeric form with cytosine in the enol form and with a hydrogen at the 7-position on guanine) are presented and compared to gas-phase IR-UV double resonance spectral data. Harmonic frequencies were obtained at the RI-MP2/cc-pVDZ, RI-MP2/TZVPP, and semiempirical PM3 levels of electronic structure theory. Anharmonic frequencies were obtained by the CC-VSCF method with improved PM3 potential surfaces; the improved PM3 potential surfaces are obtained from standard PM3 theory by coordinate scaling such that the improved PM3 harmonic frequencies are the same as those computed at the RI-MP2/cc-pVDZ level. Comparison of the data with experimental results indicates that the average absolute percentage deviation for the methods is 2.6% for harmonic RI-MP2/cc-pVDZ (3.0% with the inclusion of a 0.956 scaling factor that compensates for anharmonicity), 2.5% for harmonic RI-MP2/TZVPP (2.9% with a 0.956 anharmonicity factor included), and 2.3% for adapted PM3 CC-VSCF; the empirical scaling factor for the ab initio harmonic calculations improves the stretching frequencies but decreases the accuracy of the other mode frequencies. The agreement with experiment supports the adequacy of the improved PM3 potentials for describing the anharmonic force field of the G...C base pair in the spectroscopically probed region. These results may be useful for the prediction of the pathways of vibrational energy flow upon excitation of this system. The anharmonic calculations indicate that anharmonicity along single mode coordinates can be significant for simple stretching modes. For several other cases, coupling between different vibrational modes provides the main contribution to anharmonicity. Examples of strongly anharmonically coupled modes are the symmetric stretch and group torsion of the hydrogen-bonded NH2 group on guanine, the OH stretch and torsion of the enol group on cytosine, and the NH stretch and NH out-of-plane bend of the non-hydrogen-bonded NH group on guanine.