An NMR Study of the Conformation and Thermodynamics of the Circular Dumbbell d<pCGC-TT-GCG-TT>

Abstract

The circular DNA decamer 5′-d<pCGC-TT-GCG-TT>3′ is studied in solution by means of NMR spectroscopy. At low temperature the molecule adopts a dumbbell structure with three Watson-Crick C-G base pairs and two two-residue loops in opposite parts of the molecule. On raising the temperature another conformer appears, in which the closing C-G base pair in the 5′-GTTC-3′ loop is disrupted, whereas the opposite 5′-CTTG-3′ loop remains stable. The two conformers are in slow equilibrium over a limited temperature range.     

Solution structure of a five-adenine bulge loop within a DNA duplex.

The three-dimensional solution structure of a DNA molecule of the sequence 5'-d(GCATCGAAAAAGCTACG)-3' paired with 5'-d(CGTAGCCGATGC)-3' containing a five-adenine bulge loop (dA(5)-bulge) between two double helical stems was determined by 2D (1)H and (31)P NMR, infrared, and Raman spectroscopy. The DNA in both stems adopt a classical B-form double helical structure with Watson-Crick base pairing and C2'-endo sugar conformation. In addition, the two dG/dC base pairs framing the dA(5)-bulge loop are formed and are stable at least up to 30 degrees C. The five adenine bases of the bulge loop are localized at intrahelical positions within the double helical stems. Stacking on the double helical stem is continued for the first four 5'-adenines in the bulge loop. The total rise (the height) of these four stacked adenines roughly equals the diameter of the double helical stem. The stacking interactions are broken between the last of these four 5'-adenines and the fifth loop adenine at the 3'-end. This 3'-adenine partially stacks on the other stem. The angle between the base planes of the two nonstacking adenines (A10 and A11) in the bulge loop reflects the kinking angle of the global DNA structure. The neighboring cytosines opposite the dA(5)-bulge (being parts of the bulge flanking base pairs) do not stack on one another. This disruption of stacking is characterized by a partial shearing of these bases, such that certain sequential NOEs for this base step are preserved. In the base step opposite the loop, an extraordinary hydrogen bond is observed between the phosphate backbone of the 5'-dC and the amino proton of the 3'-dC in about two-thirds of the conformers. This hydrogen bond probably contributes to stabilizing the global DNA structure. The dA(5)-bulge induces a local kink into the DNA molecule of about 73 degrees (+/-11 degrees ). This kinking angle and the mutual orientation of the two double helical stems agree well with results from fluorescence resonance energy transfer measurements of single- and double-bulge DNA molecules.     

Slow conformational exchange in DNA minihairpin loops: A conformationalstudy of the circular dumbbell d〈pCGC-TT-GCG-TT〉

In recent years various examples of highly stable two-residue hairpin loops (miniloops) in DNA have been encountered. As the detailed structure and stability of miniloops appear to be determined not only by the nature and sequence of the two bases in the loop, but also by the closing base pair, it is desirable to carry out in-depth studies of especially designed small model DNA compounds. Therefore, a circular DNA dumbbell-like molecule is tailored to consist of a stem of three Watson–Crick base pairs, flanked on each side by a minihairpin loop. The resulting circular DNA decanter 5′-d〈pCGC- TT-GCG- TT〉 -3′ ( I ) is studied in solution by means of nmr spectroscope. At a temperature of 269 K the molecule occurs in a 50/50 mixture of two dumbbell structures (denoted L2L2 and L2L4). L2L2 contains three Watson–Crick C-G base pairs and two two-residue loops (H2-family type) in opposite parts of the molecule. On raising the temperature from 269 to 314 K. The L2L4 conformer becomes increasingly dominant (95% at 314 K). This conformer has a partially disrupted closing G-C base pair in the 5′-GTTC-3′ loop with only one remaining solvent-accessible hydrogen bond between NH_α of the cytosine C(1) and O6 of the guanine G(8), whereas the opposite 5′-CTTG-3′ loop remains stable. The disruption of the C(1)-G(8) base pair in the L2L4 form is correlated with the presence of a syn orientation for the C(1) base at the 5′-3′ loop-stem junction in the 5′-GTTC-3′ loop. The two conformers. L2L2 and L2L4, occur in slow equilibrium (2–20 s^?1). Moderate line broadening of specific ¹H, ¹³C, and ³¹P resonances of residues C(1), G(8), T(9), and T(10) at low temperatures, due to chemical exchange between L2L2 and L2L4, show that the interconversion from an anti to syn conformer in residue C(1) has a small local effect on the structure of the dumbbell. T₁ relaxation measurements, chemical-shift considerations, and complete hand-shape calculations of the exchange process of the G(8) imino proton reveal a possibility for the existence of multiconformational slates in the anti–syn equilibrium. © 1995 John Wiley & Sons, Inc.     

