Adenine Dimers Research Articles

Intermolecular hydrogen bonding delineates the stability of non-canonical adenine base pairs: a first-principles study.

Non-canonical nucleobase pairs differ from canonical Watson-Crick (WC) pairs in their hydrogen bonding patterns. This study uses density functional theory with empirical dispersion correction to examine the stability and electronic properties of free adenine dimers stabilized by hydrogen bonds along the WC, Sugar (S), and Hoogsteen (H) edges. Dispersion correction is crucial for accurate interaction energy evaluation. The most stable adenine dimer is stabilized by N-H⋯N hydrogen bonds in gas and solvent phases. Binding energy decreases by ∼10.2 kcal mol-1 for dimers with both C-H⋯N and N-H⋯N bonds, increasing the donor-acceptor distance. However, with a sugar-phosphate backbone, dimers with C-H⋯N and N-H⋯N bonds have higher binding energy in an implicit solvent, emphasizing the role of C-H⋯N interactions in stability and nucleic acid folding dynamics. This study highlights noncovalent interactions, such as hydrogen bonding and π-π stacking, within adenine pairs with potential applications in biosensing and DNA-based self-assembly on nanomaterial interfaces.

Influence of the Number, Nature and Position of Methyl Posttranscriptional Modifications on Nucleobase Stacking in RNA.

DFT calculations are employed to quantify the influence of the presence, number, nature, and position of posttranscriptional methylation on stacking strength of RNA bases. We carry out detailed potential energy scans of the variation in stacking energies with characteristic geometrical parameters in three categories of forty stacked dimers - canonical base homodimers (N||N), methylated base homodimers (mN||mN) and heterodimers of canonical bases and methylated counterparts (N||mN). Our analysis reveals that neutral methylation invariably enhances the stacking of bases. Further, N||mN stacking is stronger than mN||mN stacking and charged N||mN exhibit strongest stacking among all dimers. This indicates that methylations greatly enhance stacking when dispersed in RNA sequences containing identical bases. Comparison of stacks involving singly- and doubly-methylated purines reveal that incremental methylation enhances the stacking in neutral dimers. Although methylation at the carbon position of neutral pyrimidine dimers greatly enhances the stacking, methylation on the 5-membered ring imparts better stacking compared to methylation on the 6-membered ring in adenine dimers. However, methylation at the ring nitrogen (N1 ) provides better stacking than the amino group (N2 ) in guanine dimers. Our results thus highlight subtle structural effects of methylation on RNA base stacking and will enhance our understanding of the physicochemical principles of RNA structure and dynamics.

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.

Stacking of nucleic acid bases: optimization of the computational approach\u2014the case of adenine dimers

Stacking interactions play an important role in stabilizing DNA and RNA secondary structure. To select a computational level to study the stacking interactions, both energy and geometric criteria, as well as the time necessary to optimize the system, should be taken into account. In this work, an attempt was made to find the most optimal level of theory describing the stacking interactions in adenine dimers. The obtained results have shown that for this purpose, wB97XD/6-311G(p,d), wB97XD/aug-cc-pvdz, or B97D3/aug-cc-pvdz should be used. What is more, geometry of the most preferable arrangements of molecules was also pointed out, ensuring an optimal starting system for further analyses.

Calcium Induces Mitochondrial Oxidative Stress Because of its Binding to Adenine Nucleotide Translocase

Several studies have demonstrated that the mitochondrial membrane switches from selective to non-selective permeability because of its improved matrix Ca2+ accumulation and oxidative stress. This process, known as permeability transition, evokes severe dysfunction in mitochondria through the opening of a non-specific pore, whose chemical nature is still under discussion. There are some proposals regarding the components of the pore structure, e.g., the adenine nucleotide translocase and dimers of the F1 Fo-ATP synthase. Our results reveal that Ca2+ induces oxidative stress, which not only increases lipid peroxidation and ROS generation but also brings about both the collapse of the transmembrane potential and the membrane release of cytochrome c. Additionally, it is shown that Ca2+ increases the binding of the probe eosin-5-maleimide to adenine nucleotide translocase. Interestingly, these effects are diminished after the addition of ADP. It is suggested that pore opening is caused by the binding of Ca2+ to the adenine nucleotide translocase.

