Structural effect of the anticancer agent 6-thioguanine on duplex DNA

The incorporation of 6-thioguanine (S6G) into DNA is an essential step in the cytotoxic activity of thiopurines. However, the structural effects of this substitution on duplex DNA have not been fully characterized. Here, we present the solution structures of DNA duplexes containing S6G opposite thymine (S6G·T) and opposite cytosine (S6G·C), solved by high-resolution NMR spectroscopy and restrained molecular dynamics. The data indicate that both duplexes adopt right-handed helical conformations with all Watson–Crick hydrogen bonding in place. The S6G·T structures exhibit a wobble-type base pairing at the lesion site, with thymine shifted toward the major groove and S6G displaced toward the minor groove. Aside from the lesion site, the helices, including the flanking base pairs, are not highly perturbed by the presence of the lesion. Surprisingly, thermal dependence experiments suggest greater stability in the S6G-T mismatch than the S6G-C base pair.     

Ion-induced folding of a kink turn that departs from the conventional sequence

Kink turns (k-turns) are important structural motifs that create a sharp axial bend in RNA. Most conform to a consensus in which a three-nucleotide bulge is followed by consecutive G•A and A•G base pairs, and when these G•A pairs are modified in vitro this generally leads to a failure to adopt the k-turn conformation. Kt-23 in the 30S ribosomal subunit of Thermus thermophilus is a rare exception in which the bulge-distal A•G pair is replaced by a non-Watson–Crick A•U pair. In the context of the ribosome, Kt-23 adopts a completely conventional k-turn geometry. We show here that this sequence is induced to fold into a k-turn structure in an isolated RNA duplex by Mg²⁺ or Na⁺ ions. Therefore, the Kt-23 is intrinsically stable despite lacking the key A•G pair; its formation requires neither tertiary interactions nor protein binding. Moreover, the Kt-23 k-turn is stabilized by the same critical hydrogen-bonding interactions within the core of the structure that are found in more conventional sequences such as the near-consensus Kt-7. T. thermophilus Kt-23 has two further non-Watson–Crick base pairs within the non-canonical helix, three and four nucleotides from the bulge, and we find that the nature of these pairs influences the ability of the RNA to adopt k-turn conformation, although the base pair adjacent to the A•U pair is more important than the other.     

Free energy contributions of G.U and other terminal mismatches to helix stability

Thermodynamic parameters of helix formation were measured spectroscopically for seven hexaribonucleotides containing a GC tetramer core and G.U or other terminal mismatches. The free energies of helix formation are compared with those for the tetramer core alone and with those for the hexamer with six Watson-Crick base pairs. In 1 M NaCl, at 37 degrees C, the free energy of a terminal G.U mismatch is about equal to that of the corresponding A.U pair. Although other terminal mismatches studied add between -1.0 and -1.6 kcal/mol to delta G0 37 for helix formation, all are less stable than the corresponding Watson-Crick pairs. Comparisons of the stability increments for terminal G.U mismatches and G.C pairs suggest when stacking is weak the additional hydrogen bond in the G.C pair adds roughly -1 kcal/mol to the favorable free energy of duplex formation.     

Nearest-neighbor thermodynamics of deoxyinosine pairs in DNA duplexes

Nearest-neighbor thermodynamic parameters of the ‘universal pairing base’ deoxyinosine were determined for the pairs I·C, I·A, I·T, I·G and I·I adjacent to G·C and A·T pairs. Ultraviolet absorbance melting curves were measured and non-linear regression performed on 84 oligonucleotide duplexes with 9 or 12 bp lengths. These data were combined with data for 13 inosine containing duplexes from the literature. Multiple linear regression was used to solve for the 32 nearest-neighbor unknowns. The parameters predict the T_m for all sequences within 1.2°C on average. The general trend in decreasing stability is I·C > I·A > I·T ≈ I· G > I·I. The stability trend for the base pair 5′ of the I·X pair is G·C > C·G > A·T > T·A. The stability trend for the base pair 3′ of I·X is the same. These trends indicate a complex interplay between H-bonding, nearest-neighbor stacking, and mismatch geometry. A survey of 14 tandem inosine pairs and 8 tandem self-complementary inosine pairs is also provided. These results may be used in the design of degenerate PCR primers and for degenerate microarray probes.     