Complementary addressed modification of yeast tRNA Val 1 with alkylating derivative of d(pC-G)-A. The positions of the alkylated nucleotides and the course of the alkylation in the complex.

Yeast tRNA Val 1 alkylation with 2', 3'-O-4-(N-2-chloroethyl-N-methylamino) benzylidene d(pC-G)-A proceeds at 20 degrees - 30 degrees C in the complementary complexes which are formed by d(pC-G)-A greater than RC1 binding to 3 sequences of tRNA Val 1 : psi-C-G58 in the T loop, C-G40 at the 3'-side of the anticodon loop and C-G18 in the D loop. The reaction in the complexes results in A53, I35, and psi 13 alkylation to form beta-/N-methyl-N-(formylphenyl 17 amino/ethyl-tRNA Val 1 with the relative rate constants of the alkylation that are 3 or 2 orders of magnitude higher than that for the alkylation without a complex formation. It is the third nucleotide from the 5'-terminus of the binding site of the modifying agent that is subjected to alkylation in the t RNA Val 1. The course of the alkylation does not depend on the possible base pairing of the 3'-terminal nucleotide of the reagent. The extent of the reagent binding and the relative rate constants of the alkalytion in the complexes indicate the following order of the complex stability: (psi-C-G58) greater than (CO-G40) approximately (C-G18) at 20 degrees and (psi-C-G58) greater than (C-G40) greater than (C-G18) at 30 degrees.     

Multinuclear NMR studies of DNA hairpins. 1. Structure and dynamics of d(CGCGTTGTTCGCG)

The solution structure of the hairpin formed by d(CGCGTTGTTCGCG) has been examined in detail by a wide variety of NMR techniques. The hairpin was characterized by proton NMR to obtain interproton distances and torsion angle information. An energy-minimized model was constructed that is consistent with these data. The hairpin consists of a B-DNA stem of four C-G base pairs and a loop region consisting of five unpaired bases. Three bases in the 5' of the loop are stacked over the 3' end of the stem, and the other two bases in the 3' of the loop are stacked over the 5' end of the stem. The phosphorus NMR spectrum revealed a phosphate in the stem region with an unusual conformation, and two phosphates, P9 and P10, were found to undergo intermediate exchange between conformations. The hairpin was also synthesized with a carbon-13 label in each of the thymidine C6 carbons, and relaxation measurements were performed to determine the extent of internal motions in the loop region. The loop bases are more flexible than the stem bases and exhibit subnanosecond motions with an amplitude corresponding to diffusion in a cone of approximately 30 degrees.     

Conformational feasibility of a hairpin with two purines in the loop. 5'-d-GGTACIAGTACC-3'

Structural feasibility and conformational requirements for the sequence 5'-d-GGTACIAGTACC-3' to adopt a hairpin loop with I6 and A7 in the loop are studied. It is shown that a hairpin loop containing only two nucleotides can readily be formed without any unusual torsional angles. Stacking is continued on the 5'-side of the loop, with the I6 stacked upon C5. The base A7, on the 3'-side of the loop, can either be partially stacked with I6 or stick outside without stacking. Loop closure can be achieved for both syn and anti conformations of the glycosidic torsions for G8 while maintaining the normal Watson-Crick base pairing with the opposite C5. All torsional angles in the stem fall within the standard B-family of DNA helical structures. The phosphodiesters of the loop have trans,trans conformations. Loop formation might require the torsion about the C4'-C5' bond of G8 to be trans as opposed to the gauche+ observed in B-DNA. These results are discussed in relation to melting temperature studies [Howard et al. (1991) Biochemistry (preceding paper in this issue)] that suggest the formation of very stable hairpin structures for this sequence.     

