Stacking effects on local structure in RNA: changes in the structure of tandem GA pairs when flanking GC pairs are replaced by isoG-isoC pairs

The Watson-Crick-like isoG-isoC (iGiC) pair, with the amino and carbonyl groups transposed relative to the Watson-Crick GC pair, provides an expanded alphabet for understanding interactions that shape nucleic acid structure. Here, thermodynamic stabilities of tandem GA pairs flanked by iGiC pairs are reported along with the NMR structures of the RNA self-complementary duplexes (GCiGGAiCGCA)2 and (GGiCGAiGCCA)2. A sheared GA pairing forms in (GCiGGAiCGCA)2, and an imino GA pairing forms in (GGiCGAiGCCA)2. The structures contrast with the formation of tandem imino and sheared GA pairs flanked by GC pairs in the RNA self-complementary duplexes (GCGGACGC)2 and (GGCGAGCC)2, respectively. In both iGiC duplexes, Watson-Crick-like hydrogen bonds are formed between iG and iC, and iGiC substitutions result in less favorable loop stability. The results provide benchmarks for testing computations of molecular interactions that shape RNA three-dimensional structure.     

Cyclohexyl "base pairs" stabilize duplexes and intensify pyrene fluorescence by shielding it from natural base pairs

In this study, we investigated the stability and structure of artificial base pairs that contain cyclohexyl rings. The introduction of a single pair of isopropylcyclohexanes into the middle of DNA slightly destabilized the duplex. Interestingly, as the number of the "base pairs" increased, the duplex was remarkably stabilized. A duplex with six base pairs was even more stable than one containing six A-T pairs. Thermodynamic analysis revealed that changes in entropy and not enthalpy contributed to duplex stability, demonstrating that hydrophobic interactions between isopropyl groups facilitated the base pairing, and thus stabilized the duplex. NOESY of a duplex containing an isopropylcyclohexane-methylcyclohexane pair unambiguously demonstrated its "pairing" in the duplex because distinct NOEs between the protons of cyclohexyl moieties and imino protons of both of the neighboring natural base pairs were observed. CD spectra of duplexes tethering cyclohexyl moieties also showed a positive-negative couplet that is characteristic of the B-form DNA duplex. Taken together, these results showed that cyclohexyl moieties formed base pairs in the DNA duplex without severely disturbing the helical structure of natural DNA. Next, we introduced cyclohexyl base pairs between pyrene and nucleobases as an "insulator" that suppresses electron transfer between them. We found a massive increase in the quantum yield of pyrene due to the efficient shielding of pyrene from nucleobases. The cyclohexyl base pairs reported here have the potential to prepare highly fluorescent labeling agents by multiplying fluorophores and insulators alternately into DNA duplexes.     

Revision of AMBER Torsional Parameters for RNA Improves Free Energy Predictions for Tetramer Duplexes with GC and iGiC Base Pairs

All-atom force fields are important for predicting thermodynamic, structural, and dynamic properties of RNA. In this paper, results are reported for thermodynamic integration calculations of free energy differences of duplex formation when CG pairs in the RNA duplexes r(CCGG)(2), r(GGCC)(2), r(GCGC)(2), and r(CGCG)(2) are replaced by isocytidine-isoguanosine (iCiG) pairs. Agreement with experiment was improved when ε/ζ, α/γ, β, and χ torsional parameters in the AMBER99 force field were revised on the basis of quantum mechanical calculations. The revised force field, AMBER99TOR, brings free energy difference predictions to within 1.3, 1.4, 2.3, and 2.6 kcal/mol at 300 K, respectively, compared to experimental results for the thermodynamic cycles of CCGG → iCiCiGiG, GGCC → iGiGiCiC, GCGC → iGiCiGiC, and CGCG → iCiGiCiG. In contrast, unmodified AMBER99 predictions for GGCC → iGiGiCiC and GCGC → iGiCiGiC differ from experiment by 11.7 and 12.6 kcal/mol, respectively. In order to test the dynamic stability of the above duplexes with AMBER99TOR, four individual 50 ns molecular dynamics (MD) simulations in explicit solvent were run. All except r(CCGG)(2) retained A-form conformation for ≥82% of the time. This is consistent with NMR spectra of r(iGiGiCiC)(2), which reveal an A-form conformation. In MD simulations, r(CCGG)(2) retained A-form conformation 52% of the time, suggesting that its terminal base pairs may fray. The results indicate that revised backbone parameters improve predictions of RNA properties and that comparisons to measured sequence dependent thermodynamics provide useful benchmarks for testing force fields and computational methods.     

