Structure, stability, and dynamics of canonical and noncanonical base pairs: quantum chemical studies

The importance of non-Watson-Crick base pairs in the three-dimensional structure of RNA is now well established. The structure and stability of these noncanonical base pairs are, however, poorly understood. We have attempted to understand structural features of 33 frequently occurring base pairs using density functional theory. These are of three types, namely (i) those stabilized by two or more polar hydrogen bonds between the bases, (ii) those having one polar and another C-H...O/N type interactions, and (iii) those having one H-bond between the bases and another involving one of the sugars linked to the bases. We found that the base pairs having two polar H-bonds are very stable as compared to those having one C-H...O/N interaction. Our quantitatively analysis of structures of these optimized base pairs indicates that they possess a different amount of nonplanarity with large propeller or buckle values as also observed in the crystal structures. We further found that geometry optimization does not modify the hydrogen-bonding pattern, as values of shear and open angle of the base pairs remain conserved. The structures of initial crystal geometry and final optimized geometry of some base pairs having only one polar H-bond and a C-H...O/N interaction, however, are significantly different, indicating the weak nature of the nonpolar interaction. The base pair flexibility, as measured from normal-mode analysis, in terms of the intrinsic standard deviations of the base pair structural parameters are in conformity with those calculated from RNA crystal structures. We also noticed that deformation of a base pair along the stretch direction is impossible for all of the base pairs, and movements of the base pairs along shear and open are also quite restricted. The base pair opening mode through alteration of propeller or buckle is considerably less restricted for most of the base pairs.     

Conformational specificity of non-canonical base pairs and higher order structures in nucleic acids: crystal structure database analysis

Non-canonical base pairs contribute immensely to the structural and functional variability of RNA, which calls for a detailed characterization of their spatial conformation. Intra-base pair parameters, namely propeller, buckle, open-angle, stagger, shear and stretch describe structure of base pairs indicating planarity and proximity of association between the two bases. In order to study the conformational specificities of non-canonical base pairs occurring in RNA crystal structures, we have upgraded NUPARM software to calculate these intra-base pair parameters using a new base pairing edge specific axis system. Analysis of base pairs and base triples with the new edge specific axis system indicate the presence of specific structural signatures for different classes of non-canonical pairs and triples. Differentiating features could be identified for pairs in cis or trans orientation, as well as those involving sugar edges or C-H-mediated hydrogen bonds. It was seen that propeller for all types of base pairs in cis orientation are generally negative, while those for trans base pairs do not have any preference. Formation of a base triple is seen to reduce propeller of the associated base pair along with reduction of overall flexibility of the pairs. We noticed that base pairs involving sugar edge are generally more non-planar, with large propeller or buckle values, presumably to avoid steric clash between the bulky sugar moieties. These specific conformational signatures often provide an insight into their role in the structural and functional context of RNA.     

Vibrational dynamics of DNA. I. Vibrational basis modes and couplings

Carrying out density functional theory calculations of four DNA bases, base derivatives, Watson-Crick (WC) base pairs, and multiple-layer base pair stacks, we studied vibrational dynamics of delocalized modes with frequency ranging from 1400 to 1800 cm(-1). These modes have been found to be highly sensitive to structure fluctuation and base pair conformation of DNA. By identifying eight fundamental basis modes, it is shown that the normal modes of base pairs and multilayer base pair stacks can be described by linear combinations of these vibrational basis modes. By using the Hessian matrix reconstruction method, vibrational coupling constants between the basis modes are determined for WC base pairs and multilayer systems and are found to be most strongly affected by the hydrogen bonding interaction between bases. It is also found that the propeller twist and buckle motions do not strongly affect vibrational couplings and basis mode frequencies. Numerically simulated IR spectra of guanine-cytosine and adenine-thymine bases pairs as well as of multilayer base pair stacks are presented and described in terms of coupled basis modes. It turns out that, due to the small interlayer base-base vibrational interactions, the IR absorption spectrum of multilayer base pair system does not strongly depend on the number of base pairs.     

NMR imino proton exchange experiments on duplex DNA primarily monitor the opening of purine bases

Molecular dynamics simulations were performed to investigate GC and AT base opening events in DNA. Calculated equilibrium constants between the base open (or flipped) and closed states were shown to be in good agreement with experimental data from NMR imino proton exchange experiments. Analysis of the computed results indicates that the equilbrium constants are dominated by the opening of the A and G bases in the AT and GC base pairs, respectively. Thus, the present results predict that NMR imino proton exchange experiments of base opening are primarily monitoring the opening of purine bases.     