Thermodynamics of melting of the circular dumbbell d〈pCGC-TT-GCG-TT〉

The conformational behavior of DNA minihairpin loops is sensitive to the directionality of the base pair that closes the loop. Especially tailored circular dumbbells, consisting of a stem of three Watson–Crick base pairs capped on each side with a minihairpin loop, serve as excellent model compounds by means of which deeper insight is gained into the relative stability and melting properties of hairpin loops that differ only in directionality of the closing pair: C-G vs G-C. For this reason the thermodynamic properties of the circular DNA decamers 5′-d〈pCGC-TT-GCG-TT〉-3′( I ) and reference compounds 5′-d〈pGGC-TT-GCC-TT≤-3′( II ) and 5′-d(GCG-TC-CGC)-3′( III ) are studied by means of nmr spectroscopy. Molecules I and II adopt dumbbell structures closed on both sides by a two-membered hairpin hop. At low temperature I consists of a mixture of two slowly exchanging forms, denoted L2L2 and L2L4 . The low-temperature L2L2 form is the fully intact minihairpin structure with three Watson–Crick C-G base pairs. The high-temperature form, L2L4 ,contains a partially disrupted closing G-C base pair in the 5′-GTTC-3′ loop, with the cytosine base placed in a syn orientation. The opposite 5′-CTTG-3′ loop remains stable. A study of the noncircular hairpin structure III shows similar conformational behavior for the 5′-GTTC-3′ loop as found in I a syn orientation for C(6) and two slowly exchanging imino proton signals for G(3). The melting point T_m of II was estimated to lie above 365 K. The T_m value of the duplex stem and the 5′-CTTG-3′ loop of the L2L4 form ofIis 352 ± 2 K. The ΔH° is calculated as ?89 ± 10 kJ/mol. The T_m value determined for the individual residues of the 5′-GTTC-3′ loop lies 4°–11° lower. The enthalpy ΔH° of melting the thymine residues in the 5′-GTTC-3′ loop is calculated to be -61± 7 kJ/mol. Thermodynamic data of the equilibrium between the slowly exchanging two- and four-membered loop conformers of I reveal an upper limit for ΔH° of +30 kJ/mol in going from a two-memberedto a four-membered loop, in agreement with the enthalpy difference of +28 k.j/mol between the two loops at the T_m midpoint. For hairpin III the upper limit for ΔH° going from a two-membered to a four-membered loop amounts to ±21 kJ/mol. The mutual exchange rate between the L2 and L4 form in III is estimated as 13.6 s^?1. Our results clearly suggest that small four-way DNA junctions(model for immobilized Holliday junctions) can be designed that consist of a single DNA strandthat features -CTTG-caps on three of the four arms of the junction. © 1995 John Wiley & Sons, Inc.     

The solution structure of a d[C(TTCG)G] DNA hairpin and comparison to the unusually stable RNA analogue.

The solution structure of the DNA analogue of the unusually stable r[C(UUCG)G] RNA hairpin, 5'-d[GGA-C(TTCG)GTCC]-3', has been determined by NMR spectroscopy, and its structure has been compared to that of the RNA molecule. The RNA molecule is compact and rigid with a highly structured loop. However, the DNA molecule is much less structured. The DNA hairpin contains a B-form stem of four base pairs. The terminal base pair frays, and the 3'-terminal nucleotides, C11 and C12, are in equilibrium between 2'-endo and 3'-endo conformations. Unlike the RNA loop, the DNA loop contains no syn nucleotides, and there is no evidence for base-base or base-phosphate hydrogen bonding in the loop. The loop is flexible, and reveals no specific internucleotide interactions.     

Dynamics and stability of GCAA tetraloops with 2-aminopurine and purine substitutions

