Crystal structure of an alternating octamer r(GUAUGUA)dC with adjacent G x U wobble pairs

Structural Insights Into the 5′UG/3′GU Wobble Tandem in Complex With Ba2+ Cation

G•U wobble base pair frequently occurs in RNA structures. The unique chemical, thermodynamic, and structural properties of the G•U pair are widely exploited in RNA biology. In several RNA molecules, the G•U pair plays key roles in folding, ribozyme catalysis, and interactions with proteins. G•U may occur as a single pair or in tandem motifs with different geometries, electrostatics, and thermodynamics, further extending its biological functions. The metal binding affinity, which is essential for RNA folding, catalysis, and other interactions, differs with respect to the tandem motif type due to the different electrostatic potentials of the major grooves. In this work, we present the crystal structure of an RNA 8-mer duplex r[UCGUGCGA]₂, providing detailed structural insights into the tandem motif I (5′UG/3′GU) complexed with Ba²⁺ cation. We compare the electrostatic potential of the presented motif I major groove with previously published structures of tandem motifs I, II (5′GU/3′UG), and III (5′GG/3′UU). A local patch of a strongly negative electrostatic potential in the major groove of the presented structure forms the metal binding site with the contributions of three oxygen atoms from the tandem. These results give us a better understanding of the G•U tandem motif I as a divalent metal binder, a feature essential for RNA functions.

Modeling the vibrational couplings of nucleobases

Vibrational spectroscopy, in particular infrared spectroscopy, has been widely used to probe the three-dimensional structures and conformational dynamics of nucleic acids. As commonly used chromophores, the C=O and C=C stretch modes in the nucleobases exhibit distinct spectral features for different base pairing and stacking configurations. To elucidate the origin of their structural sensitivity, in this work, we develop transition charge coupling (TCC) models that allow one to efficiently calculate the interactions or couplings between the C=O and C=C chromophores based on the geometric arrangements of the nucleobases. To evaluate their performances, we apply the TCC models to DNA and RNA oligonucleotides with a variety of secondary and tertiary structures and demonstrate that the predicted couplings are in quantitative agreement with the reference values. We further elucidate how the interactions between the paired and stacked bases give rise to characteristic IR absorption peaks and show that the TCC models provide more reliable predictions of the coupling constants as compared to the transition dipole coupling scheme. The TCC models, together with our recently developed through-bond coupling constants and vibrational frequency maps, provide an effective theoretical strategy to model the vibrational Hamiltonian, and hence the vibrational spectra of nucleic acids in the base carbonyl stretch region directly from atomistic molecular simulations.

A novel form of RNA double helix based on G·U and C·A+ wobble base pairing

Wobble base pairs are critical in various physiological functions and have been linked to local structural perturbations in double-helical structures of nucleic acids. We report a 1.38-Å resolution crystal structure of an antiparallel octadecamer RNA double helix in overall A conformation, which includes a unique, central stretch of six consecutive wobble base pairs (W helix) with two G·U and four rare C·A+ wobble pairs. Four adenines within the W helix are N1-protonated and wobble-base-paired with the opposing cytosine through two regular hydrogen bonds. Combined with the two G·U pairs, the C·A+ base pairs facilitate formation of a half turn of W-helical RNA flanked by six regular Watson-Crick base pairs in standard A conformation on either side. RNA melting experiments monitored by differential scanning calorimetry, UV and circular dichroism spectroscopy demonstrate that the RNA octadecamer undergoes a pH-induced structural transition which is consistent with the presence of a duplex with C·A+ base pairs at acidic pH. Our crystal structure provides a first glimpse of an RNA double helix based entirely on wobble base pairs with possible applications in RNA or DNA nanotechnology and pH biosensors.

An innate twist between Crick's wobble and Watson-Crick base pairs

Non-Watson-Crick pairs like the G·U wobble are frequent in RNA duplexes. Their geometric dissimilarity (nonisostericity) with the Watson-Crick base pairs and among themselves imparts structural variations decisive for biological functions. Through a novel circular representation of base pairs, a simple and general metric scheme for quantification of base-pair nonisostericity, in terms of residual twist and radial difference that can also envisage its mechanistic effect, is proposed. The scheme is exemplified by G·U and U·G wobble pairs, and their predicable local effects on helical twist angle are validated by MD simulations. New insights into a possible rationale for contextual occurrence of G·U and other non-WC pairs, as well as the influence of a G·U pair on other non-Watson-Crick pair neighborhood and RNA-protein interactions are obtained from analysis of crystal structure data. A few instances of RNA-protein interactions along the major groove are documented in addition to the well-recognized interaction of the G·U pair along the minor groove. The nonisostericity-mediated influence of wobble pairs for facilitating helical packing through long-range interactions in ribosomal RNAs is also reviewed.