Symmetry-adapted perturbation theory potential for the adenine dimer.

A new intermolecular interaction potential for the adenine dimer has been developed with the help of a combination of symmetry-adapted perturbation theory and density functional theory (DFT-SAPT). Supermolecular intermolecular interaction energy computations on hydrogen-bonded and stacked adenine dimers at B3LYP-D, MP2, SCS-MP2, SCS-MI-MP2 and CCSD(T) levels showed that DFT-SAPT is in a very good agreement with CCSD(T). The developed ab initio intermolecular potential has been used to predict the cluster structure of adenine dimers, trimers and tetramers. These global cluster optimizations reproduced adenine dimers reported in the literature and moreover new low-energy structures were also located. For trimers and tetramers, new hydrogen-bonded and stacked low-energy structures have also been found. The current findings suggest that the new ab initio potential can further be exploited to reveal the structure and energy of much larger supramolecular adenine clusters.

Excited State Pathways Leading to Formation of Adenine Dimers

The reaction intermediate in the path leading to UV-induced formation of adenine dimers A═A and AA* is identified for the first time quantum mechanically, using PCM/TD-DFT calculations on (dA)2 (dA: 2'deoxyadenosine). In parallel, its fingerprint is detected in the absorption spectra recorded on the millisecond time-scale for the single strand (dA)20 (dA: 2'deoxyadenosine).

Thermally Controlled Adenine Dimer Chain Rotation on Cu(110): The Critical Role of van der Waals Interactions

The theory of adenine dimer chain assembly on Cu(110) surface is controversial, due in large to the lack of an adequate description of the van der Waals (vdW) interactions. Here, we show by a combined first-principles calculation and experiment that the role of the vdW interactions is to tilt the energy balance at different levels of chain hierarchy. We find that the stable chains are made of metastable dimers, whereas the metastable chains, close in energy to the stable ones, are made of stable dimers. As such, at room temperate deposition, adenine dimers exist primarily in their stable form. This leads to the formation of metastable chains. By annealing at elevated temperature, however, more dimers can exist in the metastable form. This leads to the nucleation of the stable chains at a different orientation. The thermally controlled chain rotation is expected to be of general importance to the understanding of amino acids assembly and functionalization at the most elemental level.

Dispersion corrected DFT approaches for anharmonic vibrational frequency calculations: nucleobases and their dimers.

Computational spectroscopy techniques have become in the last few years an effective means to analyze and assign infrared (IR) spectra of molecular systems of increasing dimensions and in different environments. However, transition from compilation of harmonic data to fully anharmonic simulations of spectra is still underway. The most promising results for large systems have been obtained, in our opinion, by perturbative vibrational approaches based on potential energy surfaces computed by hybrid (especially B3LYP) density functionals and medium size (e.g. SNSD) basis sets. In this framework, we are actively developing a comprehensive and robust computational protocol aimed at quantitative reproduction of the spectra of nucleic acid base complexes and their adsorption on solid supports (organic/inorganic). In this contribution we report the essential results of the first step devoted to isolated monomers and dimers. It is well known that in order to model the vibrational spectra of weakly bound molecular complexes dispersion interactions should be taken into proper account. In this work we have chosen two popular and inexpensive approaches to model dispersion interactions, namely the semi-empirical dispersion correction (D3) and pseudopotential based (DCP) methodologies both in conjunction with the B3LYP functional. These have been used for simulating fully anharmonic IR spectra of nucleobases and their dimers through generalized second order vibrational perturbation theory (GVPT2). We have studied, in particular, isolated adenine, hypoxanthine, uracil, thymine and cytosine, the hydrogen-bonded and stacked adenine and uracil dimers, and the stacked adenine-naphthalene heterodimer. Anharmonic frequencies are compared with standard B3LYP results and experimental findings, while the computed interaction energies and structures of complexes are compared to the best available theoretical estimates.

Onset of Chiral Adenine Surface Growth