The contributions of various interactions in the GGCGCAAGCC hairpin containing a GCAA tetraloop were studied by computer simulations using the substitutions of functional groups. The guanosine (G) in the first tetraloop position or in the C-G closing base pair was replaced by 2-aminopurine (AP), and the individual tetraloop's adenosines (A) were replaced by purine (PUR). These substitutions eliminated particular hydrogen bonds thought to stabilize the GCAA tetraloop. For each substitution, molecular dynamics (MD) simulations were carried out in an aqueous solution with sodium counterions, using the CHARMM27 force field. The MD simulations showed that the substitutions in the first (G-->AP) and the third (A-->PUR) position of the GCAA tetraloop did not significantly influence the conformation of the hairpin. A long-lived bridging water molecule observed in the GCAA loop was present in both modified loops. The substitutions made in the last loop position (A-->PUR) or in the C-G base pair closing the tetraloop (G-->AP) to some extent influenced the loop structure and dynamics. These loops did not display the long-lived bridging water molecules. When the second A in the GCAA loop was replaced by PUR, the first A in the loop was observed in the anti or in the syn orientation about the glycosyl bond. The G to AP substitution in C-G base pair led to a change of their arrangement from the Watson-Crick to wobble. The MD simulations of the hairpin with C-AP wobble closing base pair showed increased conformational dynamics of the hairpin. The changes of hairpin formation free energy associated with the substitutions of individual bases were calculated by the free energy perturbation method. Our theoretical estimates suggest a larger destabilization for the G to AP substitutions in GCAA loop than for the substitutions of individual A's by PUR, which is in accordance with experimental tendency. The calculations predicted a similar free energy change for G to AP substitutions in the GCAA tetraloop and in the C-G closing base pair.     

Hairpin and duplex formation of the DNA octamer d(m5C-G-m5C-G-T-G-m5C-G) in solution. An NMR study.

The partly self-complementary DNA octamer d(m5C-G-m5C-G-T-G-m5C-G) was investigated by NMR spectroscopy in solution. It is demonstrated that this peculiar DNA fragment, under suitable conditions of concentration, salt and temperature, exclusively prefers to adopt a monomeric hairpin form with a stem of three Watson-Crick type base pairs and a loop of two residues. At high single strand concentration (8 mM DNA) and low temperature (i.e. below 295 K) the hairpin occurs in slow equilibrium with a B-dimer structure. At high ionic strength (greater than or equal to 100 mM Na+) and/or in the presence of methanol a third species appears, which is assigned to a Z-like dimer. In the B form, as well as in the Z dimer, the two central base pairs form G.T wobble pairs with the bases as major tautomers.     

The consensus sequence for self-catalyzed site-specific G residue depurination in DNA

The sequence variation tolerated within the stem-loop-forming genomic consensus sequence for self-catalyzed site-specific depurination of G residues is explored. The variation in self-depurination kinetics with sequence changes in the loop residues and stem base pairs, as well as with pH, provides insights into the self-catalytic mechanism. The observations suggest that self-catalyzed depurination of the 5' G residue of the loop consensus sequence 5'-G(T/A)GG-3' probably involves formation of some intraloop hydrogen-bonded base pair with the 3'-terminal G residue; although the electronic structure of both these G residues is retained, their 2-amino substituents are not critical for that interaction. The strong dependence of the self-depurination kinetics on stem stability suggests that the lifetime of some strained form of the loop is controlled by the integrity of the stem. In addition to the effects of length and base pair sequence on stem stability, there is a base pair requirement at the base of the loop: self-depurination is suppressed by 5'-C·G-3', 5'-A·T-3', or a mismatch but is most favored by 5'T·A3' and less so by 5'-G·C-3'. The occurrence in T and G of a similarly located carbonyl capable of hydrogen-bonding to the water molecule required for glycosyl bond hydrolysis may explain this sequence requirement. In toto, the more complete definition of the consensus sequence provided by this investigation enables a more accurate estimation of their number in the human genome and their distribution among different genes.     

RecA homology search is promoted by mechanical stress along the scanned duplex DNA

A RecA-single-stranded DNA (RecA-ssDNA) filament searches a genome for sequence homology by rapidly binding and unbinding double-stranded DNA (dsDNA) until homology is found. We demonstrate that pulling on the opposite termini (3' and 5') of one of the two DNA strands in a dsDNA molecule stabilizes the normally unstable binding of that dsDNA to non-homologous RecA-ssDNA filaments, whereas pulling on the two 3', the two 5', or all four termini does not. We propose that the 'outgoing' strand in the dsDNA is extended by strong DNA-protein contacts, whereas the 'complementary' strand is extended by the tension on the base pairs that connect the 'complementary' strand to the 'outgoing' strand. The stress resulting from different levels of tension on its constitutive strands causes rapid dsDNA unbinding unless sufficient homology is present.     