Dynamic behavior of DNA base pairs containing 8-oxoguanine

The process by which DNA repair enzymes recognize and selectively excise damaged bases in duplex DNA is fundamental to our mechanistic understanding of these critical biological reactions. 8-Oxoguanine (8-oxoG) is the most common form of oxidative DNA damage; unrepaired, this lesion generates a G:C-->T:A mutation. Central to the recognition and repair of DNA damage is base extrusion, a process in which the damaged base lesion or, in some cases, its partner disengages from the helix and is bound to the enzyme's active site where base excision takes place. The conformation adopted by 8-oxoG in duplex DNA is affected by the base positioned opposite this lesion; conformational changes may also take place when the damaged base binds to its cognate repair enzyme. We performed unrestrained molecular dynamics simulations for several 13-mer DNA duplexes. Oligomers containing G:C and 8oxoG:C pairs adopted Watson-Crick geometries in stable B-form duplexes; 8oxoG showed increased local and global flexibility and a reduced barrier to base extrusion. Duplexes containing the G:A mismatch showed much larger structural fluctuations and failed to adopt a well-defined structure. For the 8oxoG:A mismatch that is recognized by the DNA glycosylase MutY, the damaged nucleoside underwent spontaneous and reproducible anti-->syn transitions. The syn conformation is thermodynamically preferred. Steric hindrance and unfavorable electrostatics associated with the 8oxoG O8 atom in the anti conformation were the major driving forces for this transition. Transition events follow two qualitatively different pathways. The overall anti-->syn transition rate and relative probability of the two transition paths were dependent on local sequence context. These simulations indicate that both the dynamic and equilibrium behavior of the duplex change as a result of oxidation; these differences may provide valuable new insight into the selective action of enzymes on damaged DNA.     

Theoretical Studies on the Intermolecular Interactions of Potentially Primordial Base‐Pair Analogues

Recent experimental studies on the Watson–Crick type base pairing of triazine and aminopyrimidine derivatives suggest that acid/base properties of the constituent bases might be related to the duplex stabilities measured in solution. Herein we use high‐level quantum chemical calculations and molecular dynamics simulations to evaluate the base pairing and stacking interactions of seven selected base pairs, which are common in that they are stabilized by two N? H???O hydrogen bonds separated by one N? H???N hydrogen bond. We show that neither the base pairing nor the base stacking interaction energies correlate with the reported pK_a data of the bases and the melting points of the duplexes. This suggests that the experimentally observed correlation between the melting point data of the duplexes and the pK_a values of the constituent bases is not rooted in the intrinsic base pairing and stacking properties. The physical chemistry origin of the observed experimental correlation thus remains unexplained and requires further investigations. In addition, since our calculations are carried out with extrapolation to the complete basis set of atomic orbitals and with inclusion of higher electron correlation effects, they provide reference data for stacking and base pairing energies of non‐natural bases.     

Sequence‐Dependent Duplex Stabilization upon Formation of a Metal‐Mediated Base Pair

An artificial nucleoside surrogate with 1H‐imidazo[4,5‐f][1,10]phenanthroline ( P ) acting as an aglycone has been introduced into DNA oligonucleotide duplexes. This nucleoside surrogate can act as a bidentate ligand, and so is useful in the context of metal‐mediated base pairs. Several duplexes involving a hetero base pair with an imidazole nucleoside have been investigated. The stability of DNA duplexes incorporating the respective Ag^I‐mediated base pairs strongly depends on the sequence context. Quantum mechanical/molecular mechanical (QM/MM) calculations have been performed in order to gain insight into the factors determining this sequence dependence. The results indicated that, in addition to the stabilizing effect that results from the formation of coordinative bonds, destabilizing effects may occur when the artificial base pair does not fit optimally into the surrounding B‐DNA duplex.     