Double-proton transfer in adenine-thymine and guanine-cytosine base pairs. A post-Hartree-Fock ab initio study

The results of a comprehensive study on the double-proton transfer in Adenine-Thymine (AT) and Guanine-Cytosine (GC) base pairs at room temperature in gas phase and with the inclusion of environmental effects are obtained. The double-proton-transfer process has been investigated in the AT and GC base pairs at the B3LYP/6-31G(d) and MP2/6-31G(d) levels of theory. It has been predicted that the hydrogen-bonded bases possess nonplanar geometries due to sp3 hybridization of nitrogen atoms and because of the soft intermolecular vibrations in the molecular complexes. An analysis of the energetic parameters of the local minima suggests that rare AT base pair conformation is not populated due to the shallowness of this minimum, which completely disappears from the Gibbs free energy surface. The stabilization of canonic or rare forms of the DNA bases by water molecules and metal cations has been predicted by calculating the optimal configuration of charges (using differential product/transition state stabilization approach) followed by calculations of the interactions between the base pair and a water/sodium cation.     

Revisiting the effects of sequence and structure on the hydrogen bonding and π-stacking interactions in nucleic acids

Calculated electron densities from PBE0/6-31+G(d,p) were analyzed with respect to the hydrogen bonding within a nucleic acid base pair and the π-stacking between sets of base pairs. From published X-ray crystallographic data, base pairs were isolated from a total of 11 DNA and RNA duplexes, and their experimental geometry was maintained throughout the analyses. Focusing solely on Watson-Crick base pairs, from the values of the electron density between interacting nuclei (at the bond critical points), we provide quantitative data on individual weak interactions. For hydrogen bonding, in addition to quantifying the scissoring effect in GC base pairs, the origin of the controversy around the relative stability of AT and AU base pairs is identified and resolved. Thus, it is illustrated how the conclusion as to their relative stability rests on the specific choice of oligonucleotides compared. For π-stacking, sequence effects for tandem AT base pairs are captured, quantified, and explained, and the greater sensitivity of GC, over AT, sequences to the rise parameter is established. The results presented here show that, from experimental geometries and their electron densities, previously determined effects of the sequence and structure of a duplex on the stabilizing interactions can be captured, quantified, and traced back to the geometry of the base pairs.     

Hydrogen bond vibrations of 2-aminopyridine.2-pyridone,a Watson-Crick analogue of adenine.uracil

The 2-aminopyridine.2-pyridone complex (2AP.2PY) is linked by antiparallel N-H.O=C and N-H.N hydrogen bonds, providing a model for the Watson-Crick hydrogen bond configuration of the adenine-thymine and adenine-uracil nucleobase pairs. Mass-selected S(1) <--> S(0) vibronic spectra of the supersonically cooled 2AP.2PY base pair analogue were measured by laser resonant two-photon ionization and emission spectroscopies. The hydrogen bond vibrations nu(2) (buckle, out-of-plane) and the three in-plane vibrations nu(3) (opening), nu(5) (shear), and nu(6) (stretch) were observed in the S(0) and S(1) states, giving detailed information on the stretching and deformation force constants of the (amide)N-H.N(pyridine) and the (amino)N-H.O=C hydrogen bonds. Density functional calculations with the B3LYP functional and the 6-311++G(d,p) and 6-311++G(2d,2p) basis sets yield ground-state hydrogen bond frequencies in close agreement with experiment.     

水环境对GC和AT碱基对质子转移的影响

采用B3LYP/DZP++的方法研究了第一水化层作用和连续化处理的水溶剂作用对鸟嘌呤-胞嘧啶(GC)碱基对和腺嘌呤-胸腺嘧啶(AT)碱基对质子转移反应的影响. GC和AT碱基对在连续化水溶剂作用下,均发生单质子转移(SPT1)和分步的双质子转移(DPT),而在第一水化层5 个水分子的作用下(GC·5H₂O,AT·5H₂O)或同时考虑第一水化层作用和连续化水溶剂作用(GC·5H₂O+PCM,AT·5H₂O+PCM)时,GC和AT碱基对的质子转移均只得到单质子转移反应(SPT1). 单质子转移过程中的活化能变化情况表明：第一水化层对GC和AT碱基对结构和质子转移影响较大,水环境对碱基对的作用主要发生在第一水化层.     

