Induced-fit recognition of DNA by small circular oligonucleotides

We have investigated the molecular interaction between cyclic and linear oligonucleotides. We have found that short cyclic oligonucleotides can induce hairpinlike structures in linear DNA fragments. By using NMR and CD spectroscopy we have studied the interaction of the cyclic oligonucleotide d with d, as well as with its two linear analogs d(GTCCCTCA) and d(CTCAGTCC). Here we report the NMR structural study of these complexes. Recognition between these oligonucleotides occurs through formation of four intermolecular Watson-Crick base pairs. The three-dimensional structure is stabilized by two tetrads, formed by facing the minor-groove side of the Watson-Crick base pairs. Overall, the structure is similar to those observed previously in other quadruplexes formed by minor-groove alignment of Watson-Crick base pairs. However, in this case the complexes are heterodimeric and are formed by two different tetrads (G:C:A:T and G:C:G:C). These complexes represent a new model of DNA recognition by small cyclic oligonucleotides, increasing the number of potential applications of these interesting molecules.     

Binding of gold clusters with DNA base pairs: a density functional study of neutral and anionic GC-Aun and AT-Aun (n = 4, 8) complexes

Binding of clusters of gold atoms (Au) with the guanine-cytosine (GC) and adenine-thymine (AT) Watson-Crick DNA base pairs was studied using the density functional theory (DFT). Geometries of the neutral GC-Au(n) and AT-Au(n) and the corresponding anionic (GC-Au(n))(-1) and (AT-Au(n))(-1) (n = 4, 8) complexes were fully optimized in different electronic states, that is, singlet and triplet states for the neutral complexes and doublet and quartet states for the anionic complexes, using the B3LYP density functional method. The 6-31+G basis set was used for all atoms except gold. For gold atoms, the Los Alamos effective core potential (ECP) basis set LanL2DZ was employed. Vibrational frequency calculations were performed to ensure that the optimized structures corresponded to potential energy surface minima. The gold clusters around the neutral GC and AT base pairs have a T-shaped structure, which satisfactorily resemble those observed experimentally and in other theoretical studies. However, in anionic GC and AT base pairs, the gold clusters have extended zigzag and T-shaped structures. We found that guanine and adenine have high affinity for Au clusters, with their N3 and N7 sites being preferentially involved in binding with the same. The calculated adiabatic electron affinities (AEAs) of the GC-Au(n)complexes (n = 4, 8) were found to be much larger than those of the isolated base pairs.     

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.     

Ultrafast deactivation processes in aminopyridine clusters: excitation energy dependence and isotope effects

Fast excited-state relaxation in H-bonded aminopyridine clusters occurs via hydrogen transfer in the excited state. We used femtosecond pump-probe spectroscopy to characterize the excited-state reaction coordinate. Considerable isotope effects for partially deuterated clusters indicate that H-transfer is the rate-limiting step and validate ab initio calculations in the literature. A nonmonotonous dependence on the excitation energy, however, disagrees with the picture of a simple barrier along the reaction coordinate. An aminopyridine dimer serves as a model for Watson-Crick base pairs, where similar reactions have been predicted by theory.     

Structural analysis of metal interactions with the dinucleotide duplex, dCG x dCG, using ion mobility mass spectrometry

The metal binding properties of the dinucleotide duplex, dCG x dCG, were analyzed in the gas phase with ion mobility mass spectrometry. Both MALDI and ESI were used to generate [M(dCG x dCG)]+ complexes. The collision cross section of each complex was measured in helium using ion mobility based methods and compared to calculated cross sections of theoretical structures. When metal cations classified as hard acids were combined with dCG x dCG, the [M(dCG x dCG)]+ complex organized into a globular structure. However, when soft acid metal cations were examined, a structure was observed where the two C-G base pairs were Watson-Crick bound.     

Cis Watson-Crick/Watson-Crick Base Pairs in DNA: A Density Functional Theory(DFT) Study

The four nucleic acid DNA bases(adenine, thymine, guanine, cytosine) and ten cis Watson-Crick/Watson-Crick(cis WC/WC) DNA base pairs were investigated by density functional theory(DFT) quantum chemical calculations. Geometry optimizations were carried out on the four bases and ten base pairs at the B3LYP level with 6-31G~(**) basis set. All the optimizations were performed within Cs symmetry. The optimum structures for the four bases and seven cis WC/WC base pairs were obtained, and Natural Bond Orbital analysis(NBO) was based on these structures. The possibilities of matches between any two of the four bases through their Watson-Crick(WC) edges were discussed. The structures of seven cis WC/WC base pairs change to a certain extent relative to these of the four bases due to the formation of hydrogen bonds. These base pairs existing in DNA have an important influence on the structural stability of the double helix. The analysis of the electronic structures and molecular orbitals for seven cis WC/WC base pairs can provide significant information about the relationship between charge transfer along the hydrogen bond and the Frontier orbitals of these base pairs.     

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.     

