Nature of Nucleic Acid−Base Stacking: Nonempirical ab Initio and Empirical Potential Characterization of 10 Stacked Base Dimers. Comparison of Stacked and H-Bonded Base Pairs

Effects of 5-[S-(2,4-dinitrophenyl)-thio]-2'-deoxyuridine analog incorporation on the structure and stability of DNA hybrids: implications for the design of nucleic acid probes

Labeled nucleic acid probes are used as diagnostic tools by detecting changes in gene expression upon hybridization to target RNAs or DNAs that are related to specific disease genes. 5-[S-(2, 4-Dinitrophenyl)-thio]-2'-deoxyuridine analog represents an excellent nucleic acid label, containing the DNP group which functions both as a probe and as a precursor for the introduction of a variety of fluorescent groups. This study describes thermal denaturation hybridization experiments with oligonucleotides containing the 5-[S-(2,4-dinitrophenyl)-thio]-2'-deoxyuridine analog. Using molecular modeling techniques, the effects of this analog on the hybrid structure and stability were examined, including (i) analog conformation, (ii) hydrogen bonding, (iii) stacking interactions and (iv) hybrid helical geometry. This analog does not prohibitively affect the hybrid thermal stability and incorporation of the analog does not compromise the structural integrity of the double helix. In particular, the sequence-dependence of the analog effects and the dependence on the modification site relative to the end(s) of the helix were investigated. Findings described here should provide guidelines in the rational design of nucleic acid probes.

Interactions of hydrated IIa and IIb group metal cations with thioguanine-cytosine DNA base pair: Ab initio and density functional theory investigation of polarization effects, differences among cations, and flexibility of the cation hydration shell

The structures and energies of the thioguanine-cytosine Watson-Crick (thioGC WC) base pair interacting with hydrated IIa (Mg2+, Ca2+, Ba2+) and IIb group (Zn2+, Cd2+, Hg2+) cations have been studied using ab initio techniques. Furthermore, complexes between guanine and thioguanine with hydrated cations have been characterized assuming various structures of the hydration shells. The complexes of the thioGC WC base pair with hydrated cations have similar properties as the previously studied GC WC base pair. There is substantial polarization stabilization of the base pairing due to cation binding which amounts to 7 - 11 kcal/mol. Soft Cd2+ and Hg2+ cations have a uniquely strong interaction with the thiogroup and induce substantial nonplanarity of the pairing. The thiogroup tends to reduce the number of water molecules in the first hydration shell of the cation. All complexes were optimized within the Hartree-Fock (HF) approximation while their energetics has been evaluated using the second-order Moller-Plesset perturbational method (MP2). All interaction energy evaluations and a substantial portion of the optimizations of the hydrated cation-(thio)guanine complexes have been repeated using Becke-3LYP Density Functional Theory method. All three approximations used (HF, Becke-3LYP, and MP2) give qualitatively the same results for the present cationic complexes. The results demonstrate specific differences among the cations and provide a set of reference structures and energies for verification and/or parametrization of empirical potentials and other theoretical methods.

Interaction and solvation energies of nonpolar DNA base analogues and their role in polymerase insertion fidelity

Although DNA polymerase fidelity has been mainly ascribed to Watson-Crick hydrogen bonds, two nonpolar isosteres for thymine (T) and adenine (A)--difluorotoluene (F) and benzimidazole (Z) --effectively mimic their natural counterparts in polymerization experiments with pol I (KF exo-) [JC Morales and ET Kool. Nature Struct Biol, 5, 950-954, 1998]. By ab initio quantum chemical gas phase methods (HF/6-31G* and MP2/6-31G**) and a solvent phase method (CPCM-HF/6-31G**), we find that the A-F interaction energy is 1/3 the A-T interaction energy in the gas phase and unstable in the solvent phase. The F-Z and T-Z interactions are very weak and T-Z is quite unstable in the solvent. Electrostatic solvation energy calculations on F, Z and toluene yield that Z is two times, and F and toluene are five times, less hydrophilic than the natural bases. Of the new "base-pairs" (F-Z, T-Z, and F-A), only F-A formed an A-T-like arrangement in unconstrained optimizations. F-Z and T-Z do not freely form planar arrangements, and constrained optimizations show that large amounts of energy are required to make these pairs fit the exact A-T geometry, suggesting that the polymerase does not require all bases to conform to the exact A-T geometry. We discuss a model for polymerase/nucleotide binding energies and investigate the forces and conformational range involved in the polymerase geometrical selection.

A study of the hydration of deoxydinucleoside monophosphates containing thymine, uracil and its 5-halogen derivatives: Monte Carlo simulation

An extensive Monte Carlo simulation of hydration of various conformations of the dinucleoside monophosphates (DNP), containing thymine, uracil and its 5-halogen derivatives has been performed. An anti-anti conformation is the most energetically stable one for each of the DNPs. In the majority of cases the energy preference is determined by water-water interaction. For other dimers conformational energy is the most important factor, or both the factors are of nearly equal importance. The introduction of the methyl group into the 5-position of uracil ring most noticeably influences the conformational energy and leads to the decrease of its stabilizing contribution to the total interaction energy. The introduction of halogen atoms increases the relative content of anti-syn and syn-anti conformations of DNPs as compared to the parent ones due to the formation of an energetically more favorable water structure around these conformations. A correlation is observed between the Monte Carlo results for the halogenated DNPs and their experimental photoproduct distribution. The data obtained demonstrates a sequence dependence in the photochemistry of the halogenated dinucleoside monophosphates.

