Nanoswitches based on DNA base pairs: why adenine-thymine is less suitable than guanine-cytosine

Effect of the Solvent and Substituent on Tautomeric Preferences of Amine-Adenine Tautomers

Adenine is one of the basic molecules of life; it is also an important building block in the synthesis of new pharmaceuticals, electrochemical (bio)sensors, or self-assembling molecular materials. Therefore, it is important to know the effects of the solvent and substituent on the electronic structure of adenine tautomers and their stability. The four most stable adenine amino tautomers (9H, 7H, 3H, and 1H), modified by substitution (C2- or C8-) of electron-withdrawing NO2 and electron-donating NH2 groups, are studied theoretically in the gas phase and in solvents of different polarities (1 ≤ ε < 109). Solvents have been modeled using the polarizable continuum model. Comparison of the stability of substituted adenine tautomers in various solvents shows that substitution can change tautomeric preferences with respect to the unsubstituted adenine. Moreover, C8 substitution results in slight energy differences between tautomers in polar solvents (<1 kcal/mol), which suggests that in aqueous solution, C8-X-substituted adenine systems may consist of a considerable amount of two tautomers-9H and 7H for X = NH2 and 3H and 9H for X = NO2. Furthermore, solvation enhances the effect of the nitro group; however, the enhancement strongly depends on the proximity effects. This enhancement for the NO2 group with two repulsive N···ON contacts can be threefold higher than that for the NO2 with one attractive NH···ON contact. The proximity effects are even more significant for the NH2 group, as the solvation may increase or decrease its electron-donating ability, depending on the type of proximity.

Effect of Alkali Metal Cations on Length and Strength of Hydrogen Bonds in DNA Base Pairs

For many years, non-covalently bonded complexes of nucleobases have attracted considerable interest. However, there is a lack of information about the nature of hydrogen bonding between nucleobases when the bonding is affected by metal coordination to one of the nucleobases, and how the individual hydrogen bonds and aromaticity of nucleobases respond to the presence of the metal cation. Here we report a DFT computational study of nucleobase pairs interacting with alkali metal cations. The metal cations contribute to the stabilization of the base pairs to varying degrees depending on their position. The energy decomposition analysis revealed that the nature of bonding between nucleobases does not change much upon metal coordination. The effect of the cations on individual hydrogen bonds were described by changes in VDD charges on frontier atoms, H-bond length, bond energy from NBO analysis, and the delocalization index from QTAIM calculations. The aromaticity changes were determined by a HOMA index.

Mutual Relations between Substituent Effect, Hydrogen Bonding, and Aromaticity in Adenine-Uracil and Adenine-Adenine Base Pairs

The electronic structure of substituted molecules is governed, to a significant extent, by the substituent effect (SE). In this paper, SEs in selected nucleic acid base pairs (Watson-Crick, Hoogsteen, adenine-adenine) are analyzed, with special emphasis on their influence on intramolecular interactions, aromaticity, and base pair hydrogen bonding. Quantum chemistry methods-DFT calculations, the natural bond orbital (NBO) approach, the Harmonic Oscillator Model of Aromaticity (HOMA) index, the charge of the substituent active region (cSAR) model, and the quantum theory of atoms in molecules (QTAIM)-are used to compare SEs acting on adenine moiety and H-bonds from various substitution positions. Comparisons of classical SEs in adenine with those observed in para- and meta-substituted benzenes allow for the better interpretation of the obtained results. Hydrogen bond stability and its other characteristics (e.g., covalency) can be significantly changed as a result of the SE, and its consequences are dependent on the substitution position. These changes allow us to investigate specific relations between H-bond parameters, leading to conclusions concerning the nature of hydrogen bonding in adenine dimers-e.g., H-bonds formed by five-membered ring nitrogen acceptor atoms have an inferior, less pronounced covalent nature as compared to those formed by six-membered ring nitrogen. The energies of individual H-bonds (obtained by the NBO method) are analyzed and compared to those predicted by the Espinosa-Molins-Lecomte (EML) model. Moreover, both SE and H-bonds can significantly affect the aromaticity of adenine rings; long-distance SEs on π-electron delocalization are also documented.

Substituted adenine quartets: interplay between substituent effect, hydrogen bonding, and aromaticity

Adenine, one of the components of DNA/RNA helices, has the ability to form self-organizing structures with cyclic hydrogen bonds (A₄), similar to guanine quartets. Here, we report a computational investigation of the effect of substituents (X = NO₂, Cl, F, H, Me, and NH₂) on the electronic structure of 9H-adenine and its quartets (A₄-N1, A₄-N3, and A₄-N7). DFT calculations were used to show the relationships between the electronic nature of the substituents, strength of H-bonds in the quartets, and aromaticity of five- and six-membered rings of adenine. We demonstrated how the remote substituent X modifies the proton-donating properties of the NH₂ group involved in the H-bonds within quartets and how the position of the substituent and its electronic nature affect the stability of the quartets. We also showed the possible changes in electronic properties of the substituent and aromaticity of adenine rings caused by tetramer formation. The results indicate that the observed relationships depend on the A₄ type. Moreover, the same substituent can both strengthen and weaken intermolecular interactions, depending on the substitution position.

Substituent effects on hydrogen bonds in adenine quartets and aromaticity of adenine rings depend on the quartet type. A₄-N3 and A₄-N7 quartets are more responsive to the electronic nature of substituents than A₄-N1.