Four-stranded DNA structure stabilized by a novel G:C:A:T tetrad

The solution structure of a cyclic oligonucleotide d has been determined by two-dimensional NMR spectroscopy and restrained molecular dynamics. Under the appropriate experimental conditions, this molecule self-associates, forming a symmetric dimer stabilized by four intermolecular Watson-Crick base pairs. The resulting four-stranded structure consists of two G:C:A:T tetrads, formed by facing the minor groove side of the Watson-Crick base-pairs. Most probably, the association of the base-pairs is stabilized by coordinating a Na(+) cation. This is the first time that this novel G:C:A:T tetrad has been found in an oligonucleotide structure. This observation increases considerably the number of sequences that may adopt a four-stranded architecture. Overall, the three-dimensional structure is similar to those observed previously in other quadruplexes formed by minor groove alignment of Watson-Crick base pairs. This resemblance strongly suggests that we may be observing a general motif for DNA-DNA recognition.     

A TDDFT study of the excited states of DNA bases and their assemblies

We present a detailed study of the optical absorption spectra of DNA bases and base pairs, carried out by means of time dependent density functional theory. The spectra for the isolated bases are compared to available theoretical and experimental data and used to assess the accuracy of the method and the quality of the exchange-correlation functional. Our approach turns out to be a reliable tool to describe the response of the nucleobases. Furthermore, we analyze in detail the impact of hydrogen bonding and pi-stacking in the calculated spectra for both Watson-Crick base pairs and Watson-Crick stacked assemblies. We show that the reduction of the UV absorption intensity (hypochromicity) for light polarized along the base-pair plane depends strongly on the type of interaction. For light polarized perpendicular to the basal plane, the hypochromicity effect is reduced, but another characteristic is found, namely a blue shift of the optical spectrum of the base-assembly compared to that of the isolated bases. The use of optical tools as fingerprints for the characterization of the structure (and type of interaction) is extensively discussed.     

Janus-type AT nucleosides: synthesis, solid and solution state structures

Novel Janus-type nucleoside analogues (1a-d) were synthesized. Their pyrimido[4,5-d]pyrimidine base moiety has one face with a bidentate Watson-Crick donor-acceptor (DA) H-bond array of adenine and the other face with an acceptor-donor (AD) H-bond array of thymine. These nucleosides may self-associate through the self-complementary base pair. Indeed, in the solid state, compound 6d displayed a honeycomb-like supramolecular structure with tetrameric membered cavities formed through the combination of reverse Watson-Crick base pairs and aromatic stacking, in which the solvent molecules were accommodated. The result of temperature-dependent CD studies showed that the free nucleosides can form higher order chiral structures in aqueous solution.     

Thermodynamic stability of Hoogsteen and Watson-Crick base pairs in the presence of histone H3-mimicking peptide

We found that Hoogsteen base pairs were stabilized by molecular crowding and a histone H3-mimicking peptide, which was not observed for Watson-Crick base pairs. Our findings demonstrate that the type of DNA base pair is critical for the interaction between DNA and histones.     

Ultrafast repair of irradiated DNA: nonadiabatic ab initio simulations of the guanine-cytosine photocycle

Nonadiabatic first-principles molecular dynamics simulations have been performed of the photoexcited Watson-Crick guanine-cytosine (GC) DNA base pair in the gas phase and in aqueous solution. An excited state coupled proton-electron transfer (CPET) from G to C along the central hydrogen bond is observed upon excitation of the pipi* state initially localized on G. In the resulting charge transfer state a conical intersection between the excited state and the ground state is easily accessible. Therefore radiationless decay is fast, of the order of 100 fs, followed by a rapid CPET back reaction retrieving the initial Watson-Crick structure. A detailed analysis of the mechanism of nonradiative decay suggests a biexponential behavior in which out-of-plane motion plays a special role for the longer decay component.     

Synthesis and DNA duplex recognition of a triplex-forming oligonucleotide with an ureido-substituted 4-phenylimidazole nucleoside

A 4-(3-n-butylureidophenyl)imidazole nucleoside was successfully incorporated into a triplex-forming oligonucleotide (TFO). Binding affinity and base pair selectivity of the TFO containing this non-natural nucleoside were studied with various duplex targets containing all four possible Watson-Crick base pairs opposite the nucleoside analog in the third strand. Triplex thermal stabilities indicate that the synthetic nucleoside acts as a universal base in binding to all four possible Watson-Crick base pairs with moderate affinity but poor selectivity. Based on an analysis of its binding thermodynamics, this can be rationalized by the absence of strong specific interactions and more favorable entropic contributions upon triplex formation.     

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.     