Use of a 3D structure data base for understanding sequence-dependent conformational aspects of DNA

The roll-twist-slide correlation in the DNA crystal structures that are collected in the Nucleic Acid Data Base is analyzed in order to obtain a general understanding of the effects of the nucleotide sequence on the 3D structure of a dinucleotide step. It is concluded that the differences between the pyrimidine bases and the purine bases in terms of their physical shapes are the major factors that determine the stereochemical characteristics of the steps through base to backbone and base to base interactions. The characteristics are further modulated by the differences between the A:T and G:C base-pairs, which can be explained by enhancement of the purine-pyrimidine asymmetry in the A:T base-pair.

Protein-nucleic acid recognition: simulation of base and "model" amino acids complexes in DMSO by the Monte Carlo method

A computer simulation of guanine (G), cytosine (C), the G-C base pair, protonated C (CH+), acetic acid in neutral (AcOH) and deprotonated (AcO-) forms, G-AcO-, C-AcOH, and CH(+)-AcO- complexes, solvated in DMSO was carried out by the Monte Carlo method. It is shown that the G-C base pair formation in DMSO is energetically favorable. The G-AcO- complex formation is comparable with the formation of G-C base pair in energetically favorability. In this case the acetate anion can replace C in the G-C base pair. The formation of the C-AcOH complex is much less favorable than the formation of the G-C pair. However proton transfer from AcOH to C leads to the formation of the CH(+)-AcO- complex, which is the most favorable of all complexes studied. Here the acetic acid can replace G in a G-C base pair. The formation of G-AcO- and CH(+)-AcO- specific complexes detected in DMSO with the help of experiment and theory is a competitive process with respect to the formation of G-C base pairs, and can be considered the primary step in the real mechanism of protein-nucleic acid recognition.

The study of the stability of Watson-Crick nucleic acid base pairs in water and dimethyl sulfoxide: computer simulation by the Monte Carlo method

An extensive computer simulation of nucleic acid bases and Watson-Crick base pairs in a water cluster and DMSO cluster is performed by the Monte Carlo method. It is demonstrated that the unfavorable energetics of pair formation in a water cluster is determined by the significant destabilizing contribution of solvent to the energy of complex formation. It is shown that the formation of coplanar base pairs in a DMSO cluster is favorable. The DMSO cluster stabilizes A-U and A-T base pairs and the insignificant destabilization of the G-C base pair by a DMSO cluster is much less than the stabilization which occurs due to the attraction between bases.

Hydrogen-bonded trimers of DNA bases and their interaction with metal cations: ab initio quantum-chemical and empirical potential study

Neutral (G.GC, A.AT, G.AT, T.AT, and C(imino).GC) and protonated (CH+.GC and AH+.GC) hydrogen-bonded trimers of nucleic acid bases were characterized by ab initio methods with the inclusion of electron correlation. In addition, the influence of metal cations on the third-strand binding in Purine-Purine-Pyrimidine (Pu.PuPy) reverse-Hoogsteen triplets has been studied. The ab initio calculations were compared with those from recently introduced force fields (AMBER4.1, CHARMM23, and CFF95). The three-body term in neutral trimers is mostly negligible, and the use of empirical potentials is justified. The only exception is the neutral G.GC Hoogsteen trimer with a three-body term of -4 kcal/mol. Protonated trimers are stabilized by molecular ion-molecular dipole attraction and the interaction within the complex is nonadditive, with the three-body term on the order of -3 kcal/mol. There is a significant induction interaction between the third-strand protonated base and guanine. The calculations indicate an enhancement of the third-strand binding in the G.GC reverse-Hoogsteen trimer due to-metal cation coordination to the N7/O6 position of the third-strand guanine. Interactions between metal cations and complexes of DNA bases are in general highly non-additive; the three-body term is above-10 kcal/mol in a complex of a divalent cation (Ca2+) with the GG reverse-Hoogsteen pair. The pairwise additive empirical potentials qualitatively underestimate the binding energy between cation and base.

Hydrogen bonding and stacking of DNA bases: a review of quantum-chemical ab initio studies

Ab initio quantum-chemical calulations with inclusion of electron correlation made since 1994 (such reliable calculations were not feasible before) significantly modified our view on interactions of nucleic acid bases. These calculations allowed to perform the first reliable comparison of the strength of stacked and hydrogen bonded pairs of nucleic acid bases, and to characterize the nature of the base-base interactions. Although hydrogen-bonded complexes of nucleobases are primarily stabilized by the electrostatic interaction, the dispersion attraction is also important. The stacked pairs are stabilized by dispersion attraction, however, the mutual orientation of stacked bases is determined rather by the electrostatic energy. Some popular theories of stacking were ruled out: The theory based on attractive interactions of polar exocyclic groups of bases with delocalized electrons of the aromatic rings (Bugg et al., Biopolymers 10, 175 (1971), and the pi-pi interactions model (C.A. Hunter, J. Mol. Biol. 230, 1025 (1993)). The calculations demonstrated that amino groups of nucleobases are very flexible and intrinsically nonplanar, allowing hydrogen-bond-like interactions which are oriented out of the plane of the nucleobase. Many H-bonded DNA base pairs are intrinsically nonplanar. Higher-level ab initio calculations provide a unique set of reliable and consistent data for parametrization and verification of empirical potentials. In this article, we present a short survey of the recent calculations, and discuss their significance and limitations. This summary is written for readers which are not experts in computational quantum chemistry.