Quantum chemistry of nucleic acids: how it could help and when it is necessary

Quantum chemical and spectroscopic investigations of 3-methyladenine

FTIR, FT-Raman and UV-Vis spectra of 3-methyladenine have been recorded and investigated using quantum chemical calculations. The molecular geometry and vibrational spectra of 3-methyladenine in the ground state are computed by using HF and DFT methods with 6-311G(d,p) basis set. VSCF, CC-VSCF methods based on 2MR-QFF and PT2 (Barone method) have been utilized for computing anharmonic vibrational frequencies. These methods yield results that are in remarkable agreement with the experimental data. The magnitudes of coupling between pair of modes have been also computed. Vibrational modes are assigned with the help of visual inspection of atomic displacements. The electronic spectra, simulated at TD-B3LYP/6-311++G(d,p) level of theory, are compared to the experiment. The global quantities: electronic chemical potential, electrophilicity index, chemical hardness and softness based on HOMO and LUMO energy eigenvalues are also computed at B3LYP/6-311++G(d,p) level of theory.

Charge transfer in DNA: effective Hamiltonian approaches

Aperiodic order plays a very significant role in biology, as it determines most informative content of genomes. Amongst the various physical, chemical or biological phenomena that might be inferred from sequence correlations, charge transfer properties deserve particular attention. Indeed, the nature of DNA-mediated charge migration has been related to the understanding of damage recognition process, protein binding, or with the task of engineering biological processes (e.g. designing nanoscale sensing of genomic mutations), opening new challenges for emerging nanobiotechnologies. Nevertheless, the solution of Schrödinger´s equation with a potential that is given by a one-dimensional array of the double-stranded DNA remains as a main open theme in solid state physics of biological macromolecules. In this contribution, I will shortly review several approaches introduced during the last few years in order to describe charge transfer migration in DNA in terms of tight-binding effective Hamiltonians.

Density functional study of isoguanine tetrad and pentad sandwich complexes with alkali metal ions

Isoguanine tetraplexes and pentaplexes contain two or more stacked polyads with intercalating metal ions. We report here the results of a density functional study of sandwiched isoguanine tetrad and pentad complexes consisting of two polyads with Na(+), K(+) and Rb(+) ions at the B3LYP level. In comparison to single polyad metal ion complexes, there is a trend towards increased non-planarity of the polyads in the sandwich complexes. In general, the pentad sandwiches have relatively planar polyad structures, whereas the tetrad complexes contain highly non-planar polyad building blocks. As in other sandwich complexes and in metal ion complexes with single polyads, the metal ion-base interaction energy plays an essential role. In iG sandwich structures, this interaction energy is slightly larger than in the corresponding guanine sandwich complexes. Because the base-base interaction energy is even more increased in passing from guanine to isoguanine, the isoguanine sandwiches are thus far the only examples where the base-base interaction energy is larger than the base-metal ion interaction energy. Stacking interactions have been studied in smaller models consisting of two bases, retaining the geometry from the complete complex structures. From the data obtained at the B3LYP and BH&H levels and with Møller-Plesset perturbation theory, one can conclude that the B3LYP method overestimates the repulsion in stacked base dimers. For the complexes studied in this work, this is only of minor importance because the direct inter-tetrad or inter-pentad interaction is supplemented by a strong metal ion-base interaction. Using a microsolvation model, the metal ion preference K(+) approximately Rb(+) > Na(+) is found for tetrad complexes. On the other hand, for pentads the ordering is Rb(+) > K(+) > Na(+). In the latter case experimental data are available that agree with this prediction.

Characteristic length scale of electric transport properties of genomes

A tight-binding model together with a statistical method are used to investigate the relation between the sequence-dependent electric transport properties and the sequences of protein-coding regions of complete genomes. A correlation parameter Omega is defined to analyze the relation. For some particular propagation length w max, the transport behaviors of the coding and noncoding sequences are very different and the correlation reaches its maximal value Omega max. w max and Omega max are characteristic values for each species. A possible reason for the difference between the features of transport properties in the coding and noncoding regions is the mechanism of DNA damage repair processes together with natural selection.

Localization properties of electronic states in a polaron model of poly(dG)-poly(dC) and poly(dA)-poly(dT) DNA polymers

