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.     

Crystal structure of an RNA octamer duplex r(CCCIUGGG)2 incorporating tandem I.U wobbles.

The crystal structure of the RNA octamer duplex r(CCCIUGGG)2has been elucidated at 2.5 A resolution. The crystals belong to the space group P21and have unit cell constants a = 33.44 A, b = 43.41 A, c = 49.39 A and beta = 104.7 degrees with three independent duplexes (duplexes 1-3) in the asymmetric unit. The structure was solved by the molecular replacement method and refined to an Rwork/Rfree of 0.185/0.243 using 3765 reflections between 8.0 and 2.5 A. This is the first report of an RNA crystal structure incorporating I.U wobbles and three molecules in the asymmetric unit. Duplex 1 displays a kink of 24 degrees between the mismatch sites, while duplexes 2 and 3 have two kinks each of 19 degrees and 27 degrees, and 24 degrees and 29 degrees, respectively, on either side of the tandem mismatches. At the I.U/U.I mismatch steps, duplex 1 has a twist angle of 33.9 degrees, close to the average for all base pair steps, but duplexes 2 and 3 are underwound, with twist angles of 24.4 degrees and 26.5 degrees, respectively. The tandem I.U wobbles show intrastrand purine-pyrimidine stacking but exhibit interstrand purine-purine stacking with the flanking C.G pairs. The three independent duplexes are stacked non-coaxially in a head-to-tail fashion to form infinite pseudo-continuous helical columns which form intercolumn hydrogen bonding interactions through the 2'-hydroxyl groups where the minor grooves come together.     

Thermodynamic characterization of tandem mismatches found in naturally occurring RNA

Although all sequence symmetric tandem mismatches and some sequence asymmetric tandem mismatches have been thermodynamically characterized and a model has been proposed to predict the stability of previously unmeasured sequence asymmetric tandem mismatches [Christiansen,M.E. and Znosko,B.M. (2008) Biochemistry, 47, 4329–4336], experimental thermodynamic data for frequently occurring tandem mismatches is lacking. Since experimental data is preferred over a predictive model, the thermodynamic parameters for 25 frequently occurring tandem mismatches were determined. These new experimental values, on average, are 1.0 kcal/mol different from the values predicted for these mismatches using the previous model. The data for the sequence asymmetric tandem mismatches reported here were then combined with the data for 72 sequence asymmetric tandem mismatches that were published previously, and the parameters used to predict the thermodynamics of previously unmeasured sequence asymmetric tandem mismatches were updated. The average absolute difference between the measured values and the values predicted using these updated parameters is 0.5 kcal/mol. This updated model improves the prediction for tandem mismatches that were predicted rather poorly by the previous model. This new experimental data and updated predictive model allow for more accurate calculations of the free energy of RNA duplexes containing tandem mismatches, and, furthermore, should allow for improved prediction of secondary structure from sequence.     

Thermodynamic characterization of the complete set of sequence symmetric tandem mismatches in RNA and an improved model for predicting the free energy contribution of sequence asymmetric tandem mismatches

Because of the availability of an abundance of RNA sequence information, the ability to rapidly and accurately predict the secondary structure of RNA from sequence is becoming increasingly important. A common method for predicting RNA secondary structure from sequence is free energy minimization. Therefore, accurate free energy contributions for every RNA secondary structure motif are necessary for accurate secondary structure predictions. Tandem mismatches are prevalent in naturally occurring sequences and are biologically important. A common method for predicting the stability of a sequence asymmetric tandem mismatch relies on the stabilities of the two corresponding sequence symmetric tandem mismatches [Mathews, D. H., Sabina, J., Zuker, M., and Turner, D. H. (1999) J. Mol. Biol. 288, 911-940]. To improve the prediction of sequence asymmetric tandem mismatches, the experimental thermodynamic parameters for the 22 previously unmeasured sequence symmetric tandem mismatches are reported. These new data, however, do not improve prediction of the free energy contributions of sequence asymmetric tandem mismatches. Therefore, a new model, independent of sequence symmetric tandem mismatch free energies, is proposed. This model consists of two penalties to account for destabilizing tandem mismatches, two bonuses to account for stabilizing tandem mismatches, and two penalties to account for A-U and G-U adjacent base pairs. This model improves the prediction of asymmetric tandem mismatch free energy contributions and is likely to improve the prediction of RNA secondary structure from sequence.     

Thermodynamic stability of DNA tandem mismatches