Crystal structure of a partly self-complementary peptide nucleic acid (PNA) oligomer showing a duplex-triplex network

The X-ray structure of a partly self-complementary peptide nucleic acid (PNA) decamer (H-GTAGATCACT-l-Lys-NH(2)) to 2.60 A resolution is reported. The structure is mainly controlled by the canonical Watson-Crick base pairs formed by the self-complementary stretch of four bases in the middle of the decamer (G(4)A(5)T(6)C(7)). One right- and one left-handed Watson-Crick duplex are formed. The two PNA units C(9)T(10) change helical handedness, so that each PNA strand contains both a right- and a left-handed section. The changed handedness in C(9)T(10) allows formation of Hoogsteen hydrogen bonding between C(9)T(10) and G(4)A(5) of a PNA strand in an adjacent Watson-Crick double helix of the same handedness. Thereby, a PNA-PNA-PNA triplex is formed. The PNA unit A(3) forms a noncanonical base pair with A(8) in a symmetry-related strand of opposite handedness; the base pair is of the A-A reverse Hoogsteen type. The structural diversity of this PNA demonstrates how the PNA backbone is able to adapt to structures governed by the stacking and hydrogen-bonding interactions between the nucleobases. The crystal structure further shows how PNA oligomers containing limited sequence complementarity may form complex hydrogen-bonding networks.     

Duplex structure of a minimal nucleic acid

The crystal structure of an 8-mer (S)-GNA duplex is presented. As a tool for phasing, the anomalous diffraction of two copper(II) ions within two artificial metallo-base pairs was employed. The duplex structure confirms a canonical Watson-Crick base pairing scheme of GNA with antiparallel strands. The duplex secondary structure is distinct from canonical A- and B-form nucleic acids and can be described as a right-handed helical ribbon wrapped around the helix axis, resulting in a large hollow core. Most intriguingly, neighboring base pairs slide strongly against each other, resulting in extensive interstrand base-base hydrophobic interactions along with unusual hydrophobic intrastrand interactions of nucleobases with their backbone. These results reveal how a minimal nucleic acid backbone can support highly stable Watson-Crick-like duplex formation.     

Hydrogen-bonding interactions in peptide nucleic acid and deoxyribonucleic acid: a comparative study

Peptide nucleic acid (PNA) is a synthetic analogue of deoxyribonucleic acid (DNA) capable of tightly binding to itself and DNA with high specificity. Using hybrid density functional methods, hydrogen-bond (H-bond) strengths have been evaluated for isolated Watson-Crick base pairs, PNA base pairs, and charged as well as neutral DNA base pairs. Heterogeneous base pairs of PNA with charged and neutral DNA have also been investigated. The competing effects of short-range H-bonding and long-range Coulombic repulsions in charged DNA base pairs have been analyzed. Polarizable continuum models have been employed to evaluate solvation effects on the binding energies.     

Photophysical characters of rationally designed hetero-ring-expanded guanine analogues and effect of cytosine pairing

We present the results of the CIS and TDB3LYP calculations of the optical absorption and emission spectra of some newly designed guanine (G) analogues and their Watson-Crick base pairs. Compared with natural G, the onset absorption peaks of these newly designed analogues are red-shifted, while the fluorescence peaks are blue-shifted. In general, the first excited singlet states (pipi*) of these analogues are nonplanar for all bases considered here. But, the Stokes shifts for the designed G-analogues are much smaller than that of natural G, suggesting that they have stronger molecular rigidity and higher fluorescence quantum yields than those of natural G. The first excited states of these Watson-Crick base pairs essentially originate from those of their isolated purine moieties, as demonstrated from the S1 geometries of their Watson-Crick base pairs. For G and its analogues, A1 and A2 (they are ring-expanded with one-bond intercalation at the C5 site), the pairing with cytosine reduces the oscillator strengths of both the first absorption peak (by 27%-60%) and the fluorescent emission (by 19%-23%), while for the analogues A3, A4, and xG in which G is ring-expanded with a two-bond intercalation at the C5 site, the pairing, in contrast, increases the oscillator strengths of both the first absorption peak (by 11%-15%) and the fluorescent emission (by 3%-20%). These observations indicate that the pairing with cytosine can quench the fluorescence for G, A1, and A2 but enhance the fluorescence quantum yields for A3, A4, and xG. The significant shifts induced by ring-expansion in the ring-expanded G with a two-bond intercalation at the C5 site reveal a possibility for their fluorescent detections.     

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.