Thermodynamics of RNA duplexes with tandem mismatches containing a uracil-uracil pair flanked by C.G/G.C or G.C/A.U closing base pairs

The thermodynamics governing the denaturation of RNA duplexes containing 8 bp and a central tandem mismatch or 10 bp were evaluated using UV absorbance melting curves. Each of the eight tandem mismatches that were examined had one U-U pair adjacent to another noncanonical base pair. They were examined in two different RNA duplex environments, one with the tandem mismatch closed by G.C base pairs and the other with G.C and A.U closing base pairs. The free energy increments (Delta Gdegrees(loop)) of the 2 x 2 loops were positive, and showed relatively small differences between the two closing base pair environments. Assuming temperature-independent enthalpy changes for the transitions, (Delta Gdegrees(loop)) for the 2 x 2 loops varied from 0.9 to 1.9 kcal/mol in 1 M Na(+) at 37 degrees C. Most values were within 0.8 kcal/mol of previously estimated values; however, a few sequences differed by 1.2-2.0 kcal/mol. Single strands employed to form the RNA duplexes exhibited small noncooperative absorbance increases with temperature or transitions indicative of partial self-complementary duplexes. One strand formed a partial self-complementary duplex that was more stable than the tandem mismatch duplexes it formed. Transitions of the RNA duplexes were analyzed using equations that included the coupled equilibrium of self-complementary duplex and non-self-complementary duplex denaturation. The average heat capacity change (DeltaC(p)) associated with the transitions of two RNA duplexes was estimated by plotting DeltaH degrees and DeltaS degrees evaluated at different strand concentrations as a function of T(m) and ln T(m), respectively. The average DeltaC(p) was 70 +/- 5 cal K(-)(1) (mol of base pairs)(-)(1). Consideration of this heat capacity change reduced the free energy of formation at 37 degrees C of the 10 bp control RNA duplexes by 0.3-0.6 kcal/mol, which may increase Delta Gdegrees(loop) values by similar amounts.     

Exploration of the conserved A+C wobble pair within the ribosomal peptidyl transferase center using affinity purified mutant ribosomes

Protein synthesis in the ribosome's large subunit occurs within an active site comprised exclusively of RNA. Mutational studies of rRNA active site residues could provide valuable insight into the mechanism of peptide bond formation, but many of these mutations cause a dominant lethal phenotype, which prevents production of the homogeneous mutant ribosomes needed for analysis. We report a general method to affinity purify in vivo assembled 50S ribosomal subunits containing lethal active site mutations via a U1A protein-binding tag inserted onto the 23S rRNA. The expected pH-dependent formation of the A2450+C2063 wobble pair has made it a potential candidate for the pH-dependent conformational change that occurs within the ribosomal active site. Using this approach, the active site A2450+C2063 pair was mutated to the isosteric, but pH-independent, G2450•U2063 wobble pair, and 50S subunits containing the mutations were affinity purified. The G•U mutation caused the adjacent A2451 to become hyper-reactive to dimethylsulfate (DMS) modification in a pH-independent manner. Furthermore, the G•U mutation decreased both the rate of peptide bond formation and the affinity of the post-translocation complex for puromycin. The reaction rate (k_pep) was reduced ~200-fold for both puromycin and the natural aminoacyl-tRNA A-site substrate. The mutations also substantially altered the pH dependence of the reaction. Mutation of this base pair has significant deleterious effects upon peptidyl transferase activity, but because G•U mutation disrupts several tertiary contacts with the wobble pair, the assignment of A2450 as the active site residue with the neutral pK_a important for the peptidyl transferase reaction cannot be fully supported or excluded based upon these data.     

Base-stacking and base-pairing contributions into thermal stability of the DNA double helix

