Stability of nucleic acid base pairs in organic solvents: molecular dynamics, molecular dynamics/quenching, and correlated ab initio study

The dynamic structure and potential energy surface of adenine...thymine and guanine...cytosine base pairs and their methylated analogues interacting with a small number (from 1 to 16 molecules) of organic solvents (methanol, dimethylsulfoxide, and chloroform) were investigated by various theoretical approaches starting from simple empirical methods employing the Cornell et al. force field to highly accurate ab initio quantum chemical calculations (MP2 and particularly CCSD(T) methods). After the simple molecular dynamics simulation, the molecular dynamics in combination with quenching technique was also used. The molecular dynamics simulations presented here have confirmed previous experimental and theoretical results from the bulk solvents showing that, whereas in chloroform the base pairs create hydrogen-bonded structures, in methanol, stacked structures are preferred. While methanol (like water) can stabilize the stacked structures of the base pairs by a higher number of hydrogen bonds than is possible in hydrogen-bonded pairs, the chloroform molecule lacks such a property, and the hydrogen-bonded structures are preferred in this solvent. The large volume of the dimethylsulfoxide molecule is an obstacle for the creation of very stable hydrogen-bonded and stacked systems, and a preference for T-shaped structures, especially for complexes of methylated adenine...thymine base pairs, was observed. These results provide clear evidence that the preference of either the stacked or the hydrogen-bonded structures of the base pairs in the solvent is not determined only by bulk properties or the solvent polarity but rather by specific interactions of the base pair with a small number of the solvent molecules. These conclusions obtained at the empirical level were verified also by high-level ab initio correlated calculations.     

Hydration and stability of nucleic acid bases and base pairs

Empirical, quantum chemical calculations and molecular dynamics simulations of the role of a solvent on tautomerism of nucleic acid bases and structure and properties of nucleic acid base pairs are summarized. Attention was paid to microhydrated (by one and two water molecules) complexes, for which structures found by scanning of empirical potential surfaces were recalculated at a correlated ab initio level. Additionally, isolated as well as mono- and dihydrated H-bonded, T-shaped and stacked structures of all possible nucleic acid base pairs were studied at the same theoretical levels. We demonstrate the strong influence of a solvent on the tautomeric equilibrium between the tautomers of bases and on the spatial arrangement of the bases in a base pair. The results provide clear evidence that the prevalence of either the stacked or hydrogen-bonded structures of the base pairs in the solvent is not determined only by its bulk properties, but rather by specific hydrophilic interactions of the base pair with a small number of solvent molecules.     

Potential energy surfaces of the microhydrated guanine...cytosine base pair and its methylated analogue.

A complete scan of the potential and free-energy surfaces of monohydrated and dihydrated guanine...cytosine and 9-methylguanine...1-methylcytosine base pairs was realized by the molecular dynamics/quenching technique using the force field of Cornell et al. implemented in the AMBER7 program. The most stable and populated structures localized were further fully reoptimized at the correlated ab initio level employing the resolution of identity M?ller-Plesset method with a large basis set. A systematic study of microhydration of these systems using a high-level correlated ab initio approach is presented for the first time. The different behavior of guanine...cytosine and adenine...thymine complexes is also discussed. These studies of nucleic acid base pairs are important for finding binding sites of water molecules around bases and for better understanding of the influence of the solvent on the stability of the structure of DNA.     

Properties of the optimal environment inducing tautomerization of complementary base pairs

The optimal environment charge configurations are predicted for the tautomerization of complementary base pairs into their corresponding rare forms, and vice versa. Results indicate that cations approaching the N₃ guanine site may induce tautomerization of the normal guanine—cytosine (G---C) base pair into its rare form. The reverse process requires that the cation approach the O₂ thymine site of the rare adenine*—thymine* pair (A*---31T*) or the O₆ guanine site of rare guanine*—cytosine* base pair (G*---C*). Possible mutagenic and antimutagenic roles of metal cations approaching base pairs are also discussed.     

Structural and electronic property changes of the nucleic acid bases upon base pair formation