The structure and stability of adenine crystals and thin layers has been studied by using scanning tunneling microscopy, X-ray diffraction, and density functional theory calculations. We have found that adenine crystals can be grown in two phases that are energetically quasi-degenerate, the structure of which can be described as a pile-up of 2D adenine planes. In each plane, the structure can be described as an aggregation of adenine dimers. Under certain conditions, kinetic effects can favor the growth of the less stable phase. These results have been used to understand the growth of adenine thin films on gold under ultra-high vacuum conditions. We have found that the grown phase corresponds to the α-phase, which is composed of stacked prochiral planes. In this way, the adenine nanocrystals exhibit a surface that is enantiopure. These results could open new insight into the applications of adenine in biological, medical, and enantioselective or pharmaceutical fields.

Evaluating How Discrete Water Molecules Affect Protein–DNA π–π and π+–π Stacking and T-Shaped Interactions: The Case of Histidine-Adenine Dimers

Changes in the magnitude of (M06-2X/6-31+G(d,p)) π-π stacking and T-shaped (nucleobase-edge and amino acid-edge) interactions between (neutral or protonated) histidine (His) and adenine (A) dimers upon microsolvation with up to four discrete water molecules were determined. A variety of histidine-water interactions were considered including conventional (N-H···O, N···H-O, C-H···O) hydrogen bonding and nonconventional (X-H···π (neutral His) or lone-pair···π (protonated His)) contacts. Overall, the effects of discrete His-H(2)O interactions on the neutral histidine-adenine π-π contacts are negligible (<3 kJ mol(-1) or 15%) regardless of the type of water binding, the number of water molecules bound, or the His-A dimer (stacked or (amino acid- or nucleobase-edge) T-shaped) configuration. This suggests that previously reported gas-phase binding strengths for a variety of neutral amino acid-nucleobase dimers are likely relevant for a wide variety of (microsolvated) environments. In contrast, the presence of water decreases the histidine-adenine π(+)-π interaction by up to 15 kJ mol(-1) (or 30%) for all water binding modes and orientations, as well as different stacked and T-shaped His(+)-A dimers. Regardless of the larger effect of discrete histidine-water interactions on the magnitude of the π(+)-π compared with π-π interactions, the π(+)-π binding strengths remain substantially larger than the corresponding π-π contacts. These findings emphasize the distinct nature of π(+)-π and π-π interactions and suggest that π(+)-π contacts can provide significant stabilization in biological systems relative to π-π contacts under many different environmental conditions.

Photofragmentation in selected tautomers of protonated adenine

The photofragmentation by UV excitation of selectively prepared 1(+) and 3(+) tautomers of protonated adenine is studied after excitation at a 266 and 263 nm wavelengths with two different experimental set-ups located in Seoul and Orsay. While the production of 1(+) tautomers with an electrospray ion source is now well accepted, calculations were used to ascribe the preparation of 3(+) tautomers from cold adenine dimers. The fragmentation patterns are rather similar for both tautomers, suggesting similar mechanisms as a statistical fragmentation in the ground electronic state after internal conversion.

Reference MP2/CBS and CCSD(T) quantum-chemical calculations on stacked adenine dimers. Comparison with DFT-D, MP2.5, SCS(MI)-MP2, M06-2X, CBS(SCS-D) and force field descriptions