The thermodynamics of nine hairpin DNAs were evaluated using UV-monitored melting curves and differential scanning calorimetry (DSC). Each DNA has the same five-base loop and a stem with 8-10 base pairs. Five of the DNAs have a tandem mismatch in the stem, while four have all base pairs. The tandem mismatches examined (ga/ga, aa/gc, ca/gc, ta/ac, and tc/tc) spanned the range of stability observed for this motif in a previous study of 28 tandem mismatches. UV-monitored melting curves were obtained in 1.0 M Na(+), 0.1 M Na(+), and 0.1 M Na(+) with 5 mM Mg(2+). DSC studies were conducted in 0.1 M Na(+). Transition T(m) values were unchanged over a 50-fold range of strand concentration. Model-independent enthalpy changes (DeltaH degrees ) evaluated by DSC were in good agreement (+/-8%) with enthalpy values determined by van't Hoff analyses of the melting curves in 0.1 M Na(+). The average heat capacity change (DeltaC(p)) associated with the hairpin to single strands transitions was estimated from plots of DeltaH degrees and DeltaS degrees with T(m) and ln T(m), respectively, and from profiles of DSC curves. The average DeltaC(p) values (113 +/- 9 and 42 +/- 27 cal x K(-1) x mol(-1) of bp), were in the range of values reported in previous studies. Consideration of DeltaC(p) produced large changes in DeltaH degrees and DeltaS degrees extrapolated from the transition region to 37 degrees C and smaller but significant changes to free energies. The loop free energy of the five tandem mismatches at 37 degrees C varied over a range of approximately 4 kcal x mol(-1) for each solvent.     

Effects of N2,N2-dimethylguanosine on RNA structure and stability: crystal structure of an RNA duplex with tandem m2 2G:A pairs

Methylation of the exocyclic amino group of guanine is a relatively common modification in rRNA and tRNA. Single methylation (N²-methylguanosine, m²G) is the second most frequently encountered nucleoside analog in Escherichia coli rRNAs. The most prominent case of dual methylation (N²,N²-dimethylguanosine, m²₂G) is found in the majority of eukaryotic tRNAs at base pair m²₂G26:A44. The latter modification eliminates the ability of the N² function to donate in hydrogen bonds and alters its pairing behavior, notably vis-à-vis C. Perhaps a less obvious consequence of the N²,N²-dimethyl modification is its role in controlling the pairing modes between G and A. We have determined the crystal structure of a 13-mer RNA duplex with central tandem m²₂G:A pairs. In the structure both pairs adopt an imino-hydrogen bonded, pseudo-Watson–Crick conformation. Thus, the sheared conformation frequently seen in tandem G:A pairs is avoided due to a potential steric clash between an N²-methyl group and the major groove edge of A. Additionally, for a series of G:A containing self-complementary RNAs we investigated how methylation affects competitive hairpin versus duplex formation based on UV melting profile analysis.     

Thermodynamics of RNA hairpins containing single internal mismatches.

Thermodynamic parameters and circular dichroism spectra are presented for RNA hairpins containing single internal mismatches in the stem regions. Three different sequence contexts for the G*U mismatch and two contexts for C*A, G*A, U*U, A*C and U*G mismatches were examined and compared with Watson-Crick base-pair stabilities. The RNA hairpins employed were a microhelix and tetraloop representing the Escherichia coli tRNAAlaacceptor stem and sequence variants that have been altered at the naturally occurring G*U mismatch site. UV melting studies were carried out under different conditions to evaluate the effects of sodium ion concentration and pH on the stability of mismatch-containing hairpins. Our main findings are that single internal mismatches exhibit a range of effects on hairpin stability. In these studies, the size and sequence of the loop and stem are shown to influence the overall stability of the RNA, and have a minor effect on the relative mismatch stabilities. The relationship of these results to RNA-ligand interactions involving mismatch base-pairs is discussed.     

Effects of GA mismatches on the structure and thermodynamics of RNA internal loops