Testing the nearest neighbor model for canonical RNA base pairs: revision of GU parameters

Thermodynamic parameters for GU pairs are important for predicting the secondary structures of RNA and for finding genomic sequences that code for structured RNA. Optical melting curves were measured for 29 RNA duplexes with GU pairs to improve nearest neighbor parameters for predicting stabilities of helixes. The updated model eliminates a prior penalty assumed for terminal GU pairs. Six additional duplexes with the 5'GG/3'UU motif were added to the single representation in the previous database. This revises the ΔG°(37) for the 5'GG/3'UU motif from an unfavorable 0.5 kcal/mol to a favorable -0.2 kcal/mol. Similarly, the ΔG°(37) for the 5'UG/3'GU motif changes from 0.3 to -0.6 kcal/mol. The correlation coefficients between predicted and experimental ΔG°(37), ΔH°, and ΔS° for the expanded database are 0.95, 0.89, and 0.87, respectively. The results should improve predictions of RNA secondary structure.

Crystal structure analysis reveals functional flexibility in the selenocysteine-specific tRNA from mouse

Background

Selenocysteine tRNAs (tRNA^Sec) exhibit a number of unique identity elements that are recognized specifically by proteins of the selenocysteine biosynthetic pathways and decoding machineries. Presently, these identity elements and the mechanisms by which they are interpreted by tRNA^Sec-interacting factors are incompletely understood.

Methodology/Principal Findings

We applied rational mutagenesis to obtain well diffracting crystals of murine tRNA^Sec. tRNA^Sec lacking the single-stranded 3′-acceptor end (^ΔGCCARNA^Sec) yielded a crystal structure at 2.0 Å resolution. The global structure of ^ΔGCCARNA^Sec resembles the structure of human tRNA^Sec determined at 3.1 Å resolution. Structural comparisons revealed flexible regions in tRNA^Sec used for induced fit binding to selenophosphate synthetase. Water molecules located in the present structure were involved in the stabilization of two alternative conformations of the anticodon stem-loop. Modeling of a 2′-O-methylated ribose at position U34 of the anticodon loop as found in a sub-population of tRNA^Secin vivo showed how this modification favors an anticodon loop conformation that is functional during decoding on the ribosome. Soaking of crystals in Mn²⁺-containing buffer revealed eight potential divalent metal ion binding sites but the located metal ions did not significantly stabilize specific structural features of tRNA^Sec.

Conclusions/Significance

We provide the most highly resolved structure of a tRNA^Sec molecule to date and assessed the influence of water molecules and metal ions on the molecule's conformation and dynamics. Our results suggest how conformational changes of tRNA^Sec support its interaction with proteins.

Synthesis and hybridization of 2'-O-methyl-RNAs incorporating 2'-O-carbamoyluridine and unique participation of the carbamoyl group in U-G base pair

2'-O-Carbamoyluridine (U(cm)) was synthesized and incorporated into DNAs and 2'-O-Me-RNAs. The oligonucleotides incorporating U(cm) formed less stable duplexes with their complementary and U(cm)-U, U(cm)-C single-base mismatched DNAs and RNAs in comparison with those without the carbamoyl group. On the contrary, the T(m) analyses revealed that the duplexes with a mismatched U(cm)-G base pair showed almost the same thermostability as the corresponding unmodified duplexes. Molecular dynamics (MD) simulations of the U(cm)-modified 2'-O-Me-RNA/RNA duplexes with U(cm)-G mismatched base pair suggested that the carbamoyl group could participate in the U(cm)-G base pair by an additional intermolecular hydrogen bond between the carbamoyl oxygen and the H2 of the guanine base.

The electrostatic characteristics of G.U wobble base pairs

G.U wobble base pairs are the most common and highly conserved non-Watson-Crick base pairs in RNA. Previous surface maps imply uniformly negative electrostatic potential at the major groove of G.U wobble base pairs embedded in RNA helices, suitable for entrapment of cationic ligands. In this work, we have used a Poisson-Boltzmann approach to gain a more detailed and accurate characterization of the electrostatic profile. We found that the major groove edge of an isolated G.U wobble displays distinctly enhanced negativity compared with standard GC or AU base pairs; however, in the context of different helical motifs, the electrostatic pattern varies. G.U wobbles with distinct widening have similar major groove electrostatic potentials to their canonical counterparts, whereas those with minimal widening exhibit significantly enhanced electronegativity, ranging from 0.8 to 2.5 kT/e, depending upon structural features. We propose that the negativity at the major groove of G.U wobble base pairs is determined by the combined effect of the base atoms and the sugar-phosphate backbone, which is impacted by stacking pattern and groove width as a result of base sequence. These findings are significant in that they provide predictive power with respect to which G.U sites in RNA are most likely to bind cationic ligands.