We have performed reference quantum-chemical calculations for about 130 structures of adenine dimers in stacked conformations, with special attention given to dimers that are either vertically compressed (parallel structures) or contain close interatomic contacts (non-parallel structures). Such geometries are sampled during thermal fluctuations of nucleic acids and contribute to the local conformational variability of these systems. Their theoretical characterization requires a good description of interaction energies in the short-range repulsion region. The reference calculations have been performed with the CBS(T) method, i.e., MP2/CBS computations corrected for higher-order electron-correlation effects using the CCSD(T) method. These benchmark data have been used to examine the performance of the DFT-D, SCS(MI)-MP2, MP2.5, M06-2X and CBS(SCS-D) quantum-mechanical methods, and of the AMBER Cornell et al. force field. The present results, as well as those of our previous study on stacked uracil dimers, confirm that the force field severely exaggerates the repulsion at short intermolecular distances. This behavior complicates the use of the force field in scans of the stacking-energy dependence on local conformational parameters in nucleic acids. Compared against the previous results obtained in the uracil dimer study, the performance of DFT-D to describe stacking at short intermolecular distances has worsened, showing for the adenine dimers a larger exaggeration of the repulsion, especially for structures where the monomers are parallel to each other. Despite these deviations, the performance of DFT-D is still reasonably good and this method provides, for example, a relatively inexpensive way to monitor stacking energies along molecular dynamics trajectories. The best performers are the MP2.5, SCS(MI)-MP2, and CBS(SCS-D) methods. In addition, the energy profiles given by the SCS(MI)-MP2 and CBS(SCS-D) methods are the ones that most closely resemble the CBS(T) data. Interestingly, the performance of the SCS(MI)-MP2 method for stacked adenine dimers is better than for stacked uracil dimers, indicating that the quality of the description may vary with the nucleobase composition. Even though the SCS(MI)-MP2 method cannot match the speed of DFT-D, the results so far render it a promising tool to study intrinsic interactions in systems of moderate size. In general, for most applications all the QM methods tested here are of sufficient accuracy.

The Role of Adenine Excimers in the Photophysics of Oligonucleotides

Energies and structures of different arrangements of the stacked adenine homodimer have been computed at the ab initio CASPT2 level of theory in isolation and in an aqueous environment. Adenine dimers are shown to form excimer singlet states with different degrees of stacking and interaction. A model for a 2-fold decay dynamics of adenine oligomers can be supported in which, after initial excitation in the middle UV range, unstacked or slightly stacked pairs of nucleobases will relax by an ultrafast internal conversion to the ground state, localizing the excitation in the monomer and through the corresponding conical intersection with the ground state. On the other hand, long-lifetime intrastrand stacked excimer singlet states will be formed in different conformations, including neutral and charge transfer dimers, which originate the red-shifted emission observed in the oligonucleotide chains and that will evolve toward the same monomer decay channel after surmounting an energy barrier. By computing the transient absorption spectra for the different structures considered and their relative stability in vacuo and in water, it is concluded that in the adenine homodimers the maximum-overlap face-to-face orientations are the most stable excimer conformations observed in experiment.

Remarkably Strong T-Shaped Interactions between Aromatic Amino Acids and Adenine: Their Increase upon Nucleobase Methylation and a Comparison to Stacking.

T-shaped geometries and interaction energies between select DNA nucleobases (adenine or 3-methyladenine) and all aromatic amino acids (histidine, phenylalanine, tyrosine, or tryptophan) were examined using BSSE-corrected MP2/6-31G*(0.25) potential energy surface scans, which determined the preferred nucleobase (face)-amino acid (edge) and nucleobase (edge)-amino acid (face) interactions. The energies of dimers with the strongest interactions were further studied at the CCSD(T)/CBS level of theory, which suggests that the T-shaped interactions in adenine dimers are very strong (up to -35 kJ mol(-1)). Nucleobase methylation to form a cationic damaged base (3-methyladenine) plays a large role in the relative monomer orientations and magnitude of the interactions, which increase by 17-125%. Most importantly, this study is the first to compare the stacking and T-shaped interactions between all aromatic amino acids and select (natural and damaged) DNA nucleobases where the differences between stacking and T-shaped interactions at the CCSD(T)/CBS level are small. Therefore, our results indicate that T-shaped interactions cannot be ignored when studying biological processes, and this manuscript discusses the importance of these interactions in the context of DNA repair.

Density Functional Study of Intradimer Proton Transfers in Hydrated Adenine Dimer Ions, A2+(H2O)n (n = 0−2)