Marked dependence on carrier-ligand bulk but not on carrier-ligand chirality of the duplex versus single-strand forms of a DNA oligonucleotide with a series of G-Pt(II)-G intrastrand cross-links modeling cisplatin-DNA adducts

The N7-Pt-N7 adjacent G,G intrastrand DNA cross-link responsible for cisplatin anticancer activity is dynamic, promotes local "melting" in long DNA, and converts many oligomer duplexes to single strands. For 5'-d(A1T2G3G4G5T6A7C8C9C10A11T12)-3' (G3), treatment of the (G3)2 duplex with five pairs of [LPt(H2O)2]2+ enantiomers (L = an asymmetric diamine) formed mixtures of LPt-G3 products (1 Pt per strand) cross-linked at G3,G4 or at G4,G5 in all cases. L chirality exerted little influence. For primary diamines L with bulk on chelate ring carbons (e.g., 1,2-diaminocyclohexane), the duplex was converted completely into single strands (G3,G4 coils and G4,G5 hairpins), exactly mirroring results for cisplatin, which lacks bulk. In sharp contrast, for secondary diamines L with bulk on chelate ring nitrogens (e.g., 2,2'-bipiperidine, Bip), unexpectedly stable duplexes having two platinated strands (even a unique G3,G4/G4,G5 heteroduplex) were formed. After enzymatic digestion of BipPt-G3 duplexes, the conformation of the relatively nondynamic G,G units was shown to be head-to-head (HH) by HPLC/mass spectrometric characterization. Because the HH conformation dominates at the G,G lesion in duplex DNA and in the BipPt-G3 duplexes, the stabilization of the duplex form only when the L nitrogen adducts possess bulk suggests that H-bonding interactions of the Pt-NH groups with the flanking DNA lead to local melting and to destabilization of oligomer duplexes. The marked dependence of adduct properties on L bulk and the minimal dependence on L chirality underscore the need for future exploration of the roles of the L periphery in affecting anticancer activity.     

New base pairing motifs. The synthesis and thermal stability of oligodeoxynucleotides containing imidazopyridopyrimidine nucleosides with the ability to form four hydrogen bonds

The synthesis and thermal stability of oligodeoxynucleotides (ODNs) containing imidazo[5',4':4,5]pyrido[2,3-d]pyrimidine nucleosides 1-4 (N(N), O(O), N(O), and O(N), respectively) with the aim of developing two sets of new base pairing motifs consisting of four hydrogen bonds (H-bonds) is described. The proposed four tricyclic nucleosides 1-4 were synthesized through the Stille coupling reaction of a 5-iodoimidazole nucleoside with an appropriate 5-stannylpyrimidine derivative, followed by an intramolecular cyclization. These nucleosides were incorporated into ODNs to investigate the H-bonding ability. When one molecule of the tricyclic nucleosides was incorporated into the center of each ODN (ODN I and II, each 17mer), no apparent specificity of base pairing was observed, and all duplexes were less stable than the duplexes containing natural G:C and A:T pairs. On the other hand, when three molecules of the tricyclic nucleosides were consecutively incorporated into the center of each ODN (ODN III and IV, each 17mer), thermal and thermodynamic stabilization of the duplexes due to the specific base pairings was observed. The melting temperature (T(m)) of the duplex containing the N(O):O(N) pairs showed the highest T(m) of 84.0 degrees C, which was 18.2 and 23.5 degrees C higher than that of the duplexes containing G:C and A:T pairs, respectively. This result implies that N(O)and O(N) form base pairs with four H-bonds when they are incorporated into ODNs. The duplex containing N(O):O(N) pairs was markedly stabilized by the assistance of the stacking ability of the imidazopyridopyrimidine bases. Thus, we developed a thermally stable new base pairing motif, which should be useful for the stabilization and regulation of a variety of DNA structures.     

Stability and structure of RNA duplexes containing isoguanosine and isocytidine.

Isoguanosine (iG) and isocytidine (iC) differ from guanosine (G) and cytidine (C), respectively, in that the amino and carbonyl groups are transposed. The thermodynamic properties of a set of iG, iC containing RNA duplexes have been measured by UV optical melting. It is found that iG-iC replacements usually stabilize duplexes, and the stabilization per iG-iC pair is sequence-dependent. The sequence dependence can be fit to a nearest-neighbor model in which the stabilities of iG--iC pairs depend on the adjacent iG--iC or G--C pairs. For 5'-CG-3'/3'-GC-5' and 5'-GG-3'/3'-CC-5' nearest neighbors, the free energy differences upon iG-iC replacement are smaller than 0.2 kcal/mol at 37 degrees C, regardless of the number of replacements. For 5'-GC-3'/3'-CG-5', however, each iG--iC replacement adds 0.6 kcal/mol stabilizing free energy at 37 degrees C. Stacking propensities of iG and iC as unpaired nucleotides at the end of a duplex are similar to those of G and C. An NMR structure is reported for r(CiGCGiCG)(2) and found to belong to the A-form family. The structure has substantial deviations from standard A-form but is similar to published NMR and/or crystal structures for r(CGCGCG)(2) and 2'-O-methyl (CGCGCG)(2). These results provide benchmarks for theoretical calculations aimed at understanding the fundamental physical basis for the thermodynamic stabilities of nucleic acid duplexes.     

The 5-Me of thyminyl (T) interaction with the neighboring nucleobases dictate the relative stability of isosequential DNA-RNA hybrid duplexes

Relative energetic contributions from the base-pairing [DeltaG(degrees)(bp)] vis-a-vis stacking [DeltaG(degrees)(stacking)] to the total free-energy of stabilization [DeltaG(degrees)(37)] for 14 pairs of isosequential hybrid DNA-RNA duplexes (taken from E. A. Lesnik and S. M. Freier, Biochemistry ,1995, 34, 10807) have been dissected in order to understand the differences in the intrinsic nature of the electrostatic forces that are responsible for the self-assembly of the heteroduplexes compared to homoduplexes. The pKa differences between the monomeric nucleotide 3'-ethylphosphates [(d/rN)pEt] as well as nucleotide 3',5'-bis-ethylphosphates [Etp(d/rN)pEt] in both 2'-deoxy (dN) and ribo (rN) series (N = A/G/C/T/U), as the model donor and acceptor (in which stacking is completely eliminated) mimicking those of the internucleotide monomer building blocks of a duplex, can be qualitatively used (P. Acharya, P. Cheruku, S. Chatterjee, S. Acharya and J. Chattopadhyaya, J. Am. Chem. Soc., 2004, 126, 2862) to understand the strength of base-pairing energies in different DNA-RNA (DR), RNA-DNA (RD), DNA-DNA (DD), and RNA-RNA (RR) duplexes. The study has led us to show the following. (1) As the number of excess %T in DR duplexes compared to the isosequential RD duplexes increase the differences in their thermal stabilization [DeltaT(m)](DR-RD) increase and vice-versa(2) The total relative stabilizations, [DeltaDeltaG(degrees)(37)](DR-RD) among the 14 pairs of isosequential DR and RD duplexes (E. A. Lesnik and S. M. Freier, Biochemistry ,1995, 34, 10807) are wholly dependent on the differences in the number of 5-Me(T) stacking interactions with the nearest-neighbors in the D strands of DR duplexes compared to that of the RD duplexes (3) In the relative stabilization of the DR or RD duplexes differences in the free-energy of stackings [DeltaDeltaG(degrees)(stacking)](DR-RD) play a more significant role than the differences in the free-energy of base-pairing, [DeltaDeltaG(degrees)(bp)](DR-RD). In contradistinction, our experimental data shows that RNA-RNA duplexes are more stable than DNA-DNA duplexes because of larger energy gain from the base-pairing in the former compared to the latter (P. Acharya, P. Cheruku, S. Chatterjee, S. Acharya and J. Chattopadhyaya, J. Am. Chem. Soc., 2004, 126, 2862).     

B-DNA Helix Stability in a Solvent-Free Environment

B-DNA is the most common DNA helix conformation under physiological conditions. However, when the amount of water in a DNA solution is decreased, B-to-A helix transitions have been observed. To understand what type of helix conformations exist in a solvent-free environment, a series of poly d(CG)(n) and mixed sequence DNA duplexes from 18 to 30 bp were examined with circular dichroism (CD), ESI-MS, ion mobility, and molecular dynamics. From the CD spectra, it was observed that all sequences had B-form helices in solution. However, the solvent-free results were more complex. For the poly d(CG)(n) series, the 18 bp duplex had an A-form helix conformation, both A- and B-helices were present for the 22 bp duplex, and only B-helices were observed for the 26 and 30 bp duplexes. Since these sequences were all present as B-DNA in solution, the observed solvent-free structures illustrate that smaller helices with fewer base pairs convert to A-DNA more easily than larger helices in the absence of solvent. A similar trend was observed for the mixed sequence duplexes where both an A- and B-helix were present for the 18 bp duplex, while only B-helices occur for the larger 22, 26, and 30 bp duplexes. Since the solvent-free B-helices appear at smaller sizes for the mixed sequences than for the pure d(CG)(n) duplexes, the pure d(CG)(n) duplexes have a greater A-philicity.     

NMR study on the photoresponsive DNA tethering an azobenzene. Assignment of the absolute configuration of two diastereomers and structure determination of their duplexes in the trans-form

Two diastereomers of a photoresponsive oligodeoxyribonucleotide tethering a trans-azobenzene, based on the chirality of the central carbon of a diol linker, were separated by reversed-phase HPLC. On the basis of 2D NMR analysis, absolute configurations of the diastereomers alpha and beta (tentatively designated from differences in their retention time) were determined as R- and S-forms, respectively. For both diastereomers, their NMR-determined duplex structure showed that trans-azobenzene intercalates between base pairs, because distinct NOEs were observed between the protons of azobenzene and those of the adjacent base pairs, such as with the imino protons and methyl protons of thymine. The melting temperatures of both duplexes were higher than that of the corresponding native duplex, which contained no azobenzene residue, due to the intercalated trans-azobenzene stabilizing the duplex by a stacking interaction. Between these two diastereomers, differences in T(m) were also found: the melting temperature of the R-form duplex (alpha-isomer) was higher than that of the S-form (beta-isomer). On the basis of the NMR-determined structure, this difference was attributed to the fact that the S-form (beta isomer) causes more stress forming the duplex than does the R-form (alpha isomer) due to disturbances of the right-hand helix.     

Bipyridyl- and biphenyl-DNA: a recognition motif based on interstrand aromatic stacking

The synthesis and incorporation into oligonucleotides of C-nucleosides containing the two aromatic, non-hydrogen-bonding nucleobase substitutes biphenyl (I) and bipyridyl (Y) are described. Their homo- and hetero-recognition properties in different sequential arrangements were then investigated via UV-melting curve analysis, gel mobility assays, CD- and NMR spectroscopy. An NMR analysis of a dodecamer duplex containing one biphenyl pair in the center, as well as CD data on duplexes with multiple insertions provide further evidence for the zipper-like interstrand stacking motif that we proposed earlier based on molecular modeling. UV-thermal melting experiments with duplexes containing one to up to seven I- or Y base pairs revealed a constant increase in T(m) in the case of I and a constant decrease for Y. Mixed I/Y base pairs lead to stabilities in between the homoseries. Insertion of alternating I/abasic site- or Y/abasic site pairs strongly decreases the thermal stability of duplexes. Asymmetric distribution of I- or Y residues on either strand of the duplex were also investigated in this context. Duplexes with three natural base pairs at both ends and 50 % of I pairs in the center are still readily formed, while duplexes with blunt ended I pairs tend to aggregate unspecifically. Duplexes with one natural overhang at the end of a I-I base pair tract can both aggregate or form ordered duplexes, depending on the nature of the natural bases in the overhang.     

Molecular modeling on recognition of sheared and normal DNA by novel metal complex Λ- and Δ-[Co(phen)2hpip]

Recognition of sheared and normal DNA by a novel metal complex [Co(phen)₂hpip]³⁺ (phen=1,10-phenanthroline, HPIP=2-(2-hydroxyphenyl)imidazole[4,5-f][1,10]phenanethroline) is studied by molecular modeling. Calculating results indicate that, this complex can specifically recognize DNA segment of sequence –MMNNMM– (M means mismatch base pairs and N means normal base pairs). Intercalating from minor groove between the middle normal duplex into the sheared DNA with the depth of 1.2 nm is of preference and enantioselectivity is observed. Comparison on the two DNA structures of optimal conformation and analysis on the interaction between DNA and the two tail ligands of the complex show that, the effect of the two neighboring mismatch duplexes on the structure of the middle normal base pairs and the steric interaction between the mismatch duplexes and the two tail ligands of the complex are the essential reason to the segment specificity. Investigation on the detailed energy terms indicate that, in effecting enantioselectivity, the electrostatic distribution of the complex is in the majority and steric interaction is at the next place. But, steric interaction is surely the only factor determining the intercalating from minor groove.     

Dramatic 5'-residue effect on conformer distribution of short oligonucleotide retro models of the cisplatin-DNA cross-link: implications for the Lippard and cross-link distorted base pair steps present in cisplatin-DNA duplex adducts

The cisplatin anticancer drug preferentially attacks the GG sequence of DNA duplexes. Virtually all DNAs containing the key G*G* lesion (G* = N7 platinated G) have large distortions in the cross-link (G*G*) base pair (bp) step and also in the adjacent Lippard (XG*) bp step, making the adducts very different from B-form DNA in the XG*G* region. The XG*G* strand in duplexes also differs in several ways from single-strand (ss) models with G*G* and XG*G* sequences. In the duplex, the X residue has an N sugar, the 5'-G* and 3'-G* bases have slight "R" canting (3'-G* H8 atom toward the 5'-G* base), and there is no or weak H-bonding by the NH3 ligands. In most XG*G* ss models, X has an S sugar, the 5'-G* base normally cants strongly toward the 3'-G* base (L canting), and the NH3 forms an H-bond. Well-defined ss models exist in the solid state, but dynamic motion obscures the properties of the ss models in solution. In this work, we employ retro models (better defined, less dynamic ss models) to understand the differences between duplex and ss models. The retro models in this study lack carrier ligand NH's, thus eliminating H-bonding. To correlate previous ss solid-state models with our solution work, we constructed hybrid molecules by overlaying parts of known structures. The combined model and experimental information indicates that the X N-pucker is not favorable in L-canted ss models, that X residue steric effects (not H-bonding) favor L canting in ss models, that X N-pucker is needed for favorable WC hydrogen bonding and stacking interactions in duplexes, and that X N-pucker minimizes X base clashes with bases in the complementary strand in duplexes. The R canting minimizing clashes between the X and G* residues of the Lippard bp step (independent of X pucker) and the repositioning of the X residue base caused by the change from S-pucker to N-pucker together lead to the unusual features of the Lippard bp step in the duplex.     

Neomycin-neomycin dimer: an all-carbohydrate scaffold with high affinity for AT-rich DNA duplexes

A dimeric neomycin-neomycin conjugate 3 with a flexible linker, 2,2'-(ethylenedioxy)bis(ethylamine), has been synthesized and characterized. Dimer 3 can selectively bind to AT-rich DNA duplexes with high affinity. Biophysical studies have been performed between 3 and different nucleic acids with varying base composition and conformation by using ITC (isothermal calorimetry), CD (circular dichroism), FID (fluorescent intercalator displacement), and UV (ultraviolet) thermal denaturation experiments. A few conclusions can be drawn from this study: (1) FID assay with 3 and polynucleotides demonstrates the preference of 3 toward AT-rich sequences over GC-rich sequences. (2) FID assay and UV thermal denaturation experiments show that 3 has a higher affinity for the poly(dA)·poly(dT) DNA duplex than for the poly(dA)·2poly(dT) DNA triplex. Contrary to neomycin, 3 destabilizes poly(dA)·2poly(dT) triplex but stabilizes poly(dA)·poly(dT) duplex, suggesting the major groove as the binding site. (3) UV thermal denaturation studies and ITC experiments show that 3 stabilizes continuous AT-tract DNA better than DNA duplexes with alternating AT bases. (4) CD and FID titration studies show a DNA binding site size of 10-12 base pairs/drug, depending upon the structure/sequence of the duplex for AT-rich DNA duplexes. (5) FID and ITC titration between 3 and an intramolecular DNA duplex [d(5'-A(12)-x-T(12)-3'), x = hexaethylene glycol linker] results in a binding stoichiometry of 1:1 with a binding constant ～10(8) M(-1) at 100 mM KCl. (6) FID assay using 3 and 512 hairpin DNA sequences that vary in their AT base content and placement also show a higher binding selectivity of 3 toward continuous AT-rich than toward DNA duplexes with alternate AT base pairs. (7) Salt-dependent studies indicate the formation of three ion pairs during binding of the DNA duplex d[5'-A(12)-x-T(12)-3'] and 3. (8) ITC-derived binding constants between 3 and DNA duplexes have the following order: AT continuous, d[5'-G(3)A(5)T(5)C(3)-3'] > AT alternate, d[5'-G(3)(AT)(5)C(3)-3'] > GC-rich d[5'-A(3)G(5)C(5)T(3)-3']. (9) 3 binds to the AT-tract-containing DNA duplex (B* DNA, d[5'-G(3)A(5)T(5)C(3)-3']) with 1 order of magnitude higher affinity than to a DNA duplex with alternating AT base pairs (B DNA, d[5'-G(3)(AT)(5)C(3)-3']) and with almost 3 orders of magnitude higher affinity than a GC-rich DNA (A-form, d[5'-A(3)G(5)C(5)T(3)-3']).     

Recognition properties of donor- and acceptor-modified biphenyl-DNA

The recognition properties of DNA duplexes containing single or triple incorporations of eight different donor-modified (OMe, NH(2)) and acceptor-modified (NO(2)) biphenyl residues as base replacements in opposite positions were probed by UV-melting and by CD and fluorescence spectroscopy. We found a remarkable dependence of duplex stability on the natures of the substituents (donor vs. acceptor). The stabilities of duplexes with one biphenyl pair increase in the order donor/donor < acceptor/donor < acceptor/acceptor substitution. The most stable biphenyl pairs stabilize duplexes by up to 6 degrees C in T(m). In duplexes with three consecutive biphenyl pairs the stability increases in the inverse order (acceptor/acceptor < donor/acceptor < donor/donor) with increases in T(m), relative to an unmodified duplex, of up to 10 degrees C. A thermodynamic analysis, combined with theoretical calculations of the physical properties of the biphenyl substituents, suggests that in duplexes with single biphenyl pairs the affinity is dominated by electrostatic forces between the biphenyl/nearest neighbor natural base pairs, whereas in the triple-modified duplexes the increase in thermal stability is predominantly determined by hydrophobic interactions of the biphenyl residues with each other. Oligonucleotides containing amino biphenyl residues are fluorescent. Their fluorescence is largely quenched when they are paired with themselves or with nitrobiphenyl-containing duplex partners.     

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.     

Effect of nucleobase sequence on the proton-transfer reaction and stability of the guanine-cytosine base pair radical anion

The formation of base pair radical anions is closely related to many fascinating research fields in biology and chemistry such as radiation damage to DNA and electron transport in DNA. However, the relevant knowledge so far mainly comes from studies on isolated base pair radical anions, and their behavior in the DNA environment is less understood. In this study, we focus on how the nucleobase sequence affects the properties of the guanine-cytosine (GC) base pair radical anion. The energetic barrier and reaction energy for the proton transfer along the N(1)(G)-H···N(3)(C) hydrogen bond and the stability of GC˙(-) (i.e., electron affinity of GC) embedded in different sequences of base-pair trimer were evaluated using density functional theory. The computational results demonstrated that the presence of neighboring base pairs has an important influence on the behavior of GC˙(-) in the gas phase. The excess electron was found to be localized on the embedded GC and the charge leakage to neighboring base pairs was very minor in all of the investigated sequences. Accordingly, the sequence behavior of the proton-transfer reaction and the stability of GC˙(-) is chiefly governed by electrostatic interactions with adjacent base pairs. However, the effect of base stacking, due to its electrostatic nature, is severely screened upon hydration, and thus, the sequence dependence of the properties of GC˙(-) in aqueous environment becomes relatively weak and less than that observed in the gas phase. The effect of geometry relaxation associated with neighboring base pairs as well as the possibility of proton transfer along the N(2)(G)-H···O(2)(C) channel have also been investigated. The implications of the present findings to the electron transport and radiation damage of DNA are discussed.