Nuclear quantum effect on the hydrogen-bonded structure of guanine�Ccytosine pair

The structure of Watson?CCrick type guanine?Ccytosine (G?CC) base pair has been studied by classical hybrid Monte Carlo (HMC) and quantum path integral hybrid Monte Carlo (PIHMC) simulations on the semiempirical PM6 potential energy surface. For the three NH?X hydrogen-bonded moieties, the intramolecular NH bonds are found systematically longer while the H?X distance shorter in the PIHMC simulation than in the HMC simulation. We found that the hydrogen bonded length N?X correlates with the H?X distance, but not with the NH distance. A correlation is also between the neighboring hydrogen bonds in the G?CC base pair.     

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.     

A study of nucleic acid base-stacking by the Monte Carlo method: Extended cluster approach

The adenine-thymine (AT), adenine-uracil (AU) and guanine-cytosine (GC) base associates in clusters containing 400 water molecules were studied using a newly implemented Metropolis Monte Carlo algorithm based on the extended cluster approach. Starting from the hydrogen-bonded Watson-Crick geometries, all three base pairs are transformed into more favorable stacked configurations during the simulation. The obtained results show, for the first time, the transition from planar base pairs to stacked base associates in the Monte Carlo framework. Analysis of the interaction energies shows that, in the water cluster, the stacked dimers are energetically preferable compared to the corresponding Watson-Crick base pairs. This is due to the larger base-water interaction in the stacked structures. The water-water interaction is one of the main factors promoting the formation of stacked dimers, and the obtained data confirm the crucial role of the water-water interactions in base stacking.     

At nonzero temperatures, stacked structures of methylated nucleic acid base pairs and microhydrated nonmethylated nucleic acid base pairs are favored over planar hydrogen-bonded structures: a molecular dynamics simulations study

The dynamic structure of all ten possible nucleic acid (NA) base pairs and methylated NA base pairs hydrated by a small number of water molecules (from 1 to 16) was determined by using molecular dynamics simulations in the NVE microcanonical and NVT canonical ensembles with the Cornell force field (W. D. Cornell, P. Cieplak, C. I. Bayly, I. R. Gould, K. M. Merz, D. M. Ferguson, D. C. Spellmeyer, T. Fox, J. E. Caldwell, P. Kollman, J. Am. Chem. Soc. 1995, 117, 5179). The presence of one water molecule does not affect the structure of any hydrogen-bonded (H-bonded) nonmethylated base pair. An equal population of H-bonded and stacked structures of adenine...adenine, adenine...guanine and adenine... thymine pairs is reached if as few as two water molecules are present, while obtaining equal populations of these structures in the case of adenine...cytosine, cytosine...thymine, guanine... guanine and guanine...thymine required the presence of four water molecules, and in the case of guanine...cytosine, six. A comparable population of planar, H-bonded and stacked structures for cytosine...cytosine and thymine... thymine base pairs was only obtained if at least eight water molecules hydrated a pair. Methylation of bases changed the situation dramatically and stacked structures were favoured over H-bonded ones even in the absence of water molecules in most cases. Only in the case of methyl cytosine...methyl cytosine, methyl guanine...methyl guanine and methyl guanine...methyl cytosine pairs were two, two or six water molecules, respectively, needed in order to obtain a comparable population of planar, H-bonded and stacked structures. We believe that these results give clear evidence that the preferred stacked structure of NA base pairs in the microhydrated environment, and also apparently in a regular solvent, is due to the hydrophilic interaction of a small number of water molecules. In the case of methylated bases, it is also due to the fact that the hydrogen atoms most suitable for the formation of H-bonds have been replaced by a methyl group. A preferred stacked structure is, thus, not due to a hydrophobic interaction between a large bulk of water molecules and the base pair, as believed.     

Nuclear quantum effect and temperature dependency on the hydrogen‐bonded structure of base pairs

The structure of Watson–Crick‐type adenine‐thymine and guanine‐cytosine pairs has been studied by hybrid Monte Carlo (HMC) and path integral hybrid Monte Carlo (PIHMC) simulations with the use of semiempirical PM6‐DH+ method in the gas phase. We elucidated the nuclear quantum effect and temperature dependency on the hydrogen‐bonded moiety of base pairs. It was shown that the contribution of nuclear quantum effect on the hydrogen‐bonded structure is significant not only at low temperature 150 K but also at temperature as high as 450 K. The relative position of hydrogen‐bonded proton between two heavy atoms and the nuclear quantum nature of the proton are also shown. Furthermore, we have applied principal component analysis to HMC and PIHMC simulations to analyze the nuclear quantum effect on intermolecular motions. We found that the ratio of Buckle mode (lowest vibrational mode from normal mode analysis) decreases due to the nuclear quantum effect, whereas that of Propeller mode (second lowest vibrational mode) increases. In addition, nonplanar structures of base pairs were found to become stable due to the nuclear quantum effect from two‐dimensional free energy landscape along Buckle and Propeller modes. © 2013 Wiley Periodicals, Inc.     