Expanding functionality of RNA: synthesis and properties of RNA containing imidazole modified tandem G-U wobble base pairs

Imidazole modification at C-5 of uridine that is part of tandem G-U wobble base pairs causes slight reduction of thermal stability (DeltaDeltaG(0)(310) < 0.4 kcal mol(-1)) and relatively small change in hydration of short RNA helices.

Conformational determinants of tandem GU mismatches in RNA: insights from molecular dynamics simulations and quantum mechanical calculations

Structure and energetic properties of base pair mismatches in duplex RNA have been the focus of numerous investigations due to their role in many important biological functions. Such efforts have contributed to the development of models for secondary structure prediction of RNA, including the nearest-neighbor model. In RNA duplexes containing GU mismatches, 5'-GU-3' tandem mismatches have a different thermodynamic stability than 5'-UG-3' mismatches. In addition, 5'-GU-3' mismatches in some sequence contexts do not follow the nearest-neighbor model for stability. To characterize the underlying atomic forces that determine the structural and thermodynamic properties of GU tandem mismatches, molecular dynamics (MD) simulations were performed on a series of 5'-GU-3' and 5'-UG-3' duplexes in different sequence contexts. Overall, the MD-derived structural models agree well with experimental data, including local deviations in base step helicoidal parameters in the region of the GU mismatches and the model where duplex stability is associated with the pattern of GU hydrogen bonding. Further analysis of the simulations, validated by data from quantum mechanical calculations, suggests that the experimentally observed differences in thermodynamic stability are dominated by GG interstrand followed by GU intrastrand base stacking interactions that dictate the one versus two hydrogen bonding scenarios for the GU pairs. In addition, the inability of 5'-GU-3' mismatches in different sequence contexts to all fit into the nearest-neighbor model is indicated to be associated with interactions of the central four base pairs with the surrounding base pairs. The results emphasize the role of GG and GU stacking interactions on the structure and thermodynamics of GU mismatches in RNA.

Structures of non-canonical tandem base pairs in RNA helices: review

The structures of tandem non-canonical base pairs, a frequently recurring motif in RNA molecules, are reviewed and analysed. The tandem non-canonical base pair motifs can be roughly divided in three groups, containing seven subgroups based on their base pairing patterns and local geometries. Structural details and helical parameters that can be used to numerically distinguish between the subgroups are tabulated. Remarkably, while the individual helical twists of the tandem and adjacent base pair steps can be substantially smaller or larger than the typical A-form value of 32.7 degrees, the average value is close to A-form. This and other striking regularities resulting from compensating geometrical adjustments, important for understanding and predicting the configurations of non-canonical base pairs geometries are discussed.

Site-specific modification of functional groups in genomic hepatitis delta virus (HDV) ribozyme

Human hepatitis delta (HDV) ribozyme is one of small ribozymes, such as hammerhead and hairpin ribozymes, etc. Its secondary structure shows pseudoknot structure composed of four stems (I to IV) and three single-stranded regions (SSrA, -B and -C). The 3D structure of 3'-cleaved product of genomic HDV ribozyme provided extensive information about tertiary hydrogen bonding interactions between nucleotide bases, phosphate oxygens and 2'OHs including new stem structure P1.1. To analyze the role of these hydrogen bond networks in the catalytic reaction, site-specific atomic-level modifications (such as deoxynucleotides, deoxyribosyl-2-aminopurine, deoxyribosylpurine, 7-deaza-ribonucleotide and inosine) were incorporated in the smallest trans-acting HDV ribozyme (47-mer). Kinetic analysis of these ribozyme variants demonstrated the importance of the two W-C base pairs of P1.1 for cleavage; in addition, the results suggest that all hydrogen bond interactions detected in the crystal structure involving 2'-OH and N7 atoms are present in the active ribozyme structure. In most of the variants, the relative reduction in kobs caused by substitution of the 2'-OH group correlated with the number of hydrogen bonds affected by the substitution. However G74 and C75 may have more than one hydrogen bond involving the 2'-OH in both the trans- and cis-acting HDV ribozyme. Moreover, in variants in which N7 was deleted, kobs was reduced 5- to 15-fold, it may suggest that N7 assists in coordinating Mg2+ ions or water molecules which bind with weak affinity in the active structure.

Hydrogen and hydration of DNA and RNA oligonucleotides