A systematic study of structures and electronic properties has been carried out for the nucleic acid bases adenine, guanine, thymine, and cytosine and for the base pairs adenine–thymine and guanine–cytosine. We focus our attention on these properties, which experience significant changes when single nucleic bases join to form base pairs. Such properties are expected to play an important role during the formation of the DNA molecule in its B conformation. All-electron calculations with inclusion of correlation effects were performed according to the local and nonlocal density functional approaches. We compare our results with previous ab initio and semiempirical values and with available experimental data. Advantages and disadvantages for these density functional-based methods are discussed. We conclude that applications of such models to investigate larger compounds of a similar nature are promising. © 1994 by John Wiley & Sons, Inc.     

Factors determining the deriving force of DNA formation: geometrical differences of base pairs, dehydration of bases, and the arginine assisting

The mechanism of the fidelity synthesis of DNA associated with the process of dGTP combination to the DNA template was explored. The exclusion of water molecules from the hydrated DNA bases can amplify the energy difference between the correct and incorrect base pairs, but the effect of the water molecules on the Gibbs free energy of formation is dependent on the binding sites for the water molecules. The water detachment from the incoming dNTP is not the only factor but the first step for the successful replication of DNA. The second step is the selection of the DNA polymerase on the DNA base pair through the comparison between the correct DNA base and the incorrect DNA base. The bonding of the Arg668 with the incoming dNTP can enlarge the Gibbs free energies of formation of the base pairs, especially the correct base pairs, thus increasing the driving force of DNA formation. When the DNA base of the primer terminus is correct, the extension of the guanine and the adenine is quicker than that of the cytosine and the thymine because of the hydrogen bonding fork formation of Arg668 with the minor groove of the primer terminus and the ring oxygen of the deoxyribose moiety of the incoming dNTP. Because of the geometry differences of the incorrect base pairs with the correct base pairs, the effect from the DNA polymerase is smaller on the incorrect base pair than on the correct base pair, and the extension of a mispair is slower than that of a correct base pair. This decreases the extension rate of the base pair and thus allows proofreading exonuclease activity to excise the incorrect base pair. Arg668 cannot prevent the extension of the GT mispair, as well as the GC correct base pair, and GA and GG mispairs. This may be attributed to the small geometry difference between the GT base pair and the correct AT base pair.     

The nature of base stacking: a Monte Carlo study

To elucidate the physical origin of the preference of nucleic acid bases for stacking over hydrogen bonding in water, Monte Carlo simulations were performed starting from Watson?CCrick structures of the adenine?Cthymine, adenine?Curacil and guanine?Ccytosine base pairs, as well as from the Hoogsteen adenine?Cthymine base pair, in clusters comprising 400 and 800 water molecules. The simulations employed a newly implemented Metropolis Monte Carlo algorithm based on the extended cluster approach. All simulations reached stacked structures, confirming that such structures are preferred over the hydrogen-bonded Watson?CCrick and Hoogsteen base pairs. The Monte Carlo simulations show the complete transition from hydrogen-bonded base pairs to stacked structures in the Monte Carlo framework. Analysis of the average energies shows that the preference of stacked over hydrogen-bonded structures is due to the increased water?Cbase interaction in these structures. This is corroborated by the increased number of water?Cbase hydrogen bonds in the stacked structures.     

Halogen bonding in DNA base pairs

Halogen bonding (R-X···Y) is a qualitative analogue of hydrogen bonding that may prove useful in the rational design of artificial proteins and nucleotides. We explore halogen-bonded DNA base pairs containing modified guanine, cytosine, adenine and thymine nucleosides. The structures and stabilities of the halogenated systems are compared to the normal hydrogen bonded base pairs. In most cases, energetically stable, coplanar structures are identified. In the most favorable cases, halogenated base pair stabilities are within 2 kcal mol(-1) of the hydrogen bonded analogues. Among the halogens X = Cl, Br, and I, bromine is best suited for inclusion in these biological systems because it possesses the best combination of polarizability and steric suitability. We find that the most stable structures result from a single substitution of a hydrogen bond for a halogen bond in dA:dT and dG:dC base pairs, which allows 1 or 2 hydrogen bonds, respectively, to complement the halogen bond.     

Stabilization energies of the hydrogen-bonded and stacked structures of nucleic acid base pairs in the crystal geometries of CG, AT, and AC DNA steps and in the NMR geometry of the 5'-d(GCGAAGC)-3' hairpin: Complete basis set calculations at the MP2 and CCSD(T) levels