Toward the exact solution of the electronic Schrödinger equation for noncovalent molecular interactions: worldwide distributed quantum monte carlo calculations

Quantum Monte Carlo (QMC) calculations on the stacked (st) and Watson/Crick (wc) bound adenine/thymine (A/T) and cytosine/guanine (C/G) DNA base pair complexes were made possible with the first large scale distributed computing project in ab initio quantum chemistry, Quantum Monte Carlo at Home (QMC@HOME). The results for the interaction energies (wc-A/T = 15.7 kcal/mol, wc-C/G = 30.2 kcal/mol, st-A/T = 13.1 kcal/mol, st-C/G = 19.6 kcal/mol) are in very good agreement with the best known coupled-cluster based estimates. The accuracy of these values is further supported by calculations on the S22 benchmark set of noncovalently bound systems, for which we obtain a small mean absolute deviation of 0.68 kcal/mol. Our results support previous claims that the stacking energies are of comparable magnitude to the interactions of the commonly discussed hydrogen-bonded motif. Furthermore, we show that QMC can serve as an advantageous alternative to conventional wave function methods for large noncovalently bound systems. We also investigated in detail all technical parameters of the QMC simulations and recommend a careful optimization procedure of the Jastrow correlation factors in order to obtain numerically stable and reliable results.     

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.     

Mismatch base pairing of the mutagen 8‐oxoguanine and its derivatives with adenine: A theoretical search for possible antimutagenic agents

Molecular geometries of 8‐oxoguanine (8OG), those of its substituted derivatives with the substitutions CH₂, CF₂, CO, CNH, O, and S in place of the N7H7 group, adenine (A), and the base pairs of 8OG and its substituted derivatives with adenine were optimized using the RHF/6‐31+G* and B3LYP/6‐31+G* methods in gas phase. All the molecules and their hydrogen‐bonded complexes were solvated in aqueous media employing the polarized continuum model (PCM) of the self‐consistent reaction field (SCRF) theory using the RHF/6‐31+G* and B3LYP/6‐31+G* methods. The optimized geometrical parameters of the 8OG‐A base pair at the RHF/6‐31+G* and B3LYP/6‐31+G* levels of theory agree satisfactorily with those of an oligonucleotide containing the base pair found from X‐ray crystallography. The pattern of hydrogen bonding in the CF₂‐ and O‐substituted 8OG‐A base pair is of Watson–Crick type and that in the unsubstituted and CH₂‐, CNH‐, and S‐substituted base pairs is of Hoogsteen type. In the CO‐substituted base pair, the hydrogen bonding pattern is of neither Watson–Crick nor Hoogsteen type. The CF₂‐substitution appears to introduce steric hindrance for stacking of DNA bases. On the basis of these results, it appears that among all the substituted 8OG molecules considered here, the O‐substituted derivative may be useful as an antimutagenic drug. It is, however, subject to experimental verification. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Switching binding sites: low-temperature NMR studies on adenosine-aspartic acid interactions

Low-temperature NMR experiments were performed on mixtures of adenine nucleosides and aspartic acid derivatives in a freonic solvent. By acquiring spectra at temperatures as low as 123 K, the regime of slow hydrogen bond exchange is reached and hydrogen-bonded complexes can be characterized in detail. With 2'-deoxyadenosine lacking a 2'-substituent, N-Boc-protected aspartic acid benzyl ester binds through its carboxylic acid side chain to the Watson-Crick site of the adenine base, forming a strong hydrogen bond with the proton located close to the center between the oxygen donor and adenine N1 nitrogen acceptor. However, in the case of 2'-O-silylated adenine ribofuranosides, noncovalent interactions of the 2'-substituent with protecting groups on the amino acid shift the binding mode toward a Hoogsteen geometry with only a moderately strong hydrogen bond involving adenine N7.     

Infrared Spectroscopic Observation of a G–C+ Hoogsteen Base Pair in the DNA:TATA‐Box Binding Protein Complex Under Solution Conditions

Hoogsteen DNA base pairs (bps) are an alternative base pairing to canonical Watson–Crick bps and are thought to play important biochemical roles. Hoogsteen bps have been reported in a handful of X‐ray structures of protein–DNA complexes. However, there are several examples of Hoogsteen bps in crystal structures that form Watson–Crick bps when examined under solution conditions. Furthermore, Hoogsteen bps can sometimes be difficult to resolve in DNA:protein complexes by X‐ray crystallography due to ambiguous electron density and by solution‐state NMR spectroscopy due to size limitations. Here, using infrared spectroscopy, we report the first direct solution‐state observation of a Hoogsteen (G–C⁺) bp in a DNA:protein complex under solution conditions with specific application to DNA‐bound TATA‐box binding protein. These results support a previous assignment of a G–C⁺ Hoogsteen bp in the complex, and indicate that Hoogsteen bps do indeed exist under solution conditions in DNA:protein complexes.     