In this paper, hydrogen bonding interaction and hydration in crystal structures of both DNA and RNA oligonucleotides are discussed. Their roles in the formation and stabilization of oligonucleotides have been covered. Details of the Watson-Crick base pairs G.C and A.U in DNA and RNA are illustrated. The geometry of the wobble (mismatched) G.U base pairs and the cis and almost trans conformations of the mismatched U.U base pairs in RNA is described. The difference in hydration of the Watson-Crick base pairs G.C, A.U and the wobble G.U in different sequences of codon-anticodon interaction in double helical molecules are indicative of the effect of hydration. The hydration patterns of the phosphate, the 2'-hydroxyl groups, the water bridges linking the phosphate group, N7 (purine) and N4 of Cs or O4 of Us in the major groove, the water bridges between the 2'-hydroxyl group and N3 (purine) and O2 (pyrimidine) in the minor groove are discussed.

Thermodynamic stabilities of internal loops with GU closing pairs in RNA

Many internal loops that form tertiary contacts in natural RNAs have GU closing pairs; examples include the tetraloop receptor and P1 helix docking site in group I introns. Thus, thermodynamic parameters of internal loops with GU closing pairs can contribute to the prediction of both secondary and tertiary structure. Oligoribonucleotide duplexes containing small internal loops with GU closing pairs were studied by optical melting, one-dimensional imino proton NMR, and one-dimensional phosphorus NMR. The thermodynamic stabilities of asymmetric internal loops with GU closing pairs relative to those of loops with GC closing pairs may be explained by hydrogen bonds. In contrast, the free energy increments for symmetric internal loops of two noncanonical pairs with GU closing pairs relative to loops with GC closing pairs show much more sequence dependence. Imino proton and phosphorus NMR spectra suggest that some GA pairs adjacent to GU closing pairs may form an overall thermodynamically stable but non-A-form conformation.

B-form to A-form conversion by a 3'-terminal ribose: crystal structure of the chimera d(CCACTAGTG)r(G)

The crystal structure of the chimerical decamer d(CCACTAGTG)r(G), bearing a 3'-terminal ribo-guanidine, has been solved and refined at 1.8 A resolution (R-factor 16.6%; free R-factor 22.8%). The decamer crystallizes in the orthorhombic space group P2(1)2(1)2(1) with unit cell constants a = 23.90 A, b = 45.76 A and c = 49.27 A. The structure was solved by molecular replacement using the coordinates of the isomorphous chimera r(GCG)d(TATACGC). The final model contains one duplex and 77 water molecules per asymmetric unit. Surprisingly, all residues adopt a conformation typical for A-form nucleic acids (C3'-endo type sugar pucker) although the all-DNA analog, d(CCACTAGTGG), has been crystallized in the B-form. Comparing circular dichroism spectra of the chimera and the corresponding all-DNA sequence reveals a similar trend of the former molecule to adopt an A-like conformation in solution. The results suggest that the preference of ribonucleotides for the A-form is communicated into the 5'-direction of an oligonucleotide strand, although direct interactions of the 2'-hydroxyl group can only be discerned with nucleotides in the 3'-direction of a C3'-endo puckered ribose. These observations imply that forces like water-mediated contacts, the concerted motions of backbone torsion angles, and stacking preferences, are responsible for such long-range influences. This bi-directional structural communication originating from a ribonucleotide can be expected to contribute to the stability of the A-form within all-RNA duplexes.

Synthesis and crystal structure of an octamer RNA r(guguuuac)/r(guaggcac) with G.G/U.U tandem wobble base pairs: comparison with other tandem G.U pairs

We have determined the crystal structure of the RNA octamer duplex r(guguuuac)/r(guaggcac) with a tandem wobble pair, G.G/U.U (motif III), to compare it with U.G/G.U (motif I) and G.U/U.G (motif II) and to better understand their relative stabilities. The crystal belongs to the rhombohedral space group R3. The hexagonal unit cell dimensions are a = b = 41.92 A, c = 56.41 A, and gamma = 120 degrees, with one duplex in the asymmetric unit. The structure was solved by the molecular replacement method at 1.9 A resolution and refined to a final R: factor of 19.9% and R(free) of 23.3% for 2862 reflections in the resolution range 10.0-1.9 A with F >/= 2sigma(F). The final model contains 335 atoms for the RNA duplex and 30 water molecules. The A-RNA stacks in the familiar head-to-tail fashion forming a pseudo-continuous helix. The uridine bases of the tandem U.G pairs have slipped towards the minor groove relative to the guanine bases and the uridine O2 atoms form bifurcated hydrogen bonds with the N1 and N2 of guanines. The N2 of guanine and O2 of uridine do not bridge the 'locked' water molecule in the minor groove, as in motifs I and II, but are bridged by water molecules in the major groove. A comparison of base stacking stabilities of motif III with motifs I and II confirms the result of thermodynamic studies, motif I > motif III > motif II.

