The post-SCF quantum chemistry characteristics of the energetic heterogeneity of stacked guanine–guanine pairs found in B-DNA and A-DNA crystals

Theoretical investigation of the coupling between hydrogen-atom transfer and stacking interaction in adenine-thymine dimers

Three different dimers of the adenine-thymine (A-T) base pair are studied to point out the changes of important properties (structure, atomic charge, energy and so on) induced by coupling between the movement of the atoms in the hydrogen bonds and the stacking interaction. The comparison of these results with those for the A-T monomer system explains the role of the stacking interaction in the hydrogen-atom transfer in this biologically important base pair. The results support the idea that this coupling depends on the exact dimer considered and is different for the N-N and N-O hydrogen bonds. In particular, the correlation between the hydrogen transfer and the stacking interaction is more relevant for the N-N bridge than for the N-O one. Also, the two different mechanisms of two-hydrogen transfer (step by step and concerted) can be modified by the stacking interaction between the base pairs.

Combined effects of π-π stacking and hydrogen bonding on the (N1) acidity of uracil and hydrolysis of 2'-deoxyuridine

M06-2X/6-31+G(d,p) is used to study the simultaneous effects of π-π stacking interactions with phenylalanine (modeled as benzene) and hydrogen bonding with small molecules (HF, H(2)O, and NH(3)) on the N1 acidity of uracil and the hydrolytic deglycosylation of 2'-deoxyuridine (dU) (facilitated by fully (OH(-)) or partially (HCOO(-)···H(2)O) activated water). When phenylalanine is complexed with isolated uracil, the proton affinity of all acceptor sites significantly increases (by up to 28 kJ mol(-1)), while the N1 acidity slightly decreases (by ~6 kJ mol(-1)). When small molecules are hydrogen bound to uracil, addition of the phenylalanine ring can increase or decrease the acidity of uracil depending on the number and nature (acidity) of the molecules bound. Furthermore, a strong correlation between the effects of π-π stacking on the acidity of U and the dU deglycosylation reaction energetics is found, where the hydrolysis barrier can increase or decrease depending on the nature and number of small molecules bound, the nucleophile considered (which dictates the negative charge on U in the transition state), and the polarity of the (bulk) environment. These findings emphasize that the catalytic (or anticatalytic) role of the active-site aromatic amino acid residues is highly dependent on the situation under consideration. In the case of uracil-DNA glycosylase (UNG), which catalyzes the hydrolytic excision of uracil from DNA, the type of discrete hydrogen-bonding interactions with U, the nature of the nucleophile, and the anticipated weak, nonpolar environment in the active site suggest that phenylalanine will be slightly anticatalytic in the chemical step, and therefore experimentally observed contributions to catalysis may entirely result from associated structural changes that occur prior to deglycosylation.

Structural and energetic consequences of oxidation of d(ApGpGpGpTpT) telomere repeat unit in complex with TRF1 protein

The configuration hyperspace of canonical and oxidized 14-mers of B-DNA comprising telomere repeat units d(ApGpGpGpTpT) was sampled over 40 ns via molecular dynamic (MD) simulations. The energetic and structural consequences of TRF1 binding to telomere B-DNA were compared with non-complexed systems. Energetic properties of analyzed pairs, di- and tri-nucleotide steps occurring in central telomere repeat unit were estimated by means of advanced quantum chemistry computations including not only BSSE corrections, electron correlation contributions but also non-negligible many-body terms. These data along with bases pair and base step parameters distributions allow for quantization of consequences of oxidation and/or TRF1 binding to telomere repeat units. Occurrence of 8-oxoguanine in central telomeric triad (CTT) is the source of high stiffness if compared to non-modified oligomer. The origin of this property comes from significantly alteration of intermolecular interactions introduced by 8-oxoguanine. The increased stability observed for base-base interactions are accumulated and characterizes also di- and tri-nucleotides. The observed changes in the intermolecular interactions originate from structural alterations imposed by TRF1 binding to canonical and oxidized telomere B-DNA. First and most direct consequence of TRF1 binding to oxidized telomere repeat unit is alteration of shift-slide correlations if compared to canonical system. This in turn leads to large differences in purine-purine overlapping in oxidized structures. Thus, oxidized telomere B-DNA double strands are sensitive to interactions with protein ligands and numerous structural and energetic changes are imposed on base pairs forming CTT.

Non-additive interactions of nucleobases in model dinucleotide steps occurring in B-DNA crystals

Non-additivity of base-base interactions in all ten possible model dinucleotide steps were analyzed on MP2/aug-cc-pvDZ quantum chemistry level. Conformations of four nucleobases exactly matched to ones occurring in B-DNA crystals. In most of thw 162 analyzed tetramers both three- and four-body contributions are negligible except for d(GpG) steps. However, in these dinucleotides both contributions are always of opposite signs and in all cases the sum of all non-additive part of intermolecular interactions do not exceed 2.6 kcal mol(-1). This stands for less than 5% of the overall binding energy of dinucleotide steps. Also replacements of guanine with 8-oxoguanine in d(GpG) systems introduces non-additivity of the same magnitude as for canonical dinucleotides. It is observed linear relationships between values of total binding energy obtained in the tetramer basis set and estimated energy exclusively in dimers basis sets with assumption of pairwise additivities. For all analyzed dinucleotides steps there are also linear correlations between amount of non-additive contributions and magnitude of pairs interactions. Based on differences in electrostatic contribution to the total binding energy of four nucleobases and polarity of dinucleotide steps three distinct classes of dinucleotide steps were identified.