Contiguous metal-mediated base pairs comprising two Ag(I) ions

The incorporation of transition‐metal ions into nucleic acids by using metal‐mediated base pairs has proved to be a promising strategy for the site‐specific functionalization of these biomolecules. We report herein the formation of Ag⁺‐mediated Hoogsteen‐type base pairs comprising 1,3‐dideaza‐2′‐deoxyadenosine and thymidine. By defunctionalizing the Watson–Crick edge of adenine, the formation of regular base pairs is prohibited. The additional substitution of the N3 nitrogen atom of adenine by a methine moiety increases the basicity of the exocyclic amino group. Hence, 1,3‐dideazaadenine and thymine are able to incorporate two Ag⁺ ions into their Hoogsteen‐type base pair (as compared with one Ag⁺ ion in base pairs with 1‐deazaadenine and thymine). We show by using a combination of experimental techniques (UV and circular dichroism (CD) spectroscopies, dynamic light scattering, and mass spectrometry) that this type of base pair is compatible with different sequence contexts and can be used contiguously in DNA double helices. The most stable duplexes were observed when using a sequence containing alternating purine and pyrimidine nucleosides. Dispersion‐corrected density functional theory calculations have been performed to provide insight into the structure, formation and stabilization of the twofold metalated base pair. They revealed that the metal ions within a base pair are separated by an Ag???Ag distance of about 2.88 Å. The Ag–Ag interaction contributes some 16 kcal mol^?1 to the overall stability of the doubly metal‐mediated base pair, with the dominant contribution to the Ag–Ag bonding resulting from a donor–acceptor interaction between silver 4d‐type and 4s orbitals. These Hoogsteen‐type base pairs enable a higher functionalization of nucleic acids with metal ions than previously reported metal‐mediated base pairs, thereby increasing the potential of DNA‐based nanotechnology.     

Stabilization by extra-helical thymines of a DNA duplex with Hoogsteen base pairs

We present the crystal structure of the DNA duplex formed by d(ATATATCT). The crystals contain seven stacked antiparallel duplexes in the asymmetric unit with A.T Hoogsteen base pairs. The terminal CT sequences bend over so that the thymines enter the minor groove and form a hydrogen bond with thymine 2 of the complementary strand in the Hoogsteen duplex. Cytosines occupy extra-helical positions; they contribute to the crystal lattice through various kinds of interactions, including a unique CAA triplet. The presence of thymine in the minor groove apparently contributes to the stability of the DNA duplex in the Hoogsteen conformation. These observations open the way toward finding under what conditions the Hoogsteen duplex may be stabilized in vivo. The present crystal structure also confirms the tendency of A.T-rich oligonucleotides to crystallize as long helical stacks of duplexes.     

Potential energy surfaces of an adenine-thymine base pair and its methylated analogue in the presence of one and two water molecules: molecular mechanics and correlated ab initio study

Potential energy surfaces of monohydrated and dihydrated adenine-thymine and 9-methyladenine-1-methylthymine base pairs were examined by the molecular dynamics/quenching technique using the Cornell et al. force field (J. Am. Chem. Soc. 1995, 117, 5179). Long runs of molecular dynamics/quenching calculations allowed us to evaluate the free energy surface. The most stable and populated structures found were fully reoptimized at the correlated ab initio level employing the resolution of identity M?ller-Plesset method. A systematic study of the base pairs' microhydration using both the empirical and the high-level correlated ab initio approaches is presented for the first time. We show that the occurrence of water molecules and their gradually increasing number as well as the methylation of the bases favor stacked structures over the planar hydrogen-bonded ones. These results based on the correlated ab initio calculations are in the excellent agreement with data obtained from our previous empirical potential molecular dynamics study (Kabelác et al. Chem.-Eur. J. 2001, 7, 2067).     

Covalently Cross-Linked Watson–Crick Base Pair Models

Structural characteristics of Watson–Crick hydrogen-bonded base pairs are displayed by methylene-bridged base pairs of type A . The shown superposition of the X-ray structure obtained for the base pair A (Rib¹=Et; Rib²=Me) over that of a C–G base pair illustrates that A occupies an area similar to that occupied by a traditional Watson–Crick hydrogen-bonded base pair. Temperature-dependent ¹H NMR studies indicate that the energy barrier for rotation along its CH₂ bridge is about 10 kcal mol⁻¹, and that it exists predominantly in one conformer at −70°C.     