Two factors are mainly responsible for the stability of the DNA double helix: base pairing between complementary strands and stacking between adjacent bases. By studying DNA molecules with solitary nicks and gaps we measure temperature and salt dependence of the stacking free energy of the DNA double helix. For the first time, DNA stacking parameters are obtained directly (without extrapolation) for temperatures from below room temperature to close to melting temperature. We also obtain DNA stacking parameters for different salt concentrations ranging from 15 to 100 mM Na⁺. From stacking parameters of individual contacts, we calculate base-stacking contribution to the stability of A•T- and G•C-containing DNA polymers. We find that temperature and salt dependences of the stacking term fully determine the temperature and the salt dependence of DNA stability parameters. For all temperatures and salt concentrations employed in present study, base-stacking is the main stabilizing factor in the DNA double helix. A•T pairing is always destabilizing and G•C pairing contributes almost no stabilization. Base-stacking interaction dominates not only in the duplex overall stability but also significantly contributes into the dependence of the duplex stability on its sequence.     

Stability and mismatch discrimination of locked nucleic acid-DNA duplexes

Locked nucleic acids (LNA; symbols of bases, +A, +C, +G, and +T) are introduced into chemically synthesized oligonucleotides to increase duplex stability and specificity. To understand these effects, we have determined thermodynamic parameters of consecutive LNA nucleotides. We present guidelines for the design of LNA oligonucleotides and introduce free online software that predicts the stability of any LNA duplex oligomer. Thermodynamic analysis shows that the single strand-duplex transition is characterized by a favorable enthalpic change and by an unfavorable loss of entropy. A single LNA modification confines the local conformation of nucleotides, causing a smaller, less unfavorable entropic loss when the single strand is restricted to the rigid duplex structure. Additional LNAs adjacent to the initial modification appear to enhance stacking and H-bonding interactions because they increase the enthalpic contributions to duplex stabilization. New nearest-neighbor parameters correctly forecast the positive and negative effects of LNAs on mismatch discrimination. Specificity is enhanced in a majority of sequences and is dependent on mismatch type and adjacent base pairs; the largest discriminatory boost occurs for the central +C·C mismatch within the +T+C+C sequence and the +A·G mismatch within the +T+A+G sequence. LNAs do not affect specificity in some sequences and even impair it for many +G·T and +C·A mismatches. The level of mismatch discrimination decreases the most for the central +G·T mismatch within the +G+G+C sequence and the +C·A mismatch within the +G+C+G sequence. We hypothesize that these discrimination changes are not unique features of LNAs but originate from the shift of the duplex conformation from B-form to A-form.     

Nearest-neighbor thermodynamics and NMR of DNA sequences with internal A.A, C.C, G.G, and T.T mismatches

Thermodynamic measurements are reported for 51 DNA duplexes with A.A, C.C, G.G, and T.T single mismatches in all possible Watson-Crick contexts. These measurements were used to test the applicability of the nearest-neighbor model and to calculate the 16 unique nearest-neighbor parameters for the 4 single like with like base mismatches next to a Watson-Crick pair. The observed trend in stabilities of mismatches at 37 degrees C is G.G > T.T approximately A.A > C.C. The observed stability trend for the closing Watson-Crick pair on the 5' side of the mismatch is G.C >/= C.G >/= A.T >/= T.A. The mismatch contribution to duplex stability ranges from -2.22 kcal/mol for GGC.GGC to +2.66 kcal/mol for ACT.ACT. The mismatch nearest-neighbor parameters predict the measured thermodynamics with average deviations of DeltaG degrees 37 = 3.3%, DeltaH degrees = 7. 4%, DeltaS degrees = 8.1%, and TM = 1.1 degrees C. The imino proton region of 1-D NMR spectra shows that G.G and T.T mismatches form hydrogen-bonded structures that vary depending on the Watson-Crick context. The data reported here combined with our previous work provide for the first time a complete set of thermodynamic parameters for molecular recognition of DNA by DNA with or without single internal mismatches. The results are useful for primer design and understanding the mechanism of triplet repeat diseases.     

i-motif solution structure and dynamics of the d(AACCCC) and d(CCCCAA) tetrahymena telomeric repeats