Stabilization energies of the H-bonded and stacked structures of a DNA base pair were studied in the crystal structures of adenine-thymine, cytosine-guanine, and adenine-cytosine steps as well as in the 5'-d(GCGAAGC)-3' hairpin (utilizing the NMR geometry). Stabilization energies were determined as the sum of the complete basis set (CBS) limit of MP2 stabilization energies and the Delta E(CCSD(T)) - Delta E(MP2) correction term evaluated with the 6-31G*(0.25) basis set. The CBS limit was determined by a two-point extrapolation using the aug-cc-pVXZ basis sets for X = D and T. While the H-bonding energies are comparable to those of base pairs in a crystal and a vacuum, the stacking energies are considerably smaller in a crystal. Despite this, the stacking is still important and accounts for a significant part of the overall stabilization. It contributes equally to the stability of DNA as does H-bonding for AT-rich DNAs, while in the case of GC-rich DNAs it forms about one-third of the total stabilization. Interstrand stacking reaches surprisingly large values, well comparable to the intrastrand ones, and thus contributes significantly to the overall stabilization. The hairpin structure is characterized by significant stacking, and both guanine...cytosine pairs possess stacking energies larger than 11.5 kcal/mol. A high portion of stabilization in the studied hairpin comes from stacking (similar to that found for AT-rich DNAs) despite the fact that it contains two GC Watson-Crick pairs having very large H-bonding stabilization. The DFT/B3LYP/6-31G** method yields satisfactory values of interaction energies for H-bonded structures, while it fails completely for stacking.     

Copolymers containing alternating sequences of nucleic acid base pairs. III. Spectroscopic studies

A study was conducted of the complementary base pair interactions between various pairs of electron-donor monomers, electron-acceptor monomers, homopolymers and alternating copolymers selected from the following group: ( 1 ) 9-(2-vinyloxyethyl)adenine; ( 2 ) 1-(2-vinyloxyethyl)thymine; ( 3 ) 1-(2-vinyloxyethyl)cytosine; ( 4 ) 9-(2-maleimidoethyl)adenine; ( 5 ) 6-chloro-9-(2-maleimidoethyl)purine; ( 6 ) 1-(2-maleimidoethyl)thymine; ( 7 ) 1-(2-maleimidoethyl)cytosine; ( 8 ) homopolymer of ( 4 ); ( 9 ) homopolymer of ( 6 ); ( 10 ) alternating copolymer of ( 2 ) and maleic anhydride; ( 11 ) alternating copolymer of ( 2 ) and ( 5 ); and ( 12 ) alternating copolymer of ( 2 ) and ( 4 ). By ¹H-NMR, in CDCL₃, the base pair interactions between ( 1 ) and ( 2 ) were shown to be hydrogen bonding, the extent of which was shown by a calculated binding constant, K = 61.81 L/mol. The nature of this interaction was conformed by IR. Neither monomer pairs ( 1 )/( 2 ) nor ( 4 )/( 6 ) exhibited hydrogen bonding in DMSO-d₆. However, hydrogen bonding interaction was observed for DMSO-d₆ solutions of homopolymers ( 8 ) and ( 9 ) and for alternating copolymer ( 12 ). On the basis of an upfield chemical shift of the 2- and 8-aromatic protons of ademine of ( 1 ) in D₂O, a partial overlap stacking interaction is proposed. No charge-transfer interactions could be observed by UV between donor-acceptor monomer pairs.     

Molecular orbital theory of the hydrogen bond : XXXII. The effect of H and Li association on the A---T and G---C pairs

Ab initio molecular orbital calculations have been performed to determine the structures and stabilization energies of the A---T and G---C base pairs and their complexes with H⁺ and Li⁺, H⁺ and Li⁺ association stabilizes the A---T pair except for Li⁺ association at O₄ in thymine. Protonation of thymine stabilizes the A---T pair to a greater extent than protonation of adenine. The association of H⁺ and Li⁺ with guanine stabilizes the G---C pair, but protonation of cytosine destabilizes G---C. Changes in the structures of the hydrogen bonds in the A---T and G---C pairs reflect changes in hydrogen bond strengths.     

Accurate interaction energies of hydrogen-bonded nucleic acid base pairs