Relative log P and solution structure for small organic solutes in the chloroform/water system using monte carlo methods

Relative log P values of dimethylether to methanol and dimethylamine to methylamine were calculated in the chloroform/water system using Monte Carlo simulations and statistical perturbation theory. Correct ordering of the calculated relative log Ps was obtained for the two pairs although the method leads to an overestimation of these values. In aqueous solution, both dimethyl ether and dimethylamine solutes are proton acceptors forming a single hydrogen bond to water. Dimethylamine forms a stable N? H-Ow hydrogen bond while the water hydrogen is poorly localized in the O? H-Ow bond to the ether. In chloroform, the solvent molecules are less ordered around the solutes than was found around methanol and methylamine.     

Heterochiral DNA with Complementary Strands with α-d and β-d Configurations: Hydrogen-Bonded and Silver-Mediated Base Pairs with Impact of 7-Deazapurines Replacing Purines

Heterochiral DNA with hydrogen-bonded and silver-mediated base pairs have been constructed using complementary strands with nucleosides with α-d or β-d configuration. Anomeric phosphoramidites were employed to assemble the oligonucleotides. According to the T_m values and thermodynamic data, the duplex stability of the heterochiral duplexes was similar to that of homochiral DNA, but mismatch discrimination was better in heterochiral DNA. Replacement of purines by 7-deazapurines resulted in stable parallel duplexes, thereby confirming Watson–Crick-type base pairing. When cytosine was facing cytosine, thymine or adenine residues, duplex DNA formed silver-mediated base pairs in the presence of silver ions. Although the CD spectra of single strands with α-d configuration display mirror-like shapes to those with the β-d configuration, the CD spectra of the hydrogen-bonded duplexes and those with a limited number of silver pairs show a B-type double helix almost indistinguishable from natural DNA. Nonmelting silver ion–DNA complexes with entirely different CD spectra were generated when the number of silver ions was equal to the number of base pairs.     

Does the A.T or G.C base-pair possess enhanced stability? Quantifying the effects of CH...O interactions and secondary interactions on base-pair stability using a phenomenological analysis and ab initio calculations

An empirically based relationship between overall complex stability (-DeltaG degrees ) and various possible component interactions is developed to probe the question of whether the A.T/U and G.C base-pairs exhibit enhanced stability relative to similarly hydrogen-bonded complexes. This phenomenological approach suggests ca. 2-2.5 kcal mol-1 in additional stability for A.T owing to a group interaction containing a CH...O contact. Pairing geometry and the role of the CH...O interaction in the A.T base-pair were also probed using MP2/6-31+G(d,p) calculations and a double mutant cycle. The ab initio studies indicated that Hoogsteen geometry is preferred over Watson-Crick geometry in A.T by ca. 1 kcal mol-1. Factors that might contribute to the preference for Hoogsteen geometry are a shorter CH...O contact, a favorable alignment of dipoles, and greater distances between secondary repulsive sites. The CH...O interaction was also investigated in model complexes of adenine with ketene and isocyanic acid. The ab initio calculations support the result of the phenomenological approach that the A.T base-pair does have enhanced stability relative to hydrogen-bonded complexes with just N-H...N and N-H...O hydrogen bonds.     

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.     

Hydration regulates thermodynamics of G-quadruplex formation under molecular crowding conditions

The effect of molecular crowding on the structure and stability of biomolecules has become a subject of increasing interest because it can clarify how biomolecules behave under cell-mimicking conditions. Here, we quantitatively analyzed the effects of molecular crowding on the thermodynamics of antiparallel G-quadruplex formation via Hoogsteen base pairs and of antiparallel hairpin-looped duplex (HP duplex) formation via Watson-Crick base pairs. The free energy change at 25 degrees C for G-quadruplex formation decreased from -3.5 to -5.5 kcal mol(-1) when the concentration of poly(ethylene glycol) 200 was increased from 0 to 40 wt %, whereas that of duplex formation increased from -9.8 to -6.9 kcal mol(-1). These results showed that the antiparallel G-quadruplex is stabilized under molecular crowding conditions, but that the HP duplex is destabilized. Moreover, plots of stability (ln K(obs)) of the DNA structures versus water activity (ln a(w)) demonstrated that the ln K(obs) for G-quadruplex formation decreased linearly as the ln a(w) increased, whereas that for duplex formation increased linearly with the increase in ln a(w), suggesting that the slope approximately equals the number of water molecules released or taken up during the formation of these structures. Thus, molecular crowding affects the thermodynamics of DNA structure formation by altering the hydration of the DNA. The stabilization of the DNA structures with Hoogsteen base pairs and destabilization of DNA structures with Watson-Crick base pairs under molecular crowding conditions lead to structural polymorphism of DNA sequences regulated by the state of hydration.