Crystal structure of a 14 bp RNA duplex with non-symmetrical tandem GxU wobble base pairs

Adjacent GxU wobble base pairs are frequently found in rRNA. Atomic structures of small RNA motifs help to provide a better understanding of the effects of various tandem mismatches on duplex structure and stability, thereby providing better rules for RNA structure prediction and validation. The crystal structure of an RNA duplex containing the sequence r(GGUAUUGC-GGUACC)2 has been solved at 2.1 A resolution using experimental phases. Novel refinement strategies were needed for building the correct solvent model. At present, this is the only short RNA duplex structure containing 5'-U-U-3'/3'-G-G-5' non-symmetric tandem GxU wobble base pairs. In the 14mer duplex, the six central base pairs are all displaced away from the helix axis, yielding significant changes in local backbone conformation, helix parameters and charge distribution that may provide specific recognition sites for biologically relevant ligand binding. The greatest deviations from A-form helix occur where the guanine of a wobble base pair stacks over a purine from the opposite strand. In this vicinity, the intra-strand phosphate distances increase significantly, and the major groove width increases up to 3 A. Structural comparisons with other short duplexes containing symmetrical tandem GxU or GxT wobble base pairs show that nearest-neighbor sequence dependencies govern helical twist and the occurrence of cross-strand purine stacks.     

Crystal structure of an RNA helix recognized by a zinc-finger protein: an 18-bp duplex at 1.6 A resolution

The crystal structure of the 19-mer RNA, 5'-GAAUGCCUGCGAGCAUCCC-3' has been determined from X-ray diffraction data to 1.6 A resolution by the multiwavelength anomalous diffraction method from crystals containing a brominated uridine. In the crystal, this RNA forms an 18-mer self-complementary double helix with the 19th nucleotide flipped out of the helix. This helix contains most of the target stem recognized by the bacteriophage Mu Com protein (control of mom), which activates translation of an unusual DNA modification enzyme, Mom. The 19-mer duplex, which contains one A.C mismatch and one A.C/G.U tandem wobble pair, was shown to bind to the Com protein by native gel electrophoresis shift assay. Comparison of the geometries and base stacking properties between Watson-Crick base pairs and the mismatches in the crystal structure suggest that both hydrogen bonding and base stacking are important for stabilizing these mismatched base pairs, and that the unusual geometry adopted by the A.C mismatch may reveal a unique structural motif required for the function of Com.     

Altered structural fluctuations in duplex RNA versus DNA: a conformational switch involving base pair opening

DNA and RNA are known to have different structural properties. In the present study, molecular dynamics (MD) simulations on a series of RNA and DNA duplexes indicate differential structural flexibility for the two classes of oligonucleotides. In duplex RNA, multiple base pairs experienced local opening events into the major groove on the nanosecond time scale, while such events were not observed in the DNA simulations. Three factors are indicated to be responsible for the base opening events in RNA: solvent-base interactions, 2'OH(n)-O4'(n+1) intra-strand hydrogen bonding, and enhanced rigid body motion of RNA at the nucleoside level. Water molecules in the major groove of RNA contribute to initiation of base pair opening. Stabilization of the base pair open state is due to a 'conformational switch' comprised of 2'OH(n)-O4'(n+1) hydrogen bonding and a rigid body motion of the nucleoside moiety in RNA. This rigid body motion is associated with decreased flexibility of the glycosyl linkage and sugar moieties in A-form structures. The observed opening rates in RNA are consistent with the imino proton exchange experiments for AU base pairs, although not for GC base pairs, while structural and flexibility changes associated with the proposed conformational switch are consistent with survey data of RNA and DNA crystal structures. The possible relevance of base pair opening events in RNA to its many biological functions is discussed.     

The crystal structure at 1.5 angstroms resolution of an RNA octamer duplex containing tandem G.U basepairs