The crystal structure of the octamer [r(guauaca)dC]2 with six Watson-Crick base-pairs and two 3' overhang residues

The crystal structure of an alternating RNA octamer, r(guauaca)dC (RNA bases are in lower case while the only DNA base is in upper case), with two 3' overhang residues one of them a terminal deoxycytosine and the other a ribose adenine, has been determined at 2.2 A resolution. The refined structure has an Rwork 18.6% and Rfree 26.8%. There are two independent duplexes (molecules I and II) in the asymmetric unit cell, a = 24.95, b = 45.25 and c = 73.67 A, with space group P2(1)2(1)2(1). Instead of forming a blunt end duplex with two a+.c mispairs and six Watson-Crick base-pairs, the strands in the duplex slide towards the 3' direction forming a two-base overhang (radC) and a six Watson-Crick base-paired duplex. The duplexes are bent (molecule I, 20 degrees; molecule II, 25 degrees) and stack head-to-head to form a right-handed superhelix. The overhang residues are looped out and the penultimate adenines of the two residues at the top end (A15) are anti and at the bottom (A7) end are syn. The syn adenine bases form minor groove A*(G.C) base triples with C8-H...N2 hydrogen bonds. The anti adenine in molecule II also forms a triple and a different C2-H...N3 hydrogen bond, while the other anti adenine in molecule I does not, it stacks on the looped out overhang base dC. The 3' terminal deoxycytosines form two stacked hemiprotonated trans d(C.C)+ base-pairs and the pseudo dyad related molecules form four consecutive deoxyribose and ribose zipper hydrogen bonds in the minor groove.

Structural Studies of Model RNA Helices with Relevance to Aminoacyl-tRNA Synthetase Specificity and HIV Reverse Transcription

Abstract We describe high-resolution crystal structures of synthetic nucleic-acid fragments determined as part of an effort to understand determinants of sequence-specific protein binding on the level of double-helix structure. In a first set of experiments, 7-base-pair RNA duplexes representing the acceptor-stem helix of Escherichia coli tRNA(Ala) and variants thereof were characterized at atomic resolution. The structures revealed a standard A-form double helix locally perturbed by a G·U wobble base pair at sequence position 3/70 of the tRNA. The G·U pair shows a characteristic hydration pattern which must be considered an integral part of the double-helix structure. It does not seem to exert a global effect on the duplex structure. A second experiment concerned the chimeric DNA-RNA hybrid structure formed transiently during initiation of minus-strand synthesis by the reverse transcriptase of HIV-1. The crystal structure of an 8-base-pair duplex with an RNA template strand derived from HIV-1 and a complementary strand representing the junction between the tRNA(Lys,3) RNA primer and the newly synthesized DNA strand was solved at a resolution of 1.9 Å. As before, the double helix was found to adopt standard A-type conformation with only local variations of backbone conformation. Based on the global helix structure as present in the crystal, it remains difficult to explain the preference of the reverse-transcriptase-associated RNAse H activity for certain sites of the template strand. Structural plasticity near the main cleavage site in suggested to govern cutting preferences. In both systems investigated, structural studies by NMR spectroscopy were carried out by others in parallel. In both cases, the solution structures are in partial disagreement with the crystallographic results by describing a significantly higher level of deviation from the canonical A-conformation.

Correlation of deformability at a tRNA recognition site and aminoacylation specificity

The fidelity of protein synthesis depends on specific tRNA aminoacylation by aminoacyl-tRNA synthetase enzymes, which in turn depends on the recognition of the identity of particular nucleotides and structural features in the substrate tRNA. These features generally reside within the acceptor helix, the anticodon stem-loop, and in some systems the variable pocket of the tRNA. In the alanine system, fidelity is ensured by a G.U wobble base pair located at the third position within the acceptor helix of alanine tRNA. We have investigated the activity of mutant alanine tRNAs to explore the mechanism of enzyme recognition. Here we show that the mismatched pair C-C is an excellent substitute for G.U in alanine-tRNA-knockout cells. A structural investigation by NMR spectroscopy of the C-C RNA acceptor end reveals that the two cytosines are intercalated into the helix, and that C-C exists in multiple conformations. Structural heterogeneity also is present in the wild-type G.U RNA, whereas inactive Watson-Crick helices are structurally rigid. The correlation between functional and structural data suggests that the G.U pair provides a distinctive structure and a point of deformability that allow the tRNA acceptor end to fit into the active site of the alanyl-tRNA synthetase. Fidelity is ensured because noncognate and inactive mutant tRNAs are bound in the active site in an incorrect conformation that reduces enzymatic activity.

Structure and thermodynamics of metal binding in the P5 helix of a group I intron ribozyme