Comparison of intrinsic stacking energies of ten unique dinucleotide steps in A-RNA and B-DNA duplexes. Can we determine correct order of stability by quantum-chemical calculations?

High level ab initio methods have been used to study stacking interactions in ten unique base pair steps both in A-RNA and in B-DNA duplexes. The protocol for selection of geometries based on molecular dynamics (MD) simulations is proposed, and its suitability is demonstrated by comparison with stacking in steps at fiber diffraction geometries. It is shown that fiber diffraction geometries are not sufficiently accurate for interaction energy calculations. In addition, the protocol for selection of geometries based on MD simulations allows for the evaluation of the variability of the intrinsic stacking energies along the MD trajectories. The uncertainty in stacking energies (difference between the most and least stable geometry) due to the dynamical nature of systems can be, in some cases, as large as 3.0 kcal x mol(-1), which is almost 50% of the actual sequence dependence of base stacking energies (the energy difference between the most and least stable sequences). Thus, assessing the relative magnitude of the gas phase stacking energy using a single geometry for each sequence is insufficient to obtain an unambiguous order of gas phase stacking energies in canonical double helices. Though the ordering of ten unique dinucleotide steps cannot be definitive, some general conclusions were drawn. The stacking energies of base pair steps in A-RNA are more evenly separated compared to B-DNA, and their ordering is less sensitive to the dynamics of the system compared to be B-DNA. The most stable step both in B-DNA and A-RNA is the GC/GC [corrected] step that is well separated from the second most stable step CG/CG. [corrected] Also the least stable step (the CC/GG step) is well separated from the rest of the structures. The calculations further show that B-DNA stacking is favorable only marginally (on average by 1.14 kcal x mol(-1) per base pair step) over A-RNA stacking, and this difference vanishes after subtracting the stabilizing van der Waals effect of the thymine 5-methyl group that is absent in RNA. Basically, no correlation between the sequence dependence of gas phase stacking energies and the sequence dependence of DeltaG degrees(37) free energies used in nearest-neighbor models was found either for B-DNA or for A-RNA. This reflects the complexity of the balance of forces that are responsible for the sequence dependence of thermodynamics stability of nucleic acids, which masks the effect of the intrinsic interactions between the stacked base pairs.

Effects of the biological backbone on DNA-protein stacking interactions

The pi-pi stacking (face-to-face) interactions between the five natural DNA or RNA nucleobases and the four aromatic amino acids were compared using three different types of dimers: (1) a truncated nucleoside (nucleobase) stacked with a truncated amino acid; (2) a truncated nucleoside (nucleobase) stacked with an extended amino acid; and (3) a nucleoside (extended nucleobase) stacked with a truncated amino acid. Systematic (MP2/6-31G*(0.25)) potential energy surface scans reveal important information about the effects of the deoxyribose sugar and protein backbone on the structure and binding energy between truncated nucleobase and amino acid models that are typically implemented in the literature. Most notably, electrostatic and steric interactions arising from the bulkiness of the biological backbones can change the preferred relative orientations of DNA and protein pi-systems. More importantly, the protein backbone can strengthen the stacking energy (by up to 10 kJ mol(-1)), while the deoxyribose moiety can strengthen or weaken the stacking interaction depending on the positioning of the amino acid relative to the sugar residue. These effects are likely due to additional interactions between the amino acid or nucleobase ring and the backbone in the extended monomer rather than significant changes in the properties of the biological pi-systems upon model extension. Since the present work reveals that all calculated DNA-protein stacking interactions are significant and approach the strength of other noncovalent interactions between biomolecules, both pi-pi and backbone-pi interactions must be considered when attempting to gain a complete picture of DNA-protein binding.

The post-SCF quantum chemistry characteristics of inter- and intra-strand stacking interactions in d(CpG) and d(GpC) steps found in B-DNA, A-DNA and Z-DNA crystals

The energies of intra- and inter-strand stacking interactions in model d(GpC) and d(CpG) two-base-pair steps were estimated by MP2/aug-cc-pVDZ single point calculations corrected for basis superposition errors. The stacked two-nucleobase pairs were constructed using experimental values of base pair and base step parameters taken from Nucleic Acid Database (http://ndbserver.rutgers.edu/). Three distinct polymorphic forms were analysed, namely A-, B- and Z-DNA. The applied methodology enables statistical analysis of structural and energetic diversities. The structural relationships between polymorphic forms are quite complex and depend on the sequence of pairs. The variability of parameters such as shift and tilt is almost the same irrespective of the polymorphic form and sequence of steps analysed. In contrast, shift and twist distributions easily discriminate all three polymorphic forms of DNA. Interestingly, despite significant structural diversities, the energies of the most frequent energy ranges are comparable irrespective of the polymorphic form and base sequence. There was observed compensation of inter- and intra-strand interactions, especially for d(GpC) and d(CpG) steps found in A- and B-DNA. Thus, among many other roles, these pairs act as a kind of energetic buffer, balancing the double helix.