The crystal structure of the RNA octamer, 5'-GGCGUGCC-3' has been determined from x-ray diffraction data to 1.5 angstroms resolution. In the crystal, this oligonucleotide forms five self-complementary double-helices in the asymmetric unit. Tandem 5'GU/3'UG basepairs comprise an internal loop in the middle of each duplex. The NMR structure of this octameric RNA sequence is also known, allowing comparison of the variation among the five crystallographic duplexes and the solution structure. The G.U pairs in the five duplexes of the crystal form two direct hydrogen bonds and are stabilized by water molecules that bridge between the base of guanine (N2) and the sugar (O2') of uracil. This contrasts with the NMR structure in which only one direct hydrogen bond is observed for the G.U pairs. The reduced stability of the r(CGUG)2 motif relative to the r(GGUC)2 motif may be explained by the lack of stacking of the uracil bases between the Watson-Crick and G.U pairs as observed in the crystal structure.     

C-C-A-G-G-C-m5C-T-G-G. Helical fine structure, hydration, and comparison with C-C-A-G-G-C-C-T-G-G.

The self-complementary DNA duplex C-C-A-G-G-C-m5C-T-G-G has been refined against 1.75-A x-ray diffraction data to an R value of 17.4%. In the crystal of space group P6, 10-base pair DNA fragments with characteristic sequence-related fine structure stack end to end to form long antiparallel B-type double helices. As shown by a structure analysis at lower resolution (Heinemann, U., and Alings, C. (1991) EMBO J. 10, 35-43), the overall geometry of C-C-A-G-G-C-m5C-T-G-G is similar to that of the unmethylated analog C-C-A-G-G-C-C-T-G-G despite a different crystal environment. The present high resolution structure analysis permits a detailed comparison of the two duplexes and their hydration spheres. Helical parameters are significantly correlated between both molecules, with the exception of the base pair propeller. Sugar pucker and backbone torsion angles alpha, gamma, delta, and chi show similar mean values, but their individual values deviate significantly between duplexes. In contrast, torsion angles beta, epsilon, and zeta change along the strands of both duplexes in much the same way. The effect of single-site methylation on DNA conformation appears to be small and limited to the base pairs directly involved. Methylation tends to push base pairs toward the minor groove of the helix. A regular minor groove hydration pattern involves dual hydrogen bonding of water molecules to O-4' and base atoms of C-C-A-G-G-C-m5C-T-G-G.     

Nucleic acid model building: the multiple backbone solutions associated with a given base morphology

A constrained model building procedure is used to generate nucleic acid structures of the familiar A-, B-, and Z-DNA duplexes. Attention is focused upon the multiple structural solutions associated with the arrangements of nucleic acid base pairs rather than the optimum sugar-phosphate structure. The glycosyl (chi) and sugar torsions (both the ring puckering and the exocyclic C5'-C4' (psi) torsion) are treated as independent variables and the resulting O3'...O5' distances are used as closure determinants. When such distances conform to the known geometry of phosphate chemical bonding, an intervening phosphorus atom with correct C-O-P valence angles can be located. Four sequential torsion angles--phi', omega', omega and phi--about the C3'-O3'-P-O5'-C5' bonds are then obtained as dependent variables. The resulting structures are categorized in terms of conformation, ranked in potential energy, and analyzed for torsional correlations. The numerical results are quite interesting with implications regarding nucleic acid models constructed to fit less than ideal experimental data. The multiple solutions to the problem are useful for comprehending the conformational complexities of the local sugar-phosphate backbone and for understanding the transitions between different helical forms. According to these studies, unique characterization of a nucleic acid duplex involves more than the determination of its base pair morphology, its sugar puckering preferences, or its groove binding features.     

Solution conformation of a bulged adenosine base in an RNA duplex by relaxation matrix refinement

Bulges are common structural motifs in RNA secondary structure and are thought to play important roles in RNA-protein and RNA-drug interactions. Adenosine bases are the most commonly occurring unpaired base in double helical RNA secondary structures. The solution conformation and dynamics of a 25-nucleotide RNA duplex containing an unpaired adenosine, r(GGCAGAGUGCCGC): r(GCGGCACCUGCC) have been studied by NMR spectroscopy and MORASS iterative relaxation matrix structural refinement. The results show that the bulged adenosine residue stacks into the RNA duplex with little perturbation around the bulged region. Most of the bases in the RNA duplex adopt C(3)'-endo conformation, exhibiting the N-type sugar pucker as found in the A form helices. The sugars of the bulged residue and the 5' flanking residue to it are found to exhibit C(2)'-endo conformation. None of the residues are in syn conformation.     

Stabilization factors affecting duplex formation of peptide nucleic acid with DNA.

Peptide nucleic acid (PNA) is an oligonucleotide analogue in which the sugar-phosphate backbone is replaced by an N-(2-aminoethyl)glycine unit to which the nucleobases are attached. We investigated the thermodynamic behavior of PNA/DNA hybrid duplexes with identical nearest neighbors but with different sequences and chain lengths (5, 6, 7, 8, 10, 12, and 16 mers) to reveal whether the nearest-neighbor model is valid for the PNA/DNA duplex stability. CD spectra of 6, 7, and 8 mer PNA/DNA duplexes showed similar signal, while 10, 12, and 16 mer duplexes did not. The average difference in Delta G degrees (37) for short PNA/DNA duplexes with identical nearest-neighbor pairs was only 3.5%, whereas that of longer duplexes (10, 12, and 16 mers) was 16.4%. Therefore, the nearest-neighbor model seems to be useful at least for the short PNA/DNA duplexes. Thermodynamics of PNA/DNA duplexes containing 1--3 bulge residues were also studied. While the stability of the 12 mer DNA/DNA duplex decreased as the number of bulge bases increases, the number of bulge bases in PNA/DNA unchanged the duplex stability. Thus, the influence of bulge insertion in the PNA/DNA duplexes is different from that of a DNA/DNA duplex. This might be due to the different base geometry in a helix which may potentially make hydrogen bonds in a base pair and stacking interaction unfavorable compared with DNA/DNA duplexes.     

Dynamics in the U6 RNA intramolecular stem-loop: a base flipping conformational change

The U6 RNA intramolecular stem-loop (ISL) structure is an essential component of the spliceosome and binds a metal ion required for pre-messenger RNA splicing. The metal binding internal loop region of the stem contains a partially protonated C67-(+)A79 base pair (pK(a) = 6.5) and an unpaired U80 nucleotide that is stacked within the helix at pH 7.0. Here, we determine that protonation occurs with an exchange lifetime of approximately 20 micros and report the solution structures of the U6 ISL at pH 5.7. The differences between pH 5.7 and 7.0 structures reveal that the pH change significantly alters the RNA conformation. At lower pH, U80 is flipped out into the major groove. Base flipping involves a purine stacking interaction of flanking nucleotides, inversion of the sugar pucker 5' to the flipped base, and phosphodiester backbone rearrangement. Analysis of residual dipolar couplings as a function of pH indicates that base flipping is not restricted to a local conformational change. Rather, base flipping alters the alignment of the upper and lower helices. The alternative conformations of the U6 ISL reveal striking structural similarities with both the NMR and crystal structures of domain 5 of self-splicing group II introns. These structures suggest that base flipping at an essential metal binding site is a conserved feature of the splicing machinery for both the spliceosome and group II self-splicing introns.     

The crystal structure analysis of d(CGCGAASSCGCG)2, a synthetic DNA dodecamer duplex containing four 4'-thio-2'-deoxythymidine nucleotides.

The crystal structure refinement of the synthetic dodecamer d(CGCGAASSCGCG), where S = 4'-thio-2'-deoxythymidine, has converged at R=0.201 for 2605 reflections with F > 2sigma(F) in the resolution range 8.0-2.4 A for a model consisting of the dodecamer duplex and 66 water molecules. A comparison of its structure with that of the native dodecamer d(CGCGAATTCGCG) has revealed that the major differences between the two structures is a change in the conformation of the sugar-phosphate backbone in the regions at and adjacent to the positions of the modified nucleosides. Examination of the fine structural parameters for each of the structures reveals that the thiosugars adopt a C3'-exo conformation in d(CGCGAASSCGCG), rather than the approximate C1'-exo conformation found for the analogous sugars in the structure of d(CGCGAATTCGCG). The observed differences in structure between the two duplexes may help to explain the enhanced resistance to nuclease digestion of synthetic oligonucleotides containing 4'-thio-2'-deoxynucleotides.     

A stable antiparallel cytosine-thymine base pair occurring only at the end of a duplex.

Ultraviolet hyperchromicity experiments indicate that in DNA duplex formation, a C-T mismatch is destabilizing in the center of a duplex, but behaves as a stable base pair at the terminus of a duplex. The C-T base pair is thought to contain two hydrogen bonds, but has thermodynamic parameters (delta Ho and delta Go of dissociation) that are similar to a G-C base pair. AMBER molecular mechanics calculations were performed to study the possible structural properties of DNA duplexes with central and terminal C-T combinations. These calculations also indicate that a central C-T pair destabilizes a duplex, while terminal C-T forms a stable base pair. Hydrogen bonding between cytosine and thymine occurs only in the energy-minimized structures when the helix diameter decreases and the propeller twist angle between the bases increases. These changes are found to occur only at the end of a duplex in the calculations, which may explain the experimental results.     

Parallel-stranded duplex DNA containing dA · dU base pairs

DNA oligonucleotides with dA and dU residues can form duplexes with trans d(A · U) base pairing and the sugar-phosphate backbone in a parallel-stranded orientation, as previously established for oligonucleotides with d(A · T) base pairs. The properties of such parallel-stranded DNA (ps-DNA) 25-mer duplexes have been characterized by absorption (uv), CD, ir, and fluorescence spectroscopy, as well as by nuclease sensitivity. Comparisons were made with duplex molecules containing (a) dT in both strands, (b) dU in one strand and dT in the second, and (c) the same base combinations in reference antiparallel-stranded (aps) structures. Thermodynamic analysis revealed that total replacement of deoxythymine by deoxyuridine was accompanied by destabilization of the ps-helix (reduction in T_m by −13°C in 2 mM MgGl₂, 10 mM Na-cacodylate). The U-containing ps-helix (U1 · U2) also melted 14°C lower than the corresponding aps-helix under the same ionic conditions; this difference was very close to that observed between ps and aps duplexes with d(A · T) base pairs. Force field minimized structures of the various ps and aps duplexes with either d(A · T) or d(A · U) base pairs ps/aps and dT/dU combinations are presented. The energy-minimized helical parameters did not differ significantly between the DNAs containing dT and dU. © 1996 John Wiley & Sons, Inc.     

Crystal structure of 2'-O-Me(CGCGCG)2, an RNA duplex at 1.30 A resolution. Hydration pattern of 2'-O-methylated RNA.

The molecular and crystal structure of 2'-O-Me (CGCGCG)2 has been determined using synchrotron radiation at near-atomic resolution (1.30 A), the highest resolution to date in the RNA field. The crystal structure is a half-turn A-type helix with some helical parameters deviating from canonical A-RNA, such as low base pair rise, elevated helical twist and inclination angles. In CG steps, inter-strand guanines are parallel while cytosines are not parallel. In steps GC this motif is reversed. This type of regularity is not seen in other RNA crystal structures. The structure includes 44 water molecules and two hydrated Mg2+ions one of which lies exactly on the crystallographic 2-fold axis. There are distinct patterns of hydration in the major and the minor grooves. The major groove is stabilised by water clusters consisting of fused five- and six-membered rings. Minor groove contains only a single row of water molecules; each water bridges either two self-parallel cytosines or two self-parallel guanines by a pair of hydrogen bonds. The structure provides the first view of the hydration scheme of 2'-O-methylated RNA duplex.     

Visualization of drug-nucleic acid interactions at atomic resolution. VII. Structure of an ethidium/dinucleoside monophosphate crystalline complex, ethidium: uridylyl(3'-5') adenosine

Ethidium forms a crystalline complex with the dinucleoside monophosphate, uridylyl (3'-5') adenosine (UpA). The complex crystallizes in the monoclinic space group P2l with unit cell dimensions, a = 13.704 A, b = 31.674 A, c = 15.131 A, beta = 113.9 degrees. This light atom structure has been solved to atomic resolution and refined by full matrix least squares to a residual of 0.12, using 3,034 observed reflections. The asymmetric unit consists of two ethidium molecules, two UpA molecules and 19 solvent molecules, a total of 145 non-hydrogen atoms. The two UpA molecules are hydrogen-bonded together by Watson-Crick type base pairing. Base-pairs in this duplex are separated by 6.7 A; this reflects intercalative binding by one of the ethidium molecules. The other ethidium molecule stacks on either side of the intercalated base-paired dinucleoside monophosphate, being related by a unit cell translation along the a axis. The conformation of the sugar-phosphate backbone accompanying intercalation has been accurately determined in this analysis, and contains the mixed sugar-puckering pattern: C3' endo (3'-5') C2' endo. This same structural feature has been observed in the ethidium-iodoUpA and ethidium-iodoCpG complexes, and exists in two additional structures containing ethidium-CpG. Taken together, these studies confirm our earlier sugar-puckering assignments and demonstrate that iodine covalently bound to the C5 position on uridine or cytosine does not alter the basic sugar-phosphate geometry or the mode of ethidium intercalation in these model studies. We have proposed this stereochemistry to explain the intercalation of ethidium (as well as other simple intercalators) into both DNA and into double-helical RNA, and discuss this aspect of our work further in this paper and in the accompanying papers.     

Sequence dependence for the energetics of dangling ends and terminal base pairs in ribonucleic acid

Stability increments of terminal unpaired nucleotides (dangling ends) and terminal base pairs on the core helixes AUGCAU and UGCGCA are reported. Enthalpy, entropy, and free energy changes of helix formation were measured spectrophotometrically for 18 oligoribonucleotides containing the core sequences. The results indicate 3' dangling purines add more stability than 3' dangling pyrimidines. In most cases, the additional stability from a 3' dangling end on an AU base pair is less than that on a GC base pair [Freier, S.M., Burger, B.J., Alkema, D., Neilson, T., & Turner, D.H. (1985) Biochemistry 22, 6198-6206]. The sequence dependence provides a test for the importance of dangling ends for various RNA interactions. Correlations are suggested with codon context effects and with the three-dimensional structure of yeast phenylalanine transfer RNA. In the latter case, all terminal unpaired nucleotides having stability increments more favorable than -1 kcal/mol are stacked on the adjacent base pair. All terminal unpaired nucleotides having stability increments less favorable than -0.3 kcal/mol are not stacked on the adjacent base pair. In several cases, this lack of stacking is associated with a turn in the sugar-phosphate backbone. This suggests stability increments measured on oligoribonucleotides may be useful for predicting tertiary structure in large RNA molecules. Comparison of the stability increments for terminal dangling ends and base pairs, and of terminal GC and AU base pairs, indicates the free energy increment associated with forming a hydrogen bond can be about -1 kcal/mol of hydrogen bond.     

An alternating sheared AA pair and elements of stability for a single sheared purine-purine pair flanked by sheared GA pairs in RNA

A previous NMR structure of the duplex 5'GGU GGA GGCU/PCCG AAG CCG5' revealed an unusually stable RNA internal loop with three consecutive sheared GA pairs. Here, we report NMR studies of two duplexes, 5'GGU GGA GGCU/PCCA AAG CCG5' (replacing the UG pair with a UA closing pair) and 5'GGU GAA GGCU/PCCG AAG CCG5' (replacing the middle GA pair with an AA pair). An unusually stable loop with three consecutive sheared GA pairs forms in the duplex 5'GGU GGA GGCU/PCCA AAG CCG5'. The structure contrasts with that reported for this loop in the crystal structure of the large ribosomal subunit of Deinococcus radiodurans [Harms, J., Schluenzen, F., Zarivach, R., Bashan, A., Gat, S., Agmon, I., Bartels, H., Franceschi, F., and Yonath, A. (2001) Cell 107, 679-688]. The middle AA pair in the duplex 5'GGU GAA GGCU/PCCG AAG CCG5' rapidly exchanges orientations, resulting in alternative base stacking and pseudosymmetry with exclusively sheared pairs. The U GAA G/G AAG C internal loop is 2.1 kcal/mol less stable than the U GGA G/G AAG C internal loop at 37 degrees C. Structural, energetic, and dynamic consequences upon functional group substitutions within related 3 x 3 and 3 x 6 internal loops are also reported.     

Inferring the conformation of RNA base pairs and triples from patterns of sequence variation.

The success of comparative analysis in resolving RNA secondary structure and numerous tertiary interactions relies on the presence of base covariations. Although the majority of base covariations in aligned sequences is associated to Watson-Crick base pairs, many involve non-canonical or restricted base pair exchanges (e.g. only G:C/A:U), reflecting more specific structural constraints. We have developed a computer program that determines potential base pairing conformations for a given set of paired nucleotides in a sequence alignment. This program (ISOPAIR) assumes that the base pair conformation is maintained through sequence variation without significantly affecting the path of the sugar-phosphate backbone. ISOPAIR identifies such 'isomorphic' structures for any set of input base pair or base triple sequences. The program was applied to base pairs and triples with known structures and sequence exchanges. In several instances, isomorphic structures were correctly identified with ISOPAIR. Thus, ISOPAIR is useful when assessing non-canonical base pair conformations in comparative analysis. ISOPAIR applications are limited to those cases where unusual base pair exchanges indeed reflect a non-canonical conformation.     

B-DNA structure is intrinsically polymorphic: even at the level of base pair positions

Increasingly exact measurement of single crystal X-ray diffraction data offers detailed characterization of DNA conformation, hydration and electrostatics. However, instead of providing a more clear and unambiguous image of DNA, highly accurate diffraction data reveal polymorphism of the DNA atomic positions and conformation and hydration. Here we describe an accurate X-ray structure of B-DNA, painstakingly fit to a multistate model that contains multiple competing positions of most of the backbone and of entire base pairs. Two of ten base-pairs of CCAGGCCTGG are in multiple states distinguished primarily by differences in slide. Similarly, all the surrounding ions are seen to fractionally occupy discrete competing and overlapping sites. And finally, the vast majority of water molecules show strong evidence of multiple competing sites. Conventional resolution appears to give a false sense of homogeneity in conformation and interactions of DNA. In addition, conventional resolution yields an average structure that is not accurate, in that it is different from any of the multiple discrete structures observed at high resolution. Because base pair positional heterogeneity has not always been incorporated into model-building, even some high and ultrahigh-resolution structures of DNA do not indicate the full extent of conformational polymorphism.     

G.T wobble base-pairing in Z-DNA at 1.0 A atomic resolution: the crystal structure of d(CGCGTG).

The DNA oligomer d(CGCGTG) crystallizes as a Z-DNA double helix containing two guanine-thymine base pair mismatches of the wobble type. The crystal diffracts to 1 A resolution and the structure has been solved and refined. At this resolution, a large amount of information is revealed about the organization of the water molecules in the lattice generally and more specifically around the wobble base pairs. By comparing this structure with the analogous high resolution structure of d(CGCGCG) we can visualize the structural changes as well as the reorganization of the solvent molecules associated with wobble base pairing. There is only a small distortion of the Z-DNA backbone resulting from introduction of the GT mismatched base pairs. The water molecules cluster around the wobble base pair taking up all of the hydrogen bonding capabilities of the bases due to wobble pairing. These bridging water molecules serve to stabilize the base-base interaction and, thus, may be generally important for base mispairing either in DNA or in RNA molecules.