The solution structure of an RNA hairpin modelling the P5 helix of a group I intron, complexed with Co(NH3)63+, has been determined by nuclear magnetic resonance. Co(NH3)63+, which possesses a geometry very close to Mg(H2O)62+, was used to identify and characterize a Mg2+binding site in the RNA. Strong and positive intermolecular nuclear Overhauser effect (NOE) cross-peaks define a specific complex in which the Co(NH3)63+molecule is in the major groove of tandem G.U base-pairs. The structure of the RNA is characterized by a very low twist angle between the two G.U base-pairs, providing a flat and narrowed major groove. The Co(NH3)63+, although highly localized, is free to rotate to hydrogen bond in several ways to the O4 atoms of the uracil bases and to N7 and O6 of the guanine bases. Negative and small NOE cross-peaks to other protons in the sequence reveal a non-specific or delocalized interaction, characterized by a high mobility of the cobalt ion. Mn2+titrations of P5 show specific broadening of protons of the G.U base-pairs that form the metal ion binding site, in agreement with the NOE data from Co(NH3)63+. Binding constants for the interaction of Co(NH3)63+and of Mg2+to P5 were determined by monitoring imino proton chemical shifts during titration of the RNA with the metal ions. Dissociation constants are on the order of 0.1 mM for Co(NH3)63+and 1 mM for Mg2+. Binding studies were done on mutants with sequences corresponding to the three orientations of tandem G.U base-pairs. The affinities of Co(NH3)63+and Mg2+for the tandem G.U base-pairs depend strongly on their sequences; the differences can be understood in terms of the different structures of the corresponding metal ion-RNA complexes. Substitution of G.C or A.U for G.U pairs also affected the binding, as expected. These structural and thermodynamic results provide systematic new information about major groove metal ion binding in RNA.

Crystal structure of an RNA duplex r(G GCGC CC)2 with non-adjacent G*U base pairs