Using NMR methods, we have resolved the i-motif structures formed by d(AACCCC) and by d(CCCCAA), two versions of the DNA sequence repeated in the telomeric regions of the C-rich strand of tetrahymena chromosomes. Both oligonucleotides form fully symmetrical i-motif tetramers built by intercalation of two hemiprotonated duplexes containing four C•C⁺ pairs. The structures are extremely stable. In the tetramer of d(AACCCC), the outermost C•C⁺ pairs are formed by the cytidines of the 5′ ends of the cytidine tracts. A2 forms an A2•A2 (H6trans–N7) pair stacked to C3•C3⁺ and cross-strand stacked to A1. At 0°C, the lifetimes of the hemiprotonated pairs range from 1 ms for the outermost pair to ~1 h for the innermost pairs. The tetramer of d(CCCCAA) adopts two distinct intercalation topologies in slow conformational exchange. One, whose outermost C•C⁺ pairs are built by the cytidines of the 5′ end and the other by those of the 3′ end. In both topologies, the adenosine bases are fairly well stacked to the adjacent C•C⁺ pairs. They are not paired but form symmetrical pseudo-pairs with their H6cis amino proton and N1 nitrogen pointing towards each other.     

Structure, internal motions and association-dissociation kinetics of the i-motif dimer of d(5mCCTCACTCC)

At slightly acidic pH, the association of two d(5mCCTCACTCC) strands results in the formation of an i-motif dimer. Using NMR methods, we investigated the structure of [d(5mCCTCACTCC)]₂, the internal motion of the base pairs stacked in the i-motif core, the dimer formation and dissociation kinetics versus pH. The excellent resolution of the ¹H and ³¹P spectra provided the determination of dihedral angles, which together with a large set of distance restraints, improve substantially the definition of the sugar-phosphate backbone by comparison with previous NMR studies of i-motif structures. [d(5mCCTCACTCC)]₂ is built by intercalation of two symmetrical hairpins held together by six symmetrical C•C⁺ pairs and by pair T7•T7. The hairpin loops that are formed by a single residue, A5, cross the narrow grooves on the same side of the i-motif core. The base pair intercalation order is C9•C9⁺/5mC1•5mC1⁺/C8•C8⁺/C2•C2⁺/T7.T7/C6•C6⁺/C4•C4⁺. The T3 bases are flipped out in the wide grooves. The core of the structure includes four long-lived pairs whose lifetimes at 15°C range from 100 s (C8•C8⁺) to 0.18 s (T7•T7). The formation rate and the lifetime of [d(5mCCTCACTCC)]₂ were measured between pH 6.8 and 4.8. The dimer formation rate is three to four magnitude orders slower than that of a B-DNA duplex. It depends on pH, as it must occur for a bimolecular process involving non cooperative association of neutral and protonated residues. In the range of pH investigated, the dimer lifetime, 500 s at 0°C, pH 6.8, varies approximately as 10^−pH.     

Thermodynamic examination of the pyrophosphate sensor helix in the thiamine pyrophosphate riboswitch

Riboswitches are functional mRNA that control gene expression. Thiamine pyrophosphate (TPP) binds to thi-box riboswitch RNA and allosterically inhibits genes that code for proteins involved in the biosynthesis and transport of thiamine. Thiamine binding to the pyrimidine sensor helix and pyrophosphate binding to the pyrophosphate sensor helix cause changes in RNA conformation that regulate gene expression. Here we examine the thermodynamic properties of the internal loop of the pyrophosphate binding domain by comparing the wild-type construct (RNA WT) with six modified 2 × 2 bulged RNA and one 2 × 2 bulged DNA. The wild-type construct retains five conserved bases of the pyrophosphate sensor domain, two of which are in the 2 × 2 bulge (C65 and G66). The RNA WT construct was among the most stable (ΔG°₃₇ = −7.7 kcal/mol) in 1 M KCl at pH 7.5. Breaking the A•G mismatch of the bulge decreases the stability of the construct ∼0.5–1 kcal/mol, but does not affect magnesium binding to the RNA WT. Guanine at position 48 is important for RNA–Mg²⁺ interactions of the TPP-binding riboswitch at pH 7.5. In the presence of 9.5 mM magnesium at pH 5.5, the bulged RNA constructs gained an average of 1.1 kcal/mol relative to 1 M salt. Formation of a single A+•C mismatch base pair contributes about 0.5 kcal/mol at pH 5.5, whereas two tandem A+•C mismatch base pairs together contribute about 2 kcal/mol.     