Hydrogen-bonded nucleic acids base pairs substantially contribute to the structure and stability of nucleic acids. The study presents reference ab initio structures and interaction energies of selected base pairs with binding energies ranging from -5 to -47 kcal/mol. The molecular structures are obtained using the RI-MP2 (resolution of identity MP2) method with extended cc-pVTZ basis set of atomic orbitals. The RI-MP2 method provides results essentially identical with the standard MP2 method. The interaction energies are calculated using the Complete Basis Set (CBS) extrapolation at the RI-MP2 level. For some base pairs, Coupled-Cluster corrections with inclusion of noniterative triple contributions (CCSD(T)) are given. The calculations are compared with selected medium quality methods. The PW91 DFT functional with the 6-31G basis set matches well the RI-MP2/CBS absolute interaction energies and reproduces the relative values of base pairing energies with a maximum relative error of 2.6 kcal/mol when applied with Becke3LYP-optimized geometries. The Becke3LYP DFT functional underestimates the interaction energies by few kcal/mol with relative error of 2.2 kcal/mol. Very good performance of nonpolarizable Cornell et al. force field is confirmed and this indirectly supports the view that H-bonded base pairs are primarily stabilized by electrostatic interactions.     

Structures of nucleobases trapped within Au triangles and its effects on hydrogen bonding in base pairs of DNA

Nucleobases (adenine (A), thymine (T), cytosine (C), and guanine (G)) trapped within two metal clusters such as Au(3) undergo expansion. Our investigation reveals that this primarily arises due to the concomitant increase in all the bond lengths in molecules. Such expansion of the molecules can be qualitatively understood on the basis of classical harmonic potentials in the bonds and loss of aromaticity in the rings. Specifically, the highly electronegative O and N elements in the base pairs anchor to Au atoms and form X-Au bonds, which leads to charge redistribution within the molecules. As a very important consequence of this, the nature of the hydrogen bonds (in Au(3)-A...T-Au(3) and in Au(3)-G...C-Au(3)) change substantially within these electrodes in comparison to gas-phase structures. These hydrogen bonds have a single-well potential energy profile (of the type N...H...O and N...H...N) instead of double-well potentials (like N-H...O or N-H...N/ N...H-N types). A detailed energy calculation along the proton movement pathway supports our conclusions.     

Copolymers containing alternating sequences of nucleic acid base pairs. II. Polymerization studies

The objective of this project was to utilize the alternating copolymerizability of electrondonor monomers with electron-acceptor monomers to selectively introduce nucleic acid bases into copolymers in a controlled sequence distribution. To this end, maleimide monomers containing the adenine, thymine, cytosine, and 6-chloropurine moieties were converted to their hompolymers. The homopolymer of 1-(vinyloxyethoxy)thymine was also prepared. Alternating copolymers of the adenine maleimide monomer and the corresponding 6-chloropurine maleimide monomer with 1-(vinyloxyethoxy)thymine were prepared. The latter copolymer was converted to the alternating adenine–thymine copolymer by reaction with ammonia. Characterization of the polymers and copolymers via spectroscopic methods and physical measurements confirm their proposed structures. Monomer syntheses and characterization, as well as studies designed to establish the extent and nature of adenine–thymine interactions in the copolymers, are reported in accompanying papers.     

Hydrogen-bonded nucleic acid base pairs containing unusual base tautomers: complete basis set calculations at the MP2 and CCSD(T) levels

The total interaction energies of altogether 15 hydrogen-bonded nucleic acid base pairs containing unusual base tautomers were calculated. The geometry properties of all selected adenine-thymine and guanine-cytosine hydrogen-bonded base pairs enable their incorporation into DNA. Unusual base pairing patterns were compared with Watson-Crick H-bonded structures of the adenine-thymine and guanine-cytosine pairs. The complete basis set (CBS) limit of the MP2 interaction energy and the CCSD(T) correction term, determined as the difference between the CCSD(T) and MP2 interaction energies, was evaluated. Extrapolation to the MP2 CBS limit was done using the aug-cc-pVDZ and aug-cc-pVTZ results, and the CCSD(T) correction term was determined with the 6-31G*(0.25) basis set. Final interaction energies were corrected while taking into account both tautomeric penalization determined at the CBS level and solvation/desolvation free energies. The situation for the adenine-thymine pairs is straightforward, and tautomeric pairs are significantly less stable than the Watson-Crick pair consisting of the canonical forms. In the case of the guanine-cytosine pair, the Watson-Crick structure made by canonical forms is again the most stable. The other two structures are, however, energetically rather similar (by 5 and 6 kcal/mol), which provides a very small but non-negligible chance of detecting these structures in the DNA double helix (1:5000). Due to the fact that DNA bases and base pairs incorporated into DNA are solvated less favorably than in isolated systems, this probability represents the very upper limit. The results clearly show how precisely the canonical building blocks of DNA molecules were chosen and how well their stability is maintained.     