Reorganization energy, activation energy, and mechanism of hole transfer process in DNA: a theoretical study

The density functional calculations with aug-cc-pVDZ basis sets on cationic guanine-cytosine (GC(+)) and adenine-thymine (AT(+)) base pairs suggest that the cationic charge is almost entirely localized on the G and A units with significant changes in the N-H and N...O distances around the H-bonded area. While the calculated intramolecular reorganization energy (lambda(v)) for a GC base pair (0.75 eV) is remarkably larger than that for an isolated G base (0.49 eV), for the AT base pairs these values (0.44 and 0.40 eV) are almost the same. The gas phase activation energies (E(a)) for GC(+)GC-->GCGC(+), AT(+)AT-->ATAT(+), and GC(+)AT-->GCAT(+) hole transfer processes are 0.19, 0.11, and 0.73 eV with rate constants of 1.69 x 10(11), 3.15 x 10(11), and 4.61(0.168) s(-1), respectively, at 298 K. An alternative mechanism of hole transfer has been proposed on the basis of energy barriers.     

Elementary steps for charge transport in DNA: thermal activation vs. tunneling

Using stacks of Watson–Crick base pairs as an important example of multichromophoric molecular assemblies, we studied charge migration in DNA with special emphasis on the mechanism of hole hopping between neighboring guanines (G) connected by the adenine–thymine (AT) bridge. The tight-binding model proposed for this elementary step shows that for short AT bridges, hole transfer between two G bases proceeds via quantum mechanical tunneling. By contrast, hopping over long bridges requires thermal activation. The condition for crossover between tunneling and thermal activation near room temperature is specified and applies to the analysis of experimental data. We show that thermal activation dominates, if the bridge between two G bases contains more than three AT pairs. Our theoretical findings predict that the replacement of AT base pairs by GC pairs increases the efficiency of hole transport only in the case of short base pair sequences. For long sequences, however, the opposite effect is expected.     

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.     

Remarkable metal counterion effect on the internucleotide J-couplings and chemical shifts of the N-H...N hydrogen bonds in the W-C base pairs

The effects of metal ion binding on the (2h) J(NN)-coupling and delta( (1)H)/Deltadelta( (15)N) chemical shifts of N-H...N H-bond units in internucleotide base pairs were explored by a combination of density functional theory calculations and molecular dynamics (MD) simulations. Results indicate that the NMR parameters vary considerably upon cation binding to the natural GC or AT base pairs, and thus can be used to identify the status of the base pairs, if cation-perturbed. The basic trend is that cation perturbation causes (2h) J(NN) to increase, Deltadelta( (15)N) to decrease, and delta( (1)H) to shift upfield for GC, and in the opposite directions for AT. The magnitudes of variation are closely related to the Lewis acidity of the metal ions. For both base pair series (M(z+)GC and M(z+)AT), these NMR parameters are linearly correlated among themselves. Their values depend strongly on the energy gaps (Delta(ELP-->sigma*)) and the second-order interaction energies ( E(2)) between the donor N lone pair (LP(N)) and the acceptor sigma* N-H localized NBO orbitals. In addition, the (2h) J NN changes are also sensitive to the amount of sigma charge transfer from LP(N) to sigma*(N-H) NBOs or from the purine to the pyrimidine moieties. The different trends are a consequence of the different H-bond patterns combined with the polarization effect of the metal ions in the cationized M(z+)AT series, M(z+) <-- A --> T, and the cationized GC series, M(z+) <-- G <-- C. The predicted cation-induced systematic trends of (2h) J(NN) and delta( (15)N, (1)H) in N-H...N H-bond units may provide a new approach to the determination of H-bond structure and strength in Watson-Crick base pairs, and provide an alternative probe of the heterogeneity of DNA sequences.     

Nucleic acid base pair dynamics: the impact of sequence and structure using free-energy calculations

Molecular dynamics free-energy calculations of base pair opening within double helical DNA and RNA are used to explain why A-tracts (oligo-adenine repeats) greatly increase the lifetimes of AT base pairs, whereas the structural and the chemical changes involved in passing from B-DNA to A-RNA have comparatively small effects.     

铂配合物与DNA碱基对间相互作用的理论研究