Base pair mismatches and carcinogen-modified bases in DNA: an NMR study of G.T and G.O4meT pairing in dodecanucleotide duplexes

High-resolution two-dimensional NMR studies have been completed on the self-complementary d(C-G-C-G-A-G-C-T-T-G-C-G) duplex (designated G.T 12-mer) and the self-complementary d(C-G-C-G-A-G-C-T-O4meT-G-C-G) duplex (designated G.O4meT 12-mer) containing G.T and G.O4meT pairs at identical positions four base pairs in from either end of the duplex. The exchangeable and nonexchangeable proton resonances have been assigned from an analysis of two-dimensional nuclear Overhauser enhancement (NOESY) spectra for the G.T 12-mer and G.O4meT 12-mer duplexes in H2O and D2O solution. The guanosine and thymidine imino protons in the G.T mismatch resonate at 10.57 and 11.98 ppm, respectively, and exhibit a strong NOE between themselves and to imino protons of flanking base pairs in the G.T 12-mer duplex. These results are consistent with wobble pairing at the G.T mismatch site involving two imino proton-carbonyl hydrogen bonds as reported previously [Hare, D. R., Shapiro, L., & Patel, D. J. (1986) Biochemistry 25, 7445-7456]. In contrast, the guanosine imino proton in the G.O4meT pair resonates at 8.67 ppm. The large upfield chemical shift of this proton relative to that of the imino proton resonance of G in the G.T mismatch or in G.C base pairs indicates that hydrogen bonding to O4meT is either very weak or absent. This guanosine imino proton has an NOE to the OCH3 group of O4meT across the pair and NOEs to the imino protons of flanking base pairs.(ABSTRACT TRUNCATED AT 250 WORDS)     

NMR structure and Mg2+ binding of an RNA segment that underlies the L7/L12 stalk in the E.coli 50S ribosomal subunit

Helix 42 of Domain II of Escherichia coli 23S ribosomal RNA underlies the L7/L12 stalk in the ribosome and may be significant in positioning this feature relative to the rest of the 50S ribosomal subunit. Unlike the Haloarcula marismortui and Deinococcus radiodurans examples, the lower portion of helix 42 in E.coli contains two consecutive G•A oppositions with both adenines on the same side of the stem. Herein, the structure of an analog of positions 1037–1043 and 1112–1118 in the helix 42 region is reported. NMR spectra and structure calculations support a cis Watson–Crick/Watson–Crick (cis W.C.) G•A conformation for the tandem (G•A)₂ in the analog and a minimally perturbed helical duplex stem. Mg²⁺ titration studies imply that the cis W.C. geometry of the tandem (G•A)₂ probably allows O6 of G20 and N1 of A4 to coordinate with a Mg²⁺ ion as indicated by the largest chemical shift changes associated with the imino group of G20 and the H8 of G20 and A4. A cross-strand bridging Mg²⁺ coordination has also been found in a different sequence context in the crystal structure of H.marismortui 23S rRNA, and therefore it may be a rare but general motif in Mg²⁺ coordination.     

New insights into Hoogsteen base pairs in DNA duplexes from a structure-based survey

Hoogsteen (HG) base pairs (bps) provide an alternative pairing geometry to Watson–Crick (WC) bps and can play unique functional roles in duplex DNA. Here, we use structural features unique to HG bps (syn purine base, HG hydrogen bonds and constricted C1′–C1′ distance across the bp) to search for HG bps in X-ray structures of DNA duplexes in the Protein Data Bank. The survey identifies 106 A•T and 34 G•C HG bps in DNA duplexes, many of which are undocumented in the literature. It also uncovers HG-like bps with syn purines lacking HG hydrogen bonds or constricted C1′–C1′ distances that are analogous to conformations that have been proposed to populate the WC-to-HG transition pathway. The survey reveals HG preferences similar to those observed for transient HG bps in solution by nuclear magnetic resonance, including stronger preferences for A•T versus G•C bps, TA versus GG steps, and also suggests enrichment at terminal ends with a preference for 5′-purine. HG bps induce small local perturbations in neighboring bps and, surprisingly, a small but significant degree of DNA bending (∼14°) directed toward the major groove. The survey provides insights into the preferences and structural consequences of HG bps in duplex DNA.     