Structures of mono- and dihydrated adenine dimers and their cations were calculated using B3LYP density functional theory with the 6-31+G(d,p) basis set, in order to help understand photofragmentation experiments of hydrated adenine dimers from the energetics point of view. Several important pathways leading to the major fragmentation product, protonated adenine ion (AH(+)), thermodynamically at minimum costs were investigated at the ground-state electronic potential surface of hydrated adenine dimer cations. Our calculations suggest that the proton transfer from one adenine moiety to the other in hydrated dimer ions readily occurs with negligible barriers in normal hydration conditions. In asymmetrically hydrated ions, however, the proton transfer to more hydrated adenine moieties is kinetically hindered due to heightened transition-state barriers, while the other way is still barrierless. Such directional preference in proton transfer may be characterized as a unique dimer ion property, stemming from the difference in basicity of the two nitrogen atoms involved in the double hydrogen bond that would be equivalent without hydration. We also found that dimer cleavage requires about 4 times larger energy than evaporation of individual water molecules, so it is likely that most solvent molecules evaporate before the eventual dimer cleavage when available internal energy is limited.

Electronic Structure of Adenine and Thymine Base Pairs Studied by Femtosecond Electron−Ion Coincidence Spectroscopy

Femtosecond pump-probe spectroscopy was combined with photoelectron-photoion coincidence detection to investigate the electronic structure and dynamics of isolated adenine (A) and thymine (T) dimers and the adenine-thymine (AT) base pair. The photoelectron spectra show that pipi* and npi* states are only weakly perturbed in the hydrogen-bound dimers as compared to the monomers. For cationic base pairs with internal energies greater than 1 eV, we observed considerable cluster fragmentation into protonated monomers. This process selectively removed signals from the npi* --> n-1 ionization channel in all dimers. The photoelectron spectra are compared to time-resolved mass spectra and confirm the assignment of short-lived pipi* and npi* populations in the adenine, thymine, and mixed AT dimers.

First-principle applications on protonated adenine dimers with N–H⋯N intermolecular coupling

The geometry structures and the electronic properties of protonated adenine dimer were described employing the first-principle hybrid density functional theory with a special emphasis on N–H⋯N hydrogen bond coupling mode complexes. Nine stable isomers with N–H⋯N type H-bonds were characterized on its potential energy surface and classified as planar and nonplanar two types. The optimized geometry and the Mulliken charge populations indicate that the spontaneous proton transfer occurs in some isomers owning to strong N–H⋯N coupling. The stabilization energies corrected with the basis set superposition error ( E BSSE) and zero-point vibrational energy (ZPVE) demonstrate the following significant points. First, the impact factors of E BSSE, such as geometry structure, the number of H-bonds, and proton situation, are proposed. Second, the dimers with the N 1 ′ – H 15 ′ ⋯ N type H-bonds are more stable than others. Third, the whole system should be included to analyze the stability of the dimer.

Experimental and theoretical identification of adenine monolayers on Ag-terminated Si(111)

We discuss the formation of ordered adenine monolayers on a Ag-Si (111)root 3 x root 3R30 degrees surface held at room temperature. Adenine forms a hexagonal network, which is stabilized by hydrogen bonding. By comparing the ordering observed using scanning tunneling microscopy with molecular superstructures calculated using ab initio density functional theory, we suggest that adenine forms an array of irregular hexagons through the combination of three distinct hydrogen-bonded adenine dimers. The dimensions of the ordered monolayer are commensurate with the surface with root 3 x root 7 periodicity relative to the Ag-Si (111) root 3 x root 3R30 degrees reconstruction, in agreement with calculated molecular dimensions. An analysis of double-domain regions leads us to conclude that there are two domains with different chiralities.

Dynamics of Photoinduced Processes in Adenine and Thymine Base Pairs

The excited-state dynamics of adenine and thymine dimers and the adenine-thymine base pair were investigated by femtosecond pump-probe ionization spectroscopy with excitation wavelengths of 250-272 nm. The base pairs showed a characteristic ultrafast decay of the initially excited pi pi* state to an n pi* state (lifetime tau(pi pi*) approximately 100 fs) followed by a slower decay of the latter with tau(n pi*) approximately 0.9 ps for (adenine)2, tau(n pi*) = 6-9 ps for (thymine)2, and tau(n pi*) approximately 2.4 ps for the adenine-thymine base pair. In the adenine dimer, a competing decay of the pi pi* state via the pi sigma* state greatly suppressed the n pi* state signals. Similarities of the excited-state decay parameters in the isolated bases and the base pairs suggest an intramonomer relaxation mechanism in the base pairs.