Synthesis and properties of oligodeoxyribonucleotides containing a mutagenic base, N4-aminocytosine.

Oligodeoxyribonucleotides containing a mutagenic base analog, N4-aminocytosine, 5'-AATTGC(am)AATT-3' and 5'-AATTAC(am)AATT-3' (C(am); N4-aminocytosine) were prepared by chemical modification of 5'-AATTGCAATT-3' and 5'-AATTACAATT-3', respectively. The values of Tm were 29 degrees C for 5'-AATTGC(am)AATT-3' and 32 degrees C for 5'-AATTGCAATT-3'. In contrast, no melting was observed for 5'-AATTAC(am)AATT-3' and 5'-AATTACAATT-3'. These data show that the stability of C(am)-purine paris is C(am)-G > C(am)-A and that C(am)-G is less stable than C-G. This property is consistent with the incorporation specificity of N4-amino-dCTP during DNA synthesis in vitro.     

Dynamics and Stability of GCAA Tetraloops with 2-Aminopurine and Purine Substitutions

Abstract

The contributions of various interactions in the GGCGCAAGCC hairpin containing a GCAA tetraloop were studied by computer simulations using the substitutions of functional groups. The guanosine (G) in the first tetraloop position or in the C-G closing base pair was replaced by 2-aminopurine (AP), and the individual tetraloop's adenosines (A) were replaced by purine (PUR). These substitutions eliminated particular hydrogen bonds thought to stabilize the GCAA tetraloop. For each substitution, molecular dynamics (MD) simulations were carried out in an aqueous solution with sodium counterions, using the CHARMM27 force field. The MD simulations showed that the substitutions in the first (G→AP) and the third (A→PUR) position of the GCAA tetraloop did not significantly influence the conformation of the hairpin. A long-lived bridging water molecule observed in the GCAA loop was present in both modified loops. The substitutions made in the last loop position (A→PUR) or in the C-G base pair closing the tetraloop (G→AP) to some extent influenced the loop structure and dynamics. These loops did not display the long- lived bridging water molecules. When the second A in the GCAA loop was replaced by PUR, the first A in the loop was observed in the anti or in the syn orientation about the gly- cosyl bond. The G to AP substitution in C-G base pair led to a change of their arrangement from the Watson-Crick to wobble. The MD simulations of the hairpin with C-AP wobble closing base pair showed increased conformational dynamics of the hairpin. The changes of hairpin formation free energy associated with the substitutions of individual bases were calculated by the free energy perturbation method. Our theoretical estimates suggest a larger destabilization for the G to AP substitutions in GCAA loop than for the substitutions of individual A's by PUR, which is in accordance with experimental tendency. The calculations predicted a similar free energy change for G to AP substitutions in the GCAA tetraloop and in the C-G closing base pair.     

Hairpin and duplex formation in DNA fragments CCAATTTTGG, CCAATTTTTTGG, and CCATTTTTGG: a proton NMR study

Three DNA fragments, CCAATTTTGG (1), CCAATTTTTTGG (2), and CCATTTTTGG (3), were studied by proton NMR spectroscopy in aqueous solution. All these oligodeoxyribonucleotides contain common sequences at the 5' and 3' ends (5'-CCA and TGG-3'). 2 as well as 3 forms only hairpin structures with four unpaired thymidylyl units, four and three base pair stems, respectively, in neutral solution under low and high NaCl concentrations. At high salt concentration the oligomer 1 forms a duplex structure with -TT- internal loop. On the other hand, the same oligomer forms a stable hairpin structure at low salt and low strand concentrations at pH 7. The hairpin structure of 1 has a stem containing only three base pairs (CCA.TGG) and a loop containing four nucleotides (-ATTT-) that includes a dissociated A.T base pair. The two secondary structures of 1 coexist in an aqueous solution containing 0.1 M NaCl, at pH 7. The equilibrium shifts to the hairpin side when the temperature is raised. The stabilities and base-stacking modes of all three oligonucleotides in two different structures are reported.     