Thermodynamics of DNA duplexes with adjacent G.A mismatches.

The sequence 5'-d(ATGAGCGAAT) forms a very stable self-complementary duplex with four G.A mismatch base pairs (underlined) out of ten total base pairs [Li et al. (1991) Proc. Natl. Acad. Sci. U.S.A. 88, 26-30]. The conformation is in the general B-family and is stabilized by base-pair hydrogen bonding of an unusual type, by favorable base dipole orientations, and by extensive purine-purine stacking at the mismatched sites. We have synthesized 13 decamers with systematic variations in the sequence above to determine how the flanking sequences, the number of G.A mismatches, and the mismatch sequence order (5'-GA-3' or 5'-AG-3') affect the duplex stability. Changing A.T to G.C base pairs in sequences flanking the mismatches stabilizes the duplexes, but only to the extent observed with B-form DNA. The sequence 5'-pyrimidine-GA-purine-3', however, is considerably more stable than 5'-purine-GA-pyrimidine-3'. The most stable sequences with two pairs of adjacent G.A mismatches have thermodynamic parameters for duplex formation that are comparable to those for fully Watson-Crick base-paired duplexes. Similar sequences with single G.A pairs are much less stable than sequences with adjacent G.A mismatches. Reversing the mismatch order from 5'-GA-3' to 5'-AG-3' results in an oligomer that does not form a duplex. These results agree with predictions from the model derived from NMR and molecular mechanics and indicate that the sequence 5'-pyrimidine-GA-purine-3' forms a stable conformational unit that fits quite well into a B-form double helix.(ABSTRACT TRUNCATED AT 250 WORDS)     

Mismatching base-pair dependence of the kinetics of DNA-DNA hybridization studied by surface plasmon fluorescence spectroscopy

Two single-stranded DNAs consisting of complementary base pairs except for one mismatching base pair (MM1) can form double-stranded DNA by molecular recognition. This type of duplex is not as stable as that formed by MM0. In order to add to a better understanding of the physical mechanism of the hybridization and dissociation processes at sensor (chip) surfaces, we studied the kinetics of the MM1 hybridization by surface plasmon fluorescence spectroscopy. Target DNA strands labelled with a fluorescent molecule Cy5 at the 5′ end and hybridizing with the surface-attached probe DNA can be excited by the strong optical field of a surface plasmon resonance mode. The emitted fluorescence can be detected with high sensitivity. The affinity of a duplex was found to depend on the chemical nature, i.e. G–G, G–T etc., and on the position of the mismatching base pair along the 15mer duplex.     

Mismatch base pair detection by fluorescence spectral change upon addition of metal cation--toward efficient analysis of single nucleotide polymorphism

Addition of mercury (II) cation to fluorescent-labeled duplex involving a T:T mismatch base pair and silver (I) cation to fluorescent-labeled duplex involving a C:C mismatch base pair significantly changed the fluorescence intensity, but no significant change in the fluorescence intensity was observed for duplexes involving the other base pairs. The fluorescence spectral change upon addition of the metal cation can discriminate T:T and C:C mismatch base pairs from the other base pairs. Our results certainly support the idea that the fluorescence spectral change upon addition of the metal cation could be a convenient strategy for the mismatch base pair detection by the heteroduplex analysis, and may eventually lead to progress in single nucleotide polymorphism genotyping.     

Studies of translational misreading in vivo show that the ribosome very efficiently discriminates against most potential errors

Protein synthesis must rapidly and repeatedly discriminate between a single correct and many incorrect aminoacyl-tRNAs. We have attempted to measure the frequencies of all possible missense errors by tRNA, tRNA and tRNA. The most frequent errors involve three types of mismatched nucleotide pairs, U•U, U•C, or U•G, all of which can form a noncanonical base pair with geometry similar to that of the canonical U•A or C•G Watson–Crick pairs. Our system is sensitive enough to measure errors at other potential mismatches that occur at frequencies as low as 1 in 500,000 codons. The ribosome appears to discriminate this efficiently against any pair with non-Watson–Crick geometry. This extreme accuracy may be necessary to allow discrimination against the errors involving near Watson–Crick pairing.