Substituent effects on the stability of the four most stable tautomers of adenine and purine

Substituent effects at the C2-, C8- and N-positions of adenine and purine in their four the most stable tautomers are studied by means of B97D3/aug-cc-pvdz computation applying substituents of varying electronic properties: NO₂, CN, CHO, Cl, F, H, Me, OMe, OH and NH₂. The substituent effect is characterized by the substituent effect stabilization energy (SESE) and substituent Hammett constant σ. For adenine, SESE is obtained with purine as the reference system. Additionally, for both adenine and purine, SESE characteristics are estimated with benzene, imidazole and amino-pyrimidine as reference systems, when possible, taking into account substitution in topologically equivalent positions. The role of a C6-NH₂ group in adenine in modifying the substitution effect is observed and discussed. Additionally, the proximity effect for some asymmetric substituents (e.g. CHO, OMe) is recognized and meticulously analyzed.

The Nature of Hydrogen Bonds: A Delineation of the Role of Different Energy Components on Hydrogen Bond Strengths and Lengths

Hydrogen bonds are a complex interplay between different energy components, and their nature is still subject of an ongoing debate. In this minireview, we therefore provide an overview of the different perspectives on hydrogen bonding. This will be done by discussing the following individual energy components: 1) electrostatic interactions, 2) charge-transfer interactions, 3) π-resonance assistance, 4) steric repulsion, 5) cooperative effects, 6) dispersion interactions and 7) secondary electrostatic interactions. We demonstrate how these energetic factors are essential in a correct description of the hydrogen bond, and discuss several examples of systems whose energetic and geometrical features are not captured by easy-to-use predictive models.

The influence of substituents and the environment on the NMR shielding constants of supramolecular complexes based on A-T and A-U base pairs

In the present study, we have theoretically analyzed supramolecular complexes based on the Watson-Crick A-T and A-U base pairs using dispersion-corrected density functional theory (DFT). Hydrogen atoms H8 and/or H6 in the natural adenine and thymine/uracil bases were replaced, respectively, by substituents X8, Y6 = NH^-, NH₂, NH₃⁺ (N series), O^-, OH, OH₂⁺ (O series), F, Cl or Br (halogen series). We examined the effect of the substituents on the hydrogen-bond lengths, strength and bonding mechanism, and the NMR shielding constants of the C2-adenine and C2-thymine/uracil atoms in the base pairs. The general belief in the literature that there is a direct connection between changes in the hydrogen-bond strength and the C2-adenine shielding constant is conclusively rejected by our computations.

Solvation effects on the static and dynamic first-order electronic and vibrational hyperpolarizabilities of uracil: a polarized continuum model investigation

Electronic (β(e)) and vibrational (β(v)) first-order hyperpolarizabilities of uracil were determined in gas and water solution using the Coulomb-attenuating Density Functional Theory level with the Dunning's correlation-consistent aug-cc-pVDZ basis set. Frequency-dependent β(e) values were computed for the Second Harmonic Generation (SHG) and Electric Optical Pockels Effect (EOPE) nonlinear optical phenomena. The Polarized Continuum Model was employed to study the solvent effects on the electronic and vibrational properties. The introduction of solvation contributions increases the β(e) (static) value by ca. 110%. In comparison, smaller enhancements are found for the β(e) (EOPE) and β(e) (SHG) data evaluated at the typical wavelength of 694 nm (by 40-50%). The gas-water hyperpolarizability difference was rationalised through a density analysis study. The magnitudes of the vibrational first-order hyperpolarizabilities are comparable to their electronic counterparts and noticeably increase in solution: β(v) (EOPE) ~ β(e) (EOPE) in aqueous phase at λ = 694 nm. Analysis of the IR and Raman spectra is useful to elucidate the most important contributing modes to the vibrational first-order hyperpolarizabilities.

Electron correlation effects and density analysis of the first-order hyperpolarizability of neutral guanine tautomers

Dipole moments (μ), charge distributions, and static electronic first-order hyperpolarizabilities (β(μ)) of the two lowest-energy keto tautomers of guanine (7H and 9H) were determined in the gas phase using Hartree-Fock, Møller-Plesset perturbation theory (MP2 and MP4), and DFT (PBE1PBE, B97-1, B3LYP, CAM-B3LYP) methods with Dunning's correlation-consistent aug-cc-pVDZ and d-aug-cc-pVDZ basis sets. The most stable isomer 7H exhibits a μ value smaller than that of the 9H form by a factor of ca. 3.5. The β μ value of the 9H tautomer is strongly dependent on the computational method employed, as it dramatically influences the β(μ) (9H)/β(μ) (7H) ratio, which at the highest correlated MP4/aug-cc-pVDZ level is predicted to be ca. 5. The Coulomb-attenuating hybrid exchange-correlation CAM-B3LYP method is superior to the conventional PBE1PBE, B3LYP, and B97-1 functionals in predicting the β(μ) values. Differences between the largest diagonal hyperpolarizability components were clarified through hyperpolarizability density analyses. Dipole moment and first-order hyperpolarizability are molecular properties that are potentially useful for distinguishing the 7H from the 9H tautomer.

Remote communication in a DNA-based nanoswitch

Nano-telegraph: It is shown quantum chemically how the hydrogen-bond strength of a DNA base pair can be remotely switched (across several nanometers) between three states (weak, intermediate, and strong) through a π-conjugated "organic wire".