The crystal structure of a self-complementary RNA duplex r(GGGCGCUCC)2with non-adjacent G*U and U*G wobble pairs separated by four Watson-Crick base pairs has been determined to 2.5 A resolution. Crystals belong to the space group R3; a = 33.09 A,alpha = 87.30 degrees with a pseudodyad related duplex in the asymmetric unit. The structure was refined to a final Rworkof 17.5% and Rfreeof 24.0%. The duplexes stack head-to-tail forming infinite columns with virtually no twist at the junction steps. The 3'-terminal cytosine nucleosides are disordered and there are no electron densities, but the 3' penultimate phosphates are observed. As expected, the wobble pairs are displaced with guanine towards the minor groove and uracil towards the major groove. The largest twist angles (37.70 and 40.57 degrees ) are at steps G1*C17/G2*U16 and U7*G11/C8*G10, while the smallest twist angles (28.24 and 27.27 degrees ) are at G2*U16/G3*C15 and C6*G12/U7*G11 and conform to the pseudo-dyad symmetry of the duplex. The molecule has two unequal kinks (17 and 11 degrees ) at the wobble sites and a third kink at the central G5 site which may be attributed to trans alpha (O5'-P), trans gamma (C4'-C5') backbone conformations. The 2'-hydroxyl groups in the minor groove form inter-column hydrogen bonding, either directly or through water molecules.

Detailed analysis of stem I and its 5' and 3' neighbor regions in the trans-acting HDV ribozyme

To determine the stem I structure of the human hepatitis delta virus (HDV) ribozyme, which is related to the substrate sequence in the trans -acting system, we kinetically studied stem I length and sequences. Stem I extension from 7 to 8 or 9 bp caused a loss of activity and a low amount of active complex with 9 bp in the trans -acting system. In a previous report, we presented cleavage in a 6 bp stem I. The observed reaction rates indicate that the original 7 bp stem I is in the most favorable location for catalytic reaction among the possible 6-8 bp stems. To test base specificity, we replaced the original GC-rich sequence in stem I with AU-rich sequences containing six AU or UA base pairs with the natural +1G.U wobble base pair at the cleavage site. The cis -acting AU-rich molecules demonstrated similar catalytic activity to that of the wild-type. In trans -acting molecules, due to stem I instability, reaction efficiency strongly depended on the concentration of the ribozyme-substrate complex and reaction temperature. Multiple turnover was observed at 37 degreesC, strongly suggesting that stem I has no base specificity and more efficient activity can be expected under multiple turnover conditions by substituting several UA or AU base pairs into stem I. We also studied the substrate damaging sequences linked to both ends of stem I for its development in therapeutic applications and confirmed the functions of the unique structure.

Hydration of RNA base pairs

The hydration patterns around the RNA Watson-Crick and non-Watson-Crick base pairs in crystals are analyzed and described. The results indicate that (i) the base pair hydration is mostly "in-plane"; (ii) eight hydration sites surround the Watson-Crick G-C and A-U base pairs, with five in the deep and three in the shallow groove, an observation which extends the characteristic isostericity of Watson-Crick pairs; (iii) while the hydration around G-C base pairs is well defined, the hydration around A-U base pairs is more diffuse; (iv) the hydration sites close to the phosphate groups are the best defined and the most recurrent ones; (v) a string of water molecules links the two shallow groove 2'-hydroxyl groups, and (vi) the water molecules fit into notches, the size and accessibility of which are almost as important as the number and strength of the hydrophilic groups lining the cavity. Residence times of water molecules at specific hydration sites, inferred from molecular dynamics simulations, are discussed in the light of present data.

Structure of a 16-mer RNA duplex r(GCAGACUUAAAUCUGC)2 with wobble C.A+ mismatches

The crystal structure of a 16-mer, the longest known RNA duplex, has been determined at 2.5 A resolution. The hexadecamer r(GCAGACUUAAAUCUGC) contains isolated C.A/A.C mismatches with two hydrogen bonds. The two hydrogen bonds in the mismatches suggests that N1 of A is protonated even though the crystallization was done at neutral pH. Therefore, the C.A mismatch is a C.A+ wobble similar to the G.U wobble. The two C.A+ pairs are isolated by four Watson-Crick pairs and flanked by five Watson-Crick base-pairs on either sides. Kinks/bends of 20 degrees are observed at the wobble sites. The Watson-Crick base-pair A5.U26 on the 5'-side of the first C6.A27(+) wobble has a twist angle of 27 degrees compared to the 3'-side U7.A28 pair of 36 degrees. The twist angles are reversed (37 degrees and 26 degrees) in the second A11(+).C22 wobble because of the approximate dyad in the molecule, the flanking base-pair sequences are A.U pairs. The wobbles expand the major groove to 7.1 A/7.3 A. The duplexes form helical columns and are tightly packed around the 31-screw axis. The minor grooves of adjacent columns in juxtaposition interact through the O2' atoms and the anionic phosphate oxygen atoms.

Crystal structure of the spliceosomal U2B"-U2A' protein complex bound to a fragment of U2 small nuclear RNA

We have determined the crystal structure at 2.4 A resolution of a ternary complex between the spliceosomal U2B"/U2A' protein complex and hairpin-loop IV of U2 small nuclear RNA. Unlike its close homologue the U1A protein, U2B" binds to its cognate RNA only in the presence of U2A', which contains leucine-rich repeats in its sequence. The concave surface of a parallel beta-sheet within the leucine-rich-repeat region of U2A' interacts with the ribonucleoprotein domain of U2B" on the surface opposite its RNA-binding surface. The basic carboxy-terminal region of U2A' interacts with the RNA stem. The crystal structure reveals how protein-protein interaction regulates RNA-binding specificity, and how replacing only a few key residues allows the U2B" and U1A proteins to discriminate between their cognate RNA hairpins by forming alternative networks of interactions.

MANIP: an interactive tool for modelling RNA

Large RNA structures can be viewed as assemblies of smaller units or modules that are usually clearly identified (helices, hairpin loops, other recurrent motifs, etc.). We have developed a program, MANIP, which allows the rapid assembly of separate motifs (each with a specified sequence) into a complex three-dimensional architecture. The already determined modules are present in a database from which they can be extracted with the appropriate sequence. Their assembly is performed in real time on the computer screen with buttons and dials that command rotation and translation of any chosen fragment with respect to the chosen pivot, or that generate all possible variations of any torsion angle within a specified segment either in the 5' or in the 3' direction. The possible in-built manipulations follow the general stereochemical rules of RNA structure. MANIP automatically recognizes and displays the allowed and nonallowed hydrogen bonds between the residues. The program is interfaced with a rapid and automatic online refinement tool of partial or full assemblies, NUCLIN-NUCLSQ. The refinement protocol incorporates canonical as well as noncanonical base pairing constraints together with restraints imposed by covalent geometry, stereochemistry, and van der Waals contacts. The computer package runs on UNIX Silicon Graphics workstations and is written in C with OpenGL and X11/Motif libraries.

RNA structure

New information concerning RNA structure is accumulating at an ever increasing rate-from short helices with mismatched bases of 5S rRNA and complex RNA aptamers. The importance of recurring structural motifs, ion binding, and the kinetics and energetics of folding in RNA structure and function is now being recognized and addressed.

The importance of tRNA backbone-mediated interactions with synthetase for aminoacylation

We have identified six new aminoacylation determinants of Escherichia coli tRNAGln in a genetic and biochemical analysis of suppressor tRNA. The new determinants occupy the interior of the acceptor stem, the inside corner of the L shape, and the anticodon loop of the molecule. They supplement the primary determinants located in the anticodon and acceptor end of tRNAGln described previously. Remarkably, the three-dimensional structure of the complex between tRNAGln and glutaminyl-tRNA synthetase shows that the enzyme interacts with the phosphate-sugar backbone but not the base of every new determinant. Moreover, a small protein motif interacts with five of these determinants, and it binds proximal to the sixth. The motif also interacts with the middle base of the anticodon and with the backbones of six other nucleotides. Our results emphasize that synthetase recognition of tRNA is more elaborate than amino acid side chains of the enzyme interacting with nucleotide bases of the tRNA. Recognition also includes synthetase interaction with tRNA backbone functionalities whose distinctive locations in three-dimensional space are exquisitely determined by the tRNA sequence.

The loop E-loop D region of Escherichia coli 5S rRNA: the solution structure reveals an unusual loop that may be important for binding ribosomal proteins

BACKGROUND

5S ribosomal RNA is the smallest rRNA. Its Watson-Crick helices were identified more than 20 years ago, but the conformations of its loops have long defied analysis. One of the three arms of 5S rRNA, residues 69-106 in Escherichia coli, contains a 14-residue internal loop called loop E. The sequence of loop E is conserved within kingdoms, and is terminated by a pyrimidine-rich loop called loop D. Loop E is the binding site for the ribosomal protein L25 in the E. coli ribosome.

RESULTS

The solution structure of a 42-nucleotide derivative of E. coli 5S rRNA that includes loops D and E has been determined by nuclear magnetic resonance spectroscopy. Formally, loop E is not a loop at all; it is a double helical structure that contains seven, consecutive non-Watson-Crick base pairs. The major groove of the molecule is narrowed in loop E, and an unusual array of hydrogen-bond donors and acceptors appear in its minor groove. Loop D, which on paper looks like a three-pyrimidine terminal loop closed by a GC, is better thought of as a five-base loop because its closing GC is not a normal Watson-Crick pair. The two pyrimidines on the 5'-side of the loop are stacked on each other, and tilt into the minor groove of the adjacent helix. The third pyrimidine is fully exposed to solvent.

CONCLUSIONS

This structure rationalizes all the biochemical and chemical protection data available for the loop E-loop D arm of intact 5S rRNA. While the molecule is double helical over its entire length, the geometry of its internal loop is highly irregular, and its irregularities may explain why the loop E-loop D arm of 5S rRNA interacts specifically with ribosomal protein L25 in E. coli.

Metals, motifs, and recognition in the crystal structure of a 5S rRNA domain

Two new RNA structures portray how non-Watson-Crick base pairs and metal ions can produce a unique RNA shape suitable for recognition by proteins. The crystal structures of a 62 nt domain of E. coli 5S ribosomal RNA and a duplex dodecamer encompassing an internal loop E have been determined at 3.0 and 1.5 A, respectively. This loop E region is distorted by three "cross-strand purine stacks" and three novel, water-mediated noncanonical base pairs and stabilized by a four metal ion zipper. These features give its minor groove a unique hydrogen-bonding surface and make the adjacent major groove wide enough to permit recognition by the ribosomal protein L25, which is expected to bind to this surface.

Crystal structure of r(GUGUGUA)dC with tandem G x U/U x G wobble pairs with strand slippage

To better understand the frequent occurrence of adjacent wobble pairs in ribosomal RNAs we have determined the crystal structure of the RNA duplex, r(GUGUGUA)dC with the 3'-terminal deoxy C residue. Two different crystal forms of the duplex were obtained and both belong to the rhombohedral space group, R3. Crystal form I has hexagonal unit cell dimensions, a = b = 40.82 A and c = 66.09 A and diffracts to 1.58 A resolution, while crystal form II has a = b = 47.11 A and c = 59.86 A, diffracting only to 2.50 A resolution. Both structures were solved by the molecular replacement method using different starting models. In spite of the large differences in the cell dimensions the overall structures in both crystals are similar. Instead of the expected blunt-end duplex with four consecutive G x U pairs, the slippage of the strands resulted in two different tandem G x U/U x G wobble pairs involving two of the central and two of the 5' overhang bases, still yielding a total of four wobble pairs. These tandem wobble pairs are flanked by two Watson-Crick pairs. The A-type duplexes stack in the familiar head-to-tail fashion forming a pseudocontinuous helix. The wobble pairs of the present motif II (G x U/U x G) structure stack with a low twist angle of 25.3 degrees in contrast to that of motif I (U x G/G x U), 38.1 degrees. The four wobble pairs are characteristically heavily hydrated in both the grooves accounting for their stability.