Experimental and computational studies of the G[UUCG]C RNA tetraloop

In prokaryotic ribosomal RNAs, most UUCG tetraloops are closed by a C-G base-pair. However, this preference is greatly reduced in eukaryotic rRNA species where many UUCG tetraloops are closed by G-C base-pairs. Here, biophysical properties of the C[UUCG]G and G[UUCG]C tetraloops are compared, using experimental and computational methods. Thermal denaturation experiments are used to derive thermodynamic parameters for the wild-type G[UUCG]C tetraloop and variants containing single deoxy substitutions in the loop. A comparison with analogous experiments on the C[UUCG]G motif shows that the two RNA species exhibit similar patterns in response to the substitutions, suggesting that their loop structures are similar. This conclusion is supported by NMR data that suggest that the essential UUCG loop structure is maintained in both tetraloops. However, NMR results show that the G[UUCG]C loop structure is disrupted prior to melting of the stem; this behavior is in contrast to the two-state behavior of the C[UUCG]G molecule. Stochastic dynamics simulations using the GB/SA continuum solvation model, run as a function of temperature, show rare conformational transitions in several G[UUCG]C simulations. These results lead to the conclusion that substitution of a G-C for a C-G closing base-pair increases the intrinsic flexibility of the UUCG loop.     

Base Pairing between U3 Small Nucleolar RNA and the 5′ End of 18S rRNA Is Required for Pre-rRNA Processing

The loop of a stem structure close to the 5' end of the 18S rRNA is complementary to the box A region of the U3 small nucleolar RNA (snoRNA). Substitution of the 18S loop nucleotides inhibited pre-rRNA cleavage at site A(1), the 5' end of the 18S rRNA, and at site A(2), located 1.9 kb away in internal transcribed spacer 1. This inhibition was largely suppressed by a compensatory mutation in U3, demonstrating functional base pairing. The U3-pre-rRNA base pairing is incompatible with the structure that forms in the mature 18S rRNA and may prevent premature folding of the pre-rRNA. In the Escherichia coli pre-rRNA the homologous region of the 16S rRNA is also sequestered, in that case by base pairing to the 5' external transcribed spacer (5' ETS). Cleavage at site A(0) in the yeast 5' ETS strictly requires base pairing between U3 and a sequence within the 5' ETS. In contrast, the U3-18S interaction is not required for A(0) cleavage. U3 therefore carries out at least two functionally distinct base pair interactions with the pre-rRNA. The nucleotide at the site of A(1) cleavage was shown to be specified by two distinct signals; one of these is the stem-loop structure within the 18S rRNA. However, in contrast to the efficiency of cleavage, the position of A(1) cleavage is not dependent on the U3-loop interaction. We conclude that the 18S stem-loop structure is recognized at least twice during pre-rRNA processing.     

The role of RNA structure in the interaction of U1A protein with U1 hairpin II RNA

The N-terminal RNA Recognition Motif (RRM1) of the spliceosomal protein U1A interacting with its target U1 hairpin II (U1hpII) has been used as a paradigm for RRM-containing proteins interacting with their RNA targets. U1A binds to U1hpII via direct interactions with a 7-nucleotide (nt) consensus binding sequence at the 5' end of a 10-nt loop, and via hydrogen bonds with the closing C-G base pair at the top of the RNA stem. Using surface plasmon resonance (Biacore), we have examined the role of structural features of U1hpII in binding to U1A RRM1. Mutational analysis of the closing base pair suggests it plays a minor role in binding and mainly prevents "breathing" of the loop. Lengthening the stem and nontarget part of the loop suggests that the increased negative charge of the RNA might slightly aid association. However, this is offset by an increase in dissociation, which may be caused by attraction of the RRM to nontarget parts of the RNA. Studies of a single stranded target and RNAs with untethered loops indicate that structure is not very relevant for association but is important for complex stability. In particular, breaking the link between the stem and the 5' side of the loop greatly increases complex dissociation, presumably by hindering simultaneous contacts between the RRM and stem and loop nucleotides. While binding of U1A to a single stranded target is much weaker than to U1hpII, it occurs with nanomolar affinity, supporting recent evidence that binding of unstructured RNA by U1A has physiological significance.