We numerically investigate localization properties of electronic states in a static model of poly(dG)-poly(dC) and poly(dA)-poly(dT) DNA polymers with realistic parameters obtained by quantum-chemical calculation. The randomness in the on-site energies caused by the electron-phonon coupling is completely correlated to the off-diagonal parts. In the single electron model, the effect of the hydrogen-bond stretchings, the twist angles between the base pairs and the finite system size effects on the energy dependence of the localization length and on the Lyapunov exponent are given. The localization length is reduced by the influence of the fluctuations in the hydrogen bond stretchings. It is also shown that the helical twist angle affects the localization length in the poly(dG)-poly(dC) DNA polymer more strongly than in the poly(dA)-poly(dT) one. Furthermore, we show resonance structures in the energy dependence of the localization length when the system size is relatively small.

Influence of DNA stiffness in protein–DNA recognition

Protein-DNA recognition plays an essential role in the regulation of gene expression. The protein-DNA binding specificity is based on direct atomic contacts between protein and DNA and/or the conformational properties of DNA. In this work, we have analyzed the influence of DNA stiffness (E) to the specificity of protein-DNA complexes. The average DNA stiffness parameters for several protein-DNA complexes have been computed using the structure based sequence dependent stiffness scale. The relationship between DNA stiffness and experimental protein-DNA binding specificity has been brought out. We have investigated the importance of DNA stiffness with the aid of experimental free energy changes (DeltaDeltaG) due to binding in several protein-DNA complexes, such as, ETS proteins, 434, lambda, Mnt and trp repressors, 434 cro protein, EcoRV endonuclease V and zinc fingers. We found a correlation in the range 0.65-0.97 between DeltaDeltaG and E in these examples. Further, we have qualitatively analyzed the effect of mutations in the target sequence of lambda repressor and we observed that the DNA stiffness could correctly identify 70% of the correct bases among the considered nine positions.

Adiabatic Electron Affinities of the Polyhydrated Adenine−Thymine Base Pair: A Density Functional Study

Adiabatic electron affinities (AEAs) of the adenine-thymine (AT) base pair surrounded by 5 and 13 water molecules have been studied by density functional theory (DFT). Geometries of neutral AT x nH2O and anionic (AT x nH2O)- complexes (n = 5 and 13) were fully optimized, and vibrational frequency analysis was performed at the B3LYP/6-31+G** level of theory. The optimized structures of the neutral (AT x nH2O) and (AT x nH2O)- complexes were found to be somewhat nonplanar. Some of the water molecules are displaced away from the AT ring plane and linked with one another by hydrogen bonds. The optimized structures of the complexes are found to be in a satisfactory agreement with the observed experimental and molecular dynamics simulation results. In the optimized anionic complexes, the thymine (T) moiety was found to be puckered, whereas the adenine (A) moiety remained almost planar. Natural population analysis (NPA) performed using the B3LYP/6-31+G** method shows that the thymine moiety in the anionic (AT x nH2O)- complexes (n = 5 and 13) has most of the excess electronic charge, i.e., approximately -0.87 and approximately -0.81 (in the unit of magnitude of the electronic charge), respectively. The zero-point energy corrected adiabatic electron affinities of the hydrated AT base pair were found to be positive both for n = 5 and 13 and have the values of 0.97 and 0.92 eV, respectively, which are almost three times the AEA of the AT base pair. The results show that the presence of water molecules appreciably enhances the EA of the base pair.

A theoretical study of structures and electron affinities of radical anions of guanine-cytosine, adenine-thymine, and hypoxanthine-cytosine base pairs

Adiabatic electron affinities (AEA) and structural perturbations due to addition of an excess electron to each of the neutral guanine-cytosine (G-C), adenine-thymine (A-T), and hypoxanthine-cytosine (HX-C) base pairs were studied using the self-consistent charge, density functional tight-binding (SCC-DFTB-D) method, augmented by the empirical London dispersion energy term. Performance of the SCC-DFTB-D method was examined by comparing the calculated results using it with those obtained from experiment as well as ab initio and other different density functional theoretical studies. An excellent agreement between the SCC-DFTB-D results and those obtained by the other calculations regarding the structural modifications, hydrogen bonding, and dissociation energies of the neutral and radical anion base pairs was found. It is shown that adiabatic electron affinity can be better predicted by considering reaction enthalpies of formation of the respective neutral and anionic base pairs from their respective molecular components instead of taking the difference between their total energies. The calculated AEAs of the base pairs were compared with those obtained by the bracketing method from Schaefer and coworkers, where a satisfactory agreement was found. It shows applicability of the SCC-DFTB-D method to study charged DNA models at a highly economical computational cost.