用量子化学方法研究一系列Pt(II)配合物作用于嘌呤碱基N7位点后对Watson-Crick碱基对AT、GC的影响. 计算结果显示铂配体与碱基对AT、GC的作用以静电作用为主，同时极化作用也是影响GC碱基对的重要因素. 静电作用极大地增强了铂化嘌呤碱基与嘧啶碱基间的相互作用，而嘌呤碱基与嘧啶碱基间作用与未铂化碱基对作用相近. Pd(II) 和 Ni(II)的相关研究得到类似的结果. 碱基对间氢键作用“二阶微扰能”分析结果与氢键强弱变化一致.     

Influence of magnesium ions on spontaneous opening of DNA base pairs

A large amount of experimental evidence is available for the effects of magnesium ions on the structure and the stability of the DNA double helix. Less is known, however, on how these ions affect the dynamics of the molecule and the stability of each individual base pair. The present work addresses these questions by a study of the DNA duplex [dCGCAGATCTGCG]2, and its interactions with magnesium ions using nuclear magnetic resonance (NMR) spectroscopy and proton exchange. Two-dimensional NMR experiments indicate that binding of magnesium to this DNA duplex does not affect its structure. However, even in the absence of structural changes, magnesium ions specifically affect the exchange properties of imino protons in the four GC/CG base pairs that are located in the interior of the double helix. These specific changes do not result from alterations in the rates of spontaneous opening of these base pairs. Instead, the changes most likely reflect an enhancement in the energetic propensity for spontaneous opening of the GC/CG base pairs that is induced by the binding of magnesium ions.     

The Protonated Guanine–Cytosine Base Pair

Protonated base pairs were recently implicated in the context of DNA proton transfer and charge migration. The effects of protonating different sites of the guanine–cytosine (GC) base pair are studied here by using the DZP++ B3LYP density functional method. Optimized structures for the protonated GC base pair are compared with those of parent GC and the neutral hydrogenated GC radical (GCH). Proton and hydrogen‐atom additions significantly disturb the structure of the GC base pair. However, the structural perturbations arising from protonation are often less than those arising from hydrogenation of GC. Protonation of the GC base pair causes significant strengthening of the interstrand hydrogen bonds and a concomitant increase in the base dissociation energies. The adiabatic ionization potentials (AIPs), vertical ionization potentials (VIPs), and proton affinities (PAs) for the different protonation sites of the GC base pair are predicted. The N7 site of guanine is the preferred site for protonation of the GC base pair.     

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.     

Intercalators. 1. Nature of stacking interactions between intercalators (ethidium, daunomycin, ellipticine, and 4',6-diaminide-2-phenylindole) and DNA base pairs. Ab initio quantum chemical, density functional theory, and empirical potential study

Properties of isolated intercalators (ethidium (E), daunomycin (D), ellipticine (EL), and 4,6'-diaminide-2-phenylindole (DAPI)) and their stacking interactions with adenine...thymine (AT) and guanine...cytosine (GC) nucleic acid base pairs were investigated by means of a nonempirical correlated ab initio method. All intercalators exhibit large charge delocalization, and none of them (including the DAPI dication) exhibits a site with dominant charge. All intercalators have large polarizability and are good electron acceptors, while base pairs are good electron donors. MP2/6-31G*(0.25) stabilization energies of intercalator...base pair complexes are large (E...AT, 22.4 kcal/mol; D...GC, 17.8 kcal/mol; EL...GC, 18.2 kcal/mol; DAPI...GC, 21.1 kcal/mol) and are well reproduced by modified AMBER potential (van der Waals radii of intercalator atoms are enlarged and their energy depths are increased). Standard AMBER potential underestimates binding, especially for DAPI-containing complexes. Because the DAPI dication is the best electron acceptor (among all intercalators studied), this difference is explained by the importance of the charge-transfer term, which is not included in the AMBER potential. For the neutral EL molecule, the standard AMBER force field provides correct results. The Hartree-Fock and DFT/B3LYP methods, not covering the dispersion energy, fail completely to reveal any energy minimum at the potential energy curve of the E...AT complex, and these methods thus cannot be recommended for a study of intercalation process. On the other hand, an approximate version of the DFT method, which was extended to cover London dispersion energy, yields for all complexes very good stabilization energies that are well comparable with referenced ab initio data. Besides the vertical dependence of the interaction, an energy twist dependence of the interaction energy was also investigated by a reference correlated ab initio method and empirical potentials. It is concluded that, despite the cationic (E +1, D +1, DAPI +2) or polar (EL) character of the intercalators investigated, it is the dispersion energy which predominantly contributes to the stability of intercalator...base pair complexes. Any procedure which does not cover dispersion energy is thus not suitable for studying the process of intercalation.     

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.