UV melting, CD and NMR studies indicate rGCGAGCG and rGCAGGCG from unusually stable duplexes of type a and b. The observed delta G degree 37's in 1 M NaCl are -6.7 and -6.3 kcal/mol, respectively. For the related duplex, c, delta G degree 37 is -4.2 kcal/mol. The predicted delta G degree 37 from nearest-neighbor parameters (formula; see text) for all three duplexes is -4.7 kcal/mol (Freier, S.M., Kierzek, R., Jaeger, J.A., Sugimoto, N., Caruthers, M.H., Neilson, T., & Turner, D.H. (1986) Proc. Natl. Acad. Sci. U.S.A. 83, 9373-9377). The results suggest a special interaction in the duplexes containing GA mismatches. Presumably, this is hydrogen bonding between G and A. While the thermodynamics for (rGCGAGCG)2 and (rGCAGGCG)2 are similar, CD and the imino region of the proton NMR spectra indicate their structures are different. In particular, (rGCAGGCG)2 exhibits a CD spectrum typical of A-form geometry with a weak negative band at 280 nm. In contrast, the CD spectrum for (rGCGAGCG)2 has an intense positive band at 285 nm. The NMR spectrum of (rGCAGGCG)2 has a resonance corresponding to a hydrogen-bonded GA mismatch, while for (rGCGAGCG)2 no hydrogen-bonded imino proton is observed for the mismatch. The glycosidic torsion angles of the bases in the GA mismatches of (rGCAGGCG)2 and (rCGCAGGCG)2 are anti. Duplexes of type d, where X is A, G, or U, are more stable than e, and the stability differences are similar to those (formula; see text) observed for f versus g. Thus, 3'-dangling ends in this system make contributions to duplex stability that are similar to contributions observed with fully paired 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 Å, c = 56.41 Å, and γ = 120°, with one duplex in the asymmetric unit. The structure was solved by the molecular replacement method at 1.9 Å 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 Å with F ≥ 2σ(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 a 26-nucleotide RNA containing a hook-turn

A crystal structure has been obtained for a 26-nucleotide RNA that contains the loop E sequence from Chromatium minutissimum. Rather than having a loop E-like conformation, it consists of an A-form helix that splits into two separate strands following a sheared A-G base pair. The backbone of the strand containing the G of the A-G pair makes a turn of almost 180 degrees in the space of two nucleotides, and then interacts with the minor groove of the helix from which it originates. Similar structures, which we call hook-turns, occur in 16S and 23S rRNAs. They are found at places where the two strands of a helix separate at an A/G juxtaposition to interact with other sequences.     

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.     

Thermodynamics of single mismatches in RNA duplexes

The thermodynamic properties and structures of single mismatches in short RNA duplexes were studied in optical melting and imino proton NMR experiments. The free energy increments at 37 degrees C measured for non-GU single mismatches range from -2.6 to 1.7 kcal/mol. These increments depend on the identity of the mismatch, adjacent base pairs, and the position in the helix. UU and AA mismatches are more stable close to a helix end, but GG mismatch stability is essentially unaffected by the position in the helix. Approximations are suggested for predicting stabilities of single mismatches in short RNA duplexes.     

PDZ tandem of human syntenin: crystal structure and functional properties

Syntenin, a 33 kDa protein, interacts with several cell membrane receptors and with merlin, the product of the causal gene for neurofibromatosis type II. We report a crystal structure of the functional fragment of human syntenin containing two canonical PDZ domains, as well as binding studies for full-length syntenin, the PDZ tandem, and isolated PDZ domains. We show that the functional properties of syntenin are a result of independent interactions with target peptides, and that each domain is able to bind peptides belonging to two different classes: PDZ1 binds peptides from classes I and III, while PDZ2 interacts with classes I and II. The independent binding of merlin by PDZ1 and syndecan-4 by PDZ2 provides direct evidence for the coupling of syndecan-mediated signaling to actin regulation by merlin.     

The active form of Xp54 RNA helicase in translational repression is an RNA-mediated oligomer

Previously, we reported that in clam oocytes, cytoplasmic polyadenylation element-binding protein (CPEB) co-immunoprecipitates with p47, a member of the highly conserved RCK family of RNA helicases which includes Drosophila Me31B and Saccharomyces cerevisiae Dhh1. Xp54, the Xenopus homologue, with helicase activity, is a component of stored mRNP. In tethered function assays in Xenopus oocytes, we showed that MS2–Xp54 represses the translation of non-adenylated firefly luciferase mRNAs and that mutations in two core helicase motifs, DEAD and HRIGR, surprisingly, activated translation. Here we show that wild-type MS2–Xp54 tethered to the reporter mRNA 3′-untranslated region (UTR) represses translation in both oocytes and eggs in an RNA-dependent complex with endogenous Xp54. Injection of mutant helicases or adenylated reporter mRNA abrogates this association. Thus Xp54 oligomerization is a hallmark of translational repression. Xp54 complexes, which also contain CPEB and eIF4E in oocytes, change during meiotic maturation. In eggs, CPEB is degraded and, while eIF4E still interacts with Xp54, this interaction becomes RNA dependent. Supporting evidence for RNA-mediated oligomerization of endogenous Xp54, and RNA-independent association with CPEB and eIF4E in oocytes was obtained by gel filtration. Altogether, our data are consistent with a model in which the active form of the Xp54 RNA helicase is an oligomer in vivo which, when tethered, via either MS2 or CPEB to the 3′UTR, represses mRNA translation, possibly by sequestering eIF4E from the translational machinery.     

The crystal and molecular structure of an osmium bispyridine adduct of thymine.

The bispyridine osmium adduct of thymine has been crystallised and subjected to an X-ray diffraction analysis. It crystallises in the triclinic space group P1, with cell dimensions a equals 7.975(3), b equals 10.381(3), c equals 11.036(3) A, alpha equals 82.73(2)degrees, beta equals 77.22(3) degrees, gamma equals 101.75(3), and with two molecules in the unit cell. The analysis has shown that the osmium reagent has added cis across the 5,6 thymine bond.     

The crystal structure of the Argonaute2 PAZ domain reveals an RNA binding motif in RNAi effector complexes

RISC, the RNA-induced silencing complex, uses short interfering RNAs (siRNAs) or micro RNAs (miRNAs) to select its targets in a sequence-dependent manner. Key RISC components are Argonaute proteins, which contain two characteristic domains, PAZ and PIWI. PAZ is highly conserved and is found only in Argonaute proteins and Dicer. We have solved the crystal structure of the PAZ domain of Drosophila Argonaute2. The PAZ domain contains a variant of the OB fold, a module that often binds single-stranded nucleic acids. PAZ domains show low-affinity nucleic acid binding, probably interacting with the 3' ends of single-stranded regions of RNA. PAZ can bind the characteristic two-base 3' overhangs of siRNAs, indicating that although PAZ may not be a primary nucleic acid binding site in Dicer or RISC, it may contribute to the specific and productive incorporation of siRNAs and miRNAs into the RNAi pathway.