Accuracy and efficiency of atomic basis set methods versus plane wave calculations with ultrasoft pseudopotentials for DNA base molecules

Computational chemistry applied to vibrational spectroscopy: A tool for characterization of nucleic acid bases and some of their 5-substituted derivatives

Computational chemistry can be applied to vibrational spectroscopy in different ways, such as for a better characterization and assignment of all the bands of the experimental spectra, as a tool in the identification of the tautomers present in the gas phase and in the solid state through their spectra and for the simulation of the solid and liquid phase of a compound and the consequent simulation and interpretation of their spectra. In the present study, as an example of the applicability of computational chemistry, the structure and spectra of cytosine and uracil nucleic acid bases and two cytosine derivatives are shown. The FTIR and Raman spectra were analysed with the support of ab initio (Hartree-Fock (HF), MP2) and density functional theory (DFT) (B3LYP, PBE, B-P, etc.) calculations using several basis sets and several scaling equations. The calculations predict an easier tautomerization of cytosine than uracil molecule, but the tautomerization is hindered in the 5-bromocytosine molecule. Thus, in the solid state, this molecule only exists in the amino-oxo tautomeric form.

The biomolecule of 5-bromocytosine: FT-IR and FT-Raman spectra and DFT calculations. Identification of the tautomers in the isolated state and simulation the spectra in the solid state

An accurate assignment of the IR spectrum in Ar matrix of 5-bromocytosine and of the IR and Raman spectra in the solid state was carried out. For this purpose Density functional calculations (DFTs) were performed to clarify wavenumber assignments of the experimental observed bands. The calculated values were scaled using scaling equations and they were compared with IR and Raman experimental data. Good reproduction of the experimental wavenumbers is obtained and the% error is very small in the majority of cases. In the isolated state all the tautomer forms of 5-bromocytosine were determined and optimized. The wavenumbers corresponding to C1 and C2b tautomers were identified and assigned in the IR experimental spectrum reported in Ar matrix. Our study confirms the existence of at least two tautomers, the amino-oxo and the amino-hydroxy in the isolated state. In the solid state the FT-IR and FT-Raman spectra of 5-bromocytosine in the powder form were recorded in the region 400-4000 cm(-1) and 50-3500 cm(-1), respectively. The unit cell found in the crystal was simulated as a tetramer form in three tautomers. Thus, it has been possible to assign all the 33 normal modes of vibration. The study indicates that the features, that are the characteristic of the vibrational spectra of cytosine, are retained by the spectra of 5-bromocytosine and it exists in the solid phase in the amino-oxo form.

A comparison between plane wave and Gaussian-type orbital basis sets for hydrogen bonded systems: formic acid as a test case

The formic acid molecule, its dimers, and its molecular crystal are adopted as test systems to compare results obtained with plane wave (PW) basis sets and norm-conserving pseudopotentials to all-electron Gaussian-type orbital (GTO) calculations. The CPMD and CRYSTAL06 codes, respectively, are applied with the PBE, PW91, and BLYP density functionals. Hydrogen bonding is the leading interaction in the dimers and the crystal. In the latter, dispersive and weak C-H...O interactions are also relevant. Irrespective of the adopted functional, for all considered structures PW and GTO results converge smoothly as a function of the quality of the adopted basis sets to the same values for structures, energies of interaction, and harmonic vibrational features. To achieve a high level of mutual agreement the use of GTO basis sets of at least of triple-zeta quality including one set of polarization functions and PW basis sets with a kinetic energy cutoff higher than 110 Ry is recommended. Pros and cons of both approaches for studying molecular crystals are also discussed.

The small planarization barriers for the amino group in the nucleic acid bases

The amino group in the nucleic acid bases frequently interacts with other bases or with other molecular systems. Thus any nonplanarity of the amino group may affect the molecular recognition of nucleic acids. Ab initio Hartree-Fock (HF) and second-order Moller-Plesset perturbation (MP2) levels of theory have been used to obtain the equilibrium geometries of the C(l) and C(s) structures for five common nucleic acid bases. The energy barriers between the C(l) and C(s) structures have also been predicted. A series of correlation consistent basis sets up to cc-pCVQZ and aug-cc-pVQZ has been used to systematically study the dependence of the amino group nonplanarity. The equilibrium geometries of the nucleic acid bases with an amino group, including adenine, guanine, and cytosine, are examined carefully. At the MP2 level of theory, larger basis sets decrease the extent of nonplanarity of the amino group, but the decrease slows down when the QZ basis sets are used, demonstrating the intrinsic property of nonplanarity for guanine. For adenine and cytosine the situation is less clear; as the HF limit is approached, these two structures become planar. Addition of core correlation effects or diffuse functions further decreases the degree of nucleic acid base nonplanarity, in comparison to the original cc-pVXZ (X=D, T, and Q) basis sets. The aug-cc-pVXZ basis shows smaller degrees of nonplanarity than the cc-pCVXZ sets. The aug-cc-pVXZ basis is less size dependent than the cc-pVXZ and cc-pCVXZ sets in the prediction of the amino-group-related bond angles and dihedral angles and energy barriers for adenine, guanine, and cytosine. The cc-pCVQZ and aug-cc-pVQZ MP2 results may be regarded as benchmark predictions for the five common bases. The predicted classical barriers to planarization are 0.02 (adenine), 0.74 (guanine), and 0.03(cytosine) kcal mol(-1).

Coulombic amino group-metal bonding: adsorption of adenine on Cu110

The interaction between molecular amino groups and metal surfaces is analyzed from first-principles calculations using the adsorption of adenine on Cu110 as a model case. The amino group nitrogens are found to adsorb on top of the surface copper atoms. However, the bonding clearly cannot be explained in terms of covalent interactions. Instead, we find it to be largely determined by mutual polarization and Coulomb interaction between substrate and adsorbate.