Cu2+/+ cation coordination to adenine--thymine base pair. Effects on intermolecular proton-transfer processes

Intermolecular proton-transfer processes in the Watson & Crick adenine-thymine Cu+ and Cu2+ cationized base pairs have been studied using the density functional theory (DFT) methods. Cationized systems subject to study are those resulting from cation coordination to the main basic sites of the base pair, N7 and N3 of adenine and O2 of thymine. For Cu+ coordinated to N7 or N3 of adenine, only the double proton-transferred product is found to be stable, similarly to the neutral system. However, when Cu+ interacts with thymine, through the O2 carbonyl atom, the single proton transfer from thymine to adenine becomes thermodynamically spontaneous, and thus rare forms of the DNA bases may spontaneously appear. For Cu2+ cation, important effects on proton-transfer processes appear due to oxidation of the base pair, which stabilizes the different single proton-transfer products. Results for hydrated systems show that the presence of the water molecules interacting with the metal cation (and their mode of coordination) can strongly influence the ability of Cu2+ to induce oxidation on the base pair.     

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

Substituted Watson–Crick guanine–cytosine (GC) base pairs were recently shown to yield robust three‐state nanoswitches. Here, we address the question: Can such supramolecular switches also be based on Watson–Crick adenine‐thymine (AT) base pairs? We have theoretically analyzed AT pairs in which purine‐C8 and/or pyrimidine‐C6 positions carry a substituent X=NH^?, NH₂, NH₃⁺ (N series), O^?, OH or OH₂⁺ (O series), using the generalized gradient approximation (GGA) of density functional theory at the BP86/TZ2P level. Thus, we explore the trend in geometrical shape and hydrogen bond strengths in AT pairs along a series of stepwise protonations of the substituents. Introducing a charge on the substituents leads to substantial and characteristic changes in the individual hydrogen bond lengths when compared to the neutral AT pair. However, the trends along the series of negative, neutral, and positive substituents are less systematic and less pronounced than for GC. In certain instances, internal proton transfer from thymine to adenine occurs. Our results suggest that AT is a less suitable candidate than GC in the quest for chemically controlled nanoswitches.     

Multiple proton-transfer reactions in DNA base pairs by coordination of pt complex

The possibility of multiple proton-transfer reactions in DNA base pairs because of coordination of cisplatin is theoretically elucidated by density functional theory (DFT) and by quantum mechanics/molecular mechanics (QM/MM) methods with an ONIOM method. From the energetics of two base pairs with the cisplatin, it is theoretically confirmed that the Pt complex is likely to bind in the order cis-(CG)-Pt-(GC), cis-(CG)-Pt-(AT), cis-(TA)-Pt-(AT), where G, C, A, and T are guanine, cytosine, adenine, and thymine, respectively, and the Pt atom bonds to the N7 site of G and A. This result supports the experimental evidence, where the structure cis-A-Pt-A is seldom observed at room temperature. The single proton-transfer reaction occurs in one of the two GC pairs. No simultaneous single proton-transfer reaction can occur in both base pairs. Two different single proton-transferred structures (cis-(CG*)(d)-Pt-(GC)(p) and cis-(CG)(d)-Pt-(G*C)(p), where the asterisk means a proton donor of G) are as stable as the original structures (CG)(d)-Pt-(GC)(p). The same tendency was observed with cis-(CG*)-Pt-(AT). In contrast to cisplatin, multiple single proton-transfer reactions occur in the system consisting of two base pairs with transplatin. The optimized structure agrees with the experimental data for Pt-G coordination except for the hydrogen-bonding length.