Common and distinctive features of GNRA tetraloops based on a GUAA tetraloop structure at 1.4 A resolution

Structural analysis provides insights into the modular organization of picornavirus IRES

Picornavirus RNA translation is driven by the internal ribosome entry site (IRES) element. The impact of RNA structure on the foot-and-mouth disease virus (FMDV) IRES activity has been analyzed using Selective 2'Hydroxyl Acylation analyzed by Primer Extension (SHAPE) and high throughput analysis of RNA conformation by antisense oligonucleotides printed on microarrays. SHAPE reactivity revealed the self-folding capacity of domain 3 and evidenced a change of RNA structure in a defective GNRA mutant. A modified RNA conformation of this mutant was also evidenced by RNA accessibility to oligonucleotides. Interestingly, comparison of nucleotide reactivity with RNA accessibility revealed that SHAPE reactive nucleotides corresponding to the GNRA motif were not accessible to their respective target oligonucleotides. The differential response was observed both in domain 3 and the entire IRES. Our results demonstrate distant effects of the GNRA motif in the domain 3 RNA conformation, and highlight the modular organization of a picornavirus IRES.

Coarse-grained model for simulation of RNA three-dimensional structures

The accurate prediction of an RNA's three-dimensional structure from its "primary structure" will have a tremendous influence on the experimental design and its interpretation and ultimately our understanding of the many functions of RNA. This paper presents a general coarse-grained (CG) potential for modeling RNA 3-D structures. Each nucleotide is represented by five pseudo atoms, two for the backbone (one for the phosphate and another for the sugar) and three for the base to represent base-stacking interactions. The CG potential has been parametrized from statistical analysis of 688 RNA experimental structures. Molecular dynamic simulations of 15 RNA molecules with the length of 12-27 nucleotides have been performed using the CG potential, with performance comparable to that from all-atom simulations. For ~75% of systems tested, simulated annealing led to native-like structures at least once out of multiple repeated runs. Furthermore, with weak distance restraints based on the knowledge of three to five canonical Watson-Crick pairs, all 15 RNAs tested are successfully folded to within 6.5 Å of native structures using the CG potential and simulated annealing. The results reveal that with a limited secondary structure model the current CG potential can reliably predict the 3-D structures for small RNA molecules. We also explored an all-atom force field to construct atomic structures from the CG simulations.

Turning limited experimental information into 3D models of RNA

Our understanding of RNA functions in the cell is evolving rapidly. As for proteins, the detailed three-dimensional (3D) structure of RNA is often key to understanding its function. Although crystallography and nuclear magnetic resonance (NMR) can determine the atomic coordinates of some RNA structures, many 3D structures present technical challenges that make these methods difficult to apply. The great flexibility of RNA, its charged backbone, dearth of specific surface features, and propensity for kinetic traps all conspire with its long folding time, to challenge in silico methods for physics-based folding. On the other hand, base-pairing interactions (either in runs to form helices or isolated tertiary contacts) and motifs are often available from relatively low-cost experiments or informatics analyses. We present RNABuilder, a novel code that uses internal coordinate mechanics to satisfy user-specified base pairing and steric forces under chemical constraints. The code recapitulates the topology and characteristic L-shape of tRNA and obtains an accurate noncrystallographic structure of the Tetrahymena ribozyme P4/P6 domain. The algorithm scales nearly linearly with molecule size, opening the door to the modeling of significantly larger structures.

The RNA backbone plays a crucial role in mediating the intrinsic stability of the GpU dinucleotide platform and the GpUpA/GpA miniduplex

The side-by-side interactions of nucleobases contribute to the organization of RNA, forming the planar building blocks of helices and mediating chain folding. Dinucleotide platforms, formed by side-by-side pairing of adjacent bases, frequently anchor helices against loops. Surprisingly, GpU steps account for over half of the dinucleotide platforms observed in RNA-containing structures. Why GpU should stand out from other dinucleotides in this respect is not clear from the single well-characterized H-bond found between the guanine N2 and the uracil O4 groups. Here, we describe how an RNA-specific H-bond between O2′(G) and O2P(U) adds to the stability of the GpU platform. Moreover, we show how this pair of oxygen atoms forms an out-of-plane backbone ‘edge’ that is specifically recognized by a non-adjacent guanine in over 90% of the cases, leading to the formation of an asymmetric miniduplex consisting of ‘complementary’ GpUpA and GpA subunits. Together, these five nucleotides constitute the conserved core of the well-known loop-E motif. The backbone-mediated intrinsic stabilities of the GpU dinucleotide platform and the GpUpA/GpA miniduplex plausibly underlie observed evolutionary constraints on base identity. We propose that they may also provide a reason for the extreme conservation of GpU observed at most 5′-splice sites.

Transposable elements in disease-associated cryptic exons

Transposable elements (TEs) make up a half of the human genome, but the extent of their contribution to cryptic exon activation that results in genetic disease is unknown. Here, a comprehensive survey of 78 mutation-induced cryptic exons previously identified in 51 disease genes revealed the presence of TEs in 40 cases (51%). Most TE-containing exons were derived from short interspersed nuclear elements (SINEs), with Alus and mammalian interspersed repeats (MIRs) covering >18 and >16% of the exonized sequences, respectively. The majority of SINE-derived cryptic exons had splice sites at the same positions of the Alu/MIR consensus as existing SINE exons and their inclusion in the mRNA was facilitated by phylogenetically conserved changes that improved both traditional and auxiliary splicing signals, thus marking intronic TEs amenable for pathogenic exonization. The overrepresentation of MIRs among TE exons is likely to result from their high average exon inclusion levels, which reflect their strong splice sites, a lack of splicing silencers and a high density of enhancers, particularly (G)AA(G) motifs. These elements were markedly depleted in antisense Alu exons, had the most prominent position on the exon-intron gradient scale and are proposed to promote exon definition through enhanced tertiary RNA interactions involving unpaired (di)adenosines. The identification of common mechanisms by which the most dynamic parts of the genome contribute both to new exon creation and genetic disease will facilitate detection of intronic mutations and the development of computational tools that predict TE hot-spots of cryptic exon activation.

Viroid pathogenicity: one process, many faces

Despite the non-coding nature of their small RNA genomes, the visible symptoms of viroid infection resemble those associated with many plant virus diseases. Recent evidence indicates that viroid-derived small RNAs acting through host RNA silencing pathways play a key role in viroid pathogenicity. Host responses to viroid infection are complex, involving signaling cascades containing host-encoded protein kinases and crosstalk between hormonal and defense-signaling pathways. Studies of viroid-host interaction in the context of entire biochemical or developmental pathways are just beginning, and many working hypotheses have yet to be critically tested.

Extensive backbone dynamics in the GCAA RNA tetraloop analyzed using 13C NMR spin relaxation and specific isotope labeling

Conformational dynamics play a key role in the properties and functions of proteins and nucleic acids. Heteronuclear NMR spin relaxation is a uniquely powerful site-specific probe of dynamics in proteins and has found increasing applications to nucleotide base side chains and anomeric sites in RNA. Applications to the nucleic acid ribose backbone, however, have been hampered by strong magnetic coupling among ring carbons in uniformly 13C-labeled samples. In this work, we apply a recently developed, metabolically directed isotope labeling scheme that places 13C with high efficiency and specificity at the nucleotide ribose C2' and C4' sites. We take advantage of this scheme to explore backbone dynamics in the well-studied GCAA RNA tetraloop. Using a combination of CPMG (Carr-Purcell-Meiboom-Gill) and R(1rho) relaxation dispersion spectroscopy to explore exchange processes on the microsecond to millisecond time scale, we find an extensive pattern of dynamic transitions connecting a set of relatively well-defined conformations. In many cases, the observed transitions appear to be linked to C3'-endo/C2'-endo sugar pucker transitions of the corresponding nucleotides, and may also be correlated across multiple nucleotides within the tetraloop. These results demonstrate the power of NMR spin relaxation based on alternate-site isotope labeling to open a new window into the dynamic properties of ribose backbone groups in RNA.

Nucleic acid folding determined by mesoscale modeling and NMR spectroscopy: solution structure of d(GCGAAAGC)

Determination of DNA solution structure is a difficult task even with the high-sensitivity method used here based on simulated annealing with 35 restraints/residue (Cryoprobe 750 MHz NMR). The conformations of both the phosphodiester linkages and the dinucleotide segment encompassing the sharp turn in single-stranded DNA are often underdetermined. To obtain higher quality structures of a DNA GNRA loop, 5'-d(GCGAAAGC)-3', we have used a mesoscopic molecular modeling approach, called Biopolymer Chain Elasticity (BCE), to provide reference conformations. By construction, these models are the least deformed hairpin loop conformation derived from canonical B-DNA at the nucleotide level. We have further explored this molecular conformation at the torsion angle level with AMBER molecular mechanics using different possible (epsilon,zeta) constraints to interpret the 31P NMR data. This combined approach yields a more accurate molecular conformation, compatible with all the NMR data, than each method taken separately, NMR/DYANA or BCE/AMBER. In agreement with the principle of minimal deformation of the backbone, the hairpin motif is stabilized by maximal base-stacking interactions on both the 5'- and 3'-sides and by a sheared G.A mismatch base pair between the first and last loop nucleotides. The sharp turn is located between the third and fourth loop nucleotides, and only two torsion angles beta6 and gamma6 deviate strongly with respect to canonical B-DNA structure. Two other torsion angle pairs epsilon3,zeta3 and epsilon5,zeta5 exhibit the newly recognized stable conformation BIIzeta+ (-70 degrees, 140 degrees). This combined approach has proven to be useful for the interpretation of an unusual 31P chemical shift in the 5'-d(GCGAAAGC)-3' hairpin.

Trans Hoogsteen/sugar edge base pairing in RNA. Structures, energies, and stabilities from quantum chemical calculations

Trans Hoogsteen/sugar edge (H/SE) RNA base pairs form one of the six families of RNA base pairs that utilize the 2'-hydroxyl group of ribose for base pairing and play key roles in stabilizing folded RNA molecules. Here, we provide a detailed quantum chemical characterization of intrinsic structures and interaction energies of this base pair family, along with the evaluation of solvent screening effects by a continuum solvent approach. We report DFT-optimized geometries and MP2 interaction energies for all 10 crystallographically identified members of the family, for a representative set of them, using complete basis set extrapolation. For 6 of the 10 base pairs, we had to apply geometric constraints to keep the geometries relevant to RNA. We confirm that the remaining, hitherto undetected, possible members of this family do not have appropriate steric features required to establish stable base pairing in the trans H/SE fashion. The interaction patterns in the trans H/SE family are highly diverse, with gas-phase interaction energies in the range from -1 to -17 kcal/mol. Except for the C/rC and G/rG trans H/SE base pairs, the interaction energy is roughly evenly distributed between the HF and correlation components. Thus, in the trans H/SE base pairs, the relative importance of electron correlation is noticeably smaller than in the cis WC/SE or cis and trans SE/SE base pairs, but still larger than in canonical base pairs. The trans H/SE A/rG base pair is the intrinsically most stable member of this family. This base pair is also known as the sheared AG base pair and belongs to the most prominent set of RNA base pairs utilized in molecular building blocks of functional RNAs. For all trans H/SE base pairs that we identified, in addition to conventional base pairing, viable alternative structures were stabilized by amino-acceptor interactions. In the QM calculations, these amino-acceptor complexes appear to be equally as stable as those with common H-bonds, and more importantly, the switch to amino-acceptor interaction does not require any significant geometrical rearrangement of the base pairs. Such interactions are worthy of further investigations, as X-ray crystallography cannot unambiguously distinguish between conventional and amino-acceptor interactions involving the 2'-hydroxyl group, formation of such interactions is usually not considered, and molecular modeling force fields do not include such interactions properly as a result of neglect of aminogroup pyramidalization.

Finding 3D motifs in ribosomal RNA structures

The identification of small structural motifs and their organization into larger subassemblies is of fundamental interest in the analysis, prediction and design of 3D structures of large RNAs. This problem has been studied only sparsely, as most of the existing work is limited to the characterization and discovery of motifs in RNA secondary structures. We present a novel geometric method for the characterization and identification of structural motifs in 3D rRNA molecules. This method enables the efficient recognition of known 3D motifs, such as tetraloops, E-loops, kink-turns and others. Furthermore, it provides a new way of characterizing complex 3D motifs, notably junctions, that have been defined and identified in the secondary structure but have not been analyzed and classified in three dimensions. We demonstrate the relevance and utility of our approach by applying it to the Haloarcula marismortui large ribosomal unit. Pending the implementation of a dedicated web server, the code accompanying this article, written in JAVA, is available upon request from the contact author.

5' cis elements direct nodavirus RNA1 recruitment to mitochondrial sites of replication complex formation

Positive-strand RNA viruses replicate their genomes on intracellular membranes, usually in conjunction with virus-induced membrane rearrangements. For the nodavirus flock house virus (FHV), we recently showed that multifunctional FHV replicase protein A induces viral RNA template recruitment to a membrane-associated state, but the site(s) and function of this recruitment were not determined. By tagging viral RNA with green fluorescent protein, we show here in Drosophila cells that protein A recruits FHV RNA specifically to the outer mitochondrial membrane sites of RNA replication complex formation. Using Drosophila cells and yeast cells, which also support FHV replication, we also defined the cis-acting regions that direct replication and template recruitment for FHV genomic RNA1. RNA1 nucleotides 68 to 205 were required for RNA replication and directed efficient protein A-mediated RNA recruitment in both cell types. RNA secondary structure prediction, structure probing, and phylogenetic comparisons in this region identified two stable, conserved stem-loops with nearly identical loop sequences. Further mutational analysis showed that both stem-loops and certain flanking sequences were required for RNA1 recruitment, negative-strand synthesis, and subsequent positive-strand amplification in yeast and Drosophila cells. Thus, we have shown that protein A recruits RNA1 templates to mitochondria, as expected for RNA replication, and identified a new RNA1 cis element that is necessary and sufficient for RNA1 template recognition and recruitment to these mitochondrial membranes for negative-strand RNA1 synthesis. These results establish RNA recruitment to the sites of replication complex formation as an essential, distinct, and selective early step in nodavirus replication.

The GANC tetraloop: a novel motif in the group IIC intron structure

Tetraloops are a common building block for RNA tertiary structure, and most tetraloops fall into one of three well-characterized classes: GNRA, UNCG, and CUYG. Here, we present the sequence and structure of a fourth highly conserved class of tetraloop that occurs only within the zeta-zeta' interaction of group IIC introns. This GANC tetraloop was identified, along with an unusual cognate receptor, in the crystal structure of the group IIC intron and through phylogenetic analysis of intron RNA sequence alignments. Unlike conventional tetraloop-receptor interactions, which are stabilized by extensive hydrogen-bonding interactions, the GANC-receptor interaction is limited to a single base stack between the conserved adenosine of the tetraloop and a single purine of the receptor, which consists of a one- to three-nucleotide bulge and does not contain an A-platform. Unlike GNRA tetraloops, the GANC tetraloop forms a sharp angle relative to the adjacent helix, bending by approximately 45 degrees toward the major groove side of the helix. These structural attributes allow GANC tetraloops to fit precisely within the group IIC intron core, thereby demonstrating that structural motifs can adapt to function in a specific niche.

Electrostatic interactions guide the active site face of a structure-specific ribonuclease to its RNA substrate

Restrictocin, a member of the α-sarcin family of site-specific endoribonucleases, uses electrostatic interactions to bind to the ribosome and to RNA oligonucleotides, including the minimal specific substrate, the sarcin/ricin loop (SRL) of 23S−28S rRNA. Restrictocin binds to the SRL by forming a ground-state E:S complex that is stabilized predominantly by Coulomb interactions and depends on neither the sequence nor structure of the RNA, suggesting a nonspecific complex. The 22 cationic residues of restrictocin are dispersed throughout this protein surface, complicating a priori identification of a Coulomb interacting surface. Structural studies have identified an enzyme−substrate interface, which is expected to overlap with the electrostatic E:S interface. Here, we identified restrictocin residues that contribute to binding in the E:S complex by determining the salt dependence [∂ log(k₂/K_1/2)/∂ log[KCl]] of cleavage of the minimal SRL substrate for eight point mutants within the protein designed to disrupt contacts in the crystallographically defined interface. Relative to the wild-type salt dependence of −4.1, a subset of the mutants clustering near the active site shows significant changes in salt dependence, with differences of magnitude being ≥0.4. This same subset was identified using calculated salt dependencies for each mutant derived from solutions to the nonlinear Poisson−Boltzmann equation. Our findings support a mechanism in which specific residues on the active site face of restrictocin (primarily K110, K111, and K113) contribute to formation of the E:S complex, thereby positioning the SRL substrate for site-specific cleavage. The same restrictocin residues are expected to facilitate targeting of the SRL on the surface of the ribosome.

Bacteriophage P22 antitermination boxB sequence requirements are complex and overlap with those of lambda

Transcription antitermination in phages lambda and P22 uses N proteins that bind to similar boxB RNA hairpins in regulated transcripts. In contrast to the lambda N-boxB interaction, the P22 N-boxB interaction has not been extensively studied. A nuclear magnetic resonance structure of the P22 N peptide boxB(left) complex and limited mutagenesis have been reported but do not reveal a consensus sequence for boxB. We have used a plasmid-based antitermination system to screen boxBs with random loops and to test boxB mutants. We find that P22 N requires boxB to have a GNRA-like loop with no simple requirements on the remaining sequences in the loop or stem. U:A or A:U base pairs are strongly preferred adjacent to the loop and appear to modulate N binding in cooperation with the loop and distal stem. A few GNRA-like hexaloops have moderate activity. Some boxB mutants bind P22 and lambda N, indicating that the requirements imposed on boxB by P22 N overlap those imposed by lambda N. Point mutations can dramatically alter boxB specificity between P22 and lambda N. A boxB specific for P22 N can be mutated to lambda N specificity by a series of single mutations via a bifunctional intermediate, as predicted by neutral theories of evolution.

The impact of RNA structure on picornavirus IRES activity

Internal ribosome entry site (IRES) elements consist of cis-acting regions that recruit the translation machinery to an internal position in the mRNA. The biological relevance of RNA structure-mediated mechanisms involved in internal ribosome recruitment is now emerging from the structural and functional analysis of viral IRES elements. However, because IRES elements found in genetically distant mRNAs seem to be organized in different RNA structures, the definition of the structural requirements for IRES activity is challenging and demands multidisciplinary approaches. This review discusses the latest reports that establish a relationship between RNA structure and IRES function in picornavirus genomes, the first RNAs described to contain these specialized regulatory elements.

DNA mimicry by a high-affinity anti-NF-kappaB RNA aptamer

The binding of RNA molecules to proteins or other ligands can require extensive RNA folding to create an induced fit. Understanding the generality of this principle involves comparing structures of RNA before and after complex formation. Here we report the NMR solution structure of a 29-nt RNA aptamer whose crystal structure had previously been determined in complex with its transcription factor target, the p50₂ form of NF-κB. The RNA aptamer internal loop structure has pre-organized features that are also found in the complex, including non-canonical base pairing and cross-strand base stacking. Remarkably, the free RNA aptamer structure possesses a major groove that more closely resembles B-form DNA than RNA. Upon protein binding, changes in RNA structure include the kinking of the internal loop and distortion of the terminal tetraloop. Thus, complex formation involves both pre-formed and induced fit binding interactions. The high affinity of the NF-κB transcription factor for this RNA aptamer may largely be due to the structural pre-organization of the RNA that results in its ability to mimic DNA.

Comprehensive features of natural and in vitro selected GNRA tetraloop-binding receptors

Specific recognitions of GNRA tetraloops by small helical receptors are among the most widespread long-range packing interactions in large ribozymes. However, in contrast to GYRA and GAAA tetraloops, very few GNRA/receptor interactions have yet been identified to involve GGAA tetraloops in nature. A novel in vitro selection scheme based on a rigid self-assembling tectoRNA scaffold designed for isolation of intermolecular interactions with A-minor motifs has yielded new GGAA tetraloop-binding receptors with affinity in the nanomolar range. One of the selected receptors is a novel 12 nt RNA motif, (CCUGUG … AUCUGG), that recognizes GGAA tetraloop hairpin with a remarkable specificity and affinity. Its physical and chemical characteristics are comparable to those of the well-studied ‘11nt’ GAAA tetraloop receptor motif. A second less specific motif (CCCAGCCC … GAUAGGG) binds GGRA tetraloops and appears to be related to group IC3 tetraloop receptors. Mutational, thermodynamic and comparative structural analysis suggests that natural and in vitro selected GNRA receptors can essentially be grouped in two major classes of GNRA binders. New insights about the evolution, recognition and structural modularity of GNRA and A-minor RNA–RNA interactions are proposed.

Probing the structural hierarchy and energy landscape of an RNA T-loop hairpin

The T-loop motif is an important recurrent RNA structural building block consisting of a U-turn sub-motif and a UA trans Watson–Crick/Hoogsteen base pair. In the presence of a hairpin stem, the UA non-canonical base pair becomes part of the UA-handle motif. To probe the hierarchical organization and energy landscape of the T-loop, we performed replica exchange molecular dynamics (REMD) simulations of the T-loop in isolation and as part of a hairpin. Our simulations reveal that the isolated T-loop adopts coil conformers stabilized by base stacking. The T-loop hairpin shows a highly rugged energy landscape featuring multiple local minima with a transition state for folding consisting of partially zipped states. The U-turn displays a high conformational flexibility both when the T-loop is in isolation and as part of a hairpin. On the other hand, the stability of the UA non-canonical base pair is enhanced in the presence of the UA-handle. This motif is apparently a key component for stabilizing the T-loop, while the U-turn is mostly involved in long-range interaction. Our results suggest that the stability and folding of small RNA motifs are highly dependent on local context.

Identification of an aminoacridine derivative that binds to RNA tetraloops

RNA folds into diverse structures that form unique targets for small molecules and thus provide significant potential for controlling biological processes involving RNA with small-molecule ligands. We are investigating molecular recognition of tetraloop RNA by small molecules. RNA tetraloops are four-nucleotide stem loops with unusual stability that are involved in biological processes involving RNA by forming binding sites for proteins and other RNAs. We have sequentially used the docking programs DOCK and AutoDock to screen 1990 small molecules in the NCI diversity set to identify molecules selective for RNA tetraloops over double-stranded RNA. The compounds predicted to bind to tetraloop RNA were evaluated for binding RNA tetraloops using 1H NMR spectroscopy and fluorescence techniques. An aminoacridine derivative (AD2) was identified that binds to a GAAA tetraloop in a 2:1 ratio with dissociation constants of 1.0 and 4.0 microM. AD2 binds with approximately 20-fold and 9-fold higher affinity to tetraloop RNA than to double- and single-stranded RNAs, respectively.

Fungal ribotoxins: molecular dissection of a family of natural killers

RNase T1 is the best known representative of a large family of ribonucleolytic proteins secreted by fungi, mostly Aspergillus and Penicillium species. Ribotoxins stand out among them by their cytotoxic character. They exert their toxic action by first entering the cells and then cleaving a single phosphodiester bond located within a universally conserved sequence of the large rRNA gene, known as the sarcin-ricin loop. This cleavage leads to inhibition of protein biosynthesis, followed by cellular death by apoptosis. Although no protein receptor has been found for ribotoxins, they preferentially kill cells showing altered membrane permeability, such as those that are infected with virus or transformed. Many steps of the cytotoxic process have been elucidated at the molecular level by means of a variety of methodological approaches and the construction and purification of different mutant versions of these ribotoxins. Ribotoxins have been used for the construction of immunotoxins, because of their cytotoxicity. Besides this activity, Aspf1, a ribotoxin produced by Aspergillus fumigatus, has been shown to be one of the major allergens involved in allergic aspergillosis-related pathologies. Protein engineering and peptide synthesis have been used in order to understand the basis of these pathogenic mechanisms as well as to produce hypoallergenic proteins with potential diagnostic and immunotherapeutic applications.

The UAA/GAN internal loop motif: a new RNA structural element that forms a cross-strand AAA stack and long-range tertiary interactions

Analysis of aligned RNA sequences and high-resolution crystal structures has revealed a new RNA structural element, termed the UAA/GAN motif. Found in internal loops of the 23 S rRNA, as well as in RNase P RNA and group I and II introns, this six-nucleotide motif adopts a distinctive local structure that includes two base-pairs with non-canonical conformations and three conserved adenine bases, which form a cross-strand AAA stack in the minor groove. Most importantly, the motif invariably forms long-range tertiary contacts, as the AAA stack typically forms A-minor interactions and the flipped-out N nucleotide forms additional contacts that are specific to the structural context of each loop. The widespread presence of this motif and its propensity to form long-range contacts suggest that it plays a critical role in defining the architectures of structured RNAs.

Single nucleotide RNA choreography

New structural analysis methods, and a tree formalism re-define and expand the RNA motif concept, unifying what previously appeared to be disparate groups of structures. We find RNA tetraloops at high frequencies, in new contexts, with unexpected lengths, and in novel topologies. The results, with broad implications for RNA structure in general, show that even at this most elementary level of organization, RNA tolerates astounding variation in conformation, length, sequence and context. However the variation is not random; it is well-described by four distinct modes, which are 3-2 switches (backbone topology variations), insertions, deletions and strand clips.

Molecular dynamics simulations of sarcin-ricin rRNA motif

Explicit solvent molecular dynamics (MD) simulations were carried out for sarcin-ricin domain (SRD) motifs from 23S (Escherichia coli) and 28S (rat) rRNAs. The SRD motif consists of GAGA tetraloop, G-bulged cross-strand A-stack, flexible region and duplex part. Detailed analysis of the overall dynamics, base pairing, hydration, cation binding and other SRD features is presented. The SRD is surprisingly static in multiple 25 ns long simulations and lacks any non-local motions, with root mean square deviation (r.m.s.d.) values between averaged MD and high-resolution X-ray structures of 1-1.4 A. Modest dynamics is observed in the tetraloop, namely, rotation of adenine in its apex and subtle reversible shift of the tetraloop with respect to the adjacent base pair. The deformed flexible region in low-resolution rat X-ray structure is repaired by simulations. The simulations reveal few backbone flips, which do not affect positions of bases and do not indicate a force field imbalance. Non-Watson-Crick base pairs are rigid and mediated by long-residency water molecules while there are several modest cation-binding sites around SRD. In summary, SRD is an unusually stiff rRNA building block. Its intrinsic structural and dynamical signatures seen in simulations are strikingly distinct from other rRNA motifs such as Loop E and Kink-turns.

Evidence of reciprocal tertiary interactions between conserved motifs involved in organizing RNA structure essential for internal initiation of translation

Internal ribosome entry site (IRES) elements consist of highly structured RNA regions that determine internal initiation of translation. We have previously shown that the foot-and-mouth disease virus (FMDV) IRES contains a GNRA tetraloop spanning residues G178UAA181. Here we show that tertiary RNA interactions dependent on the GNRA motif determine the structural organization of the central domain. By using mutational analysis in combination with RNA probing, we have identified distant reciprocal interactions between the GNRA motif and the invariant region G240CACG244, termed motif A. Mutations in motif A caused a decrease in IRES activity as severe as the GUAG substitution in the GNRA motif. Substitutions in either GNRA or motif A sequences induced a common reorganization around the conserved R199AAA202 stem-loop, suggesting that the latter contributes to stabilize the GNRA-motif A interaction. This finding was also consistent with a significant increase in the efficiency of RNA-RNA interactions determined in gel shift assays using as probe the hairpin that contains the GNRA motif compared to a transcript encompassing the entire apical region of the central domain. Thus, we propose that the central domain of the FMDV IRES contains a structural conformation essential for IRES activity stabilized by a tertiary contact involving residues in the GNRA tetraloop and motif A conserved sequences.

Assessment of adenyl residue reactivity within model nucleic acids by surface enhanced Raman spectroscopy

We rank the reactivity of the adenyl residues (A) of model DNA and RNA molecules with electropositive subnano size [Ag]n+ sites as a function of nucleic acid primary sequences and secondary structures and in the presence of biological amounts of Cl- and Na+ or Mg2+ ions. In these conditions A is markedly more reactive than any other nucleic acid bases. A reactivity is higher in ribo (r) than in deoxyribo (d) species [pA>pdA and (pA)n>>(pdA)n]. Base pairing decreases A reactivity in corresponding duplexes but much less in r than in d. In linear single and paired dCAG or dGAC loci, base stacking inhibits A reactivity even if A is bulged or mispaired (A.A). dA tracts are highly reactive only when dilution prevents self-association and duplex structures. In d hairpins the solvent-exposed A residues are reactive in CAG and GAC triloops and even more in ATC loops. Among the eight rG1N2R3A4 loops, those bearing a single A (A4) are the least reactive. The solvent-exposed A2 is reactive, but synergistic structural transitions make the initially stacked A residues of any rGNAA loop much more reactive. Mg2+ cross-bridging single strands via phosphates may screen A reactivity. In contrast d duplexes cross-bridging enables "A flipping" much more in rA.U pairs than in dA.T. Mg2+ promotes A reactivity in unpaired strands. For hairpins Mg2+ binding stabilizes the stems, but according to A position in the loops, A reactivity may be abolished, reduced, or enhanced. It is emphasized that not only accessibility but also local flexibility, concerted docking, and cation and anion binding control A reactivity.

Predicting RNA folding thermodynamics with a reduced chain representation model

Based on the virtual bond representation for the nucleotide backbone, we develop a reduced conformational model for RNA. We use the experimentally measured atomic coordinates to model the helices and use the self-avoiding walks in a diamond lattice to model the loop conformations. The atomic coordinates of the helices and the lattice representation for the loops are matched at the loop-helix junction, where steric viability is accounted for. Unlike the previous simplified lattice-based models, the present virtual bond model can account for the atomic details of realistic three-dimensional RNA structures. Based on the model, we develop a statistical mechanical theory for RNA folding energy landscapes and folding thermodynamics. Tests against experiments show that the theory can give much more improved predictions for the native structures, the thermal denaturation curves, and the equilibrium folding/unfolding pathways than the previous models. The application of the model to the P5abc region of Tetrahymena group I ribozyme reveals the misfolded intermediates as well as the native-like intermediates in the equilibrium folding process. Moreover, based on the free energy landscape analysis for each and every loop mutation, the model predicts five lethal mutations that can completely alter the free energy landscape and the folding stability of the molecule.

Exploring the complex folding kinetics of RNA hairpins: II. Effect of sequence, length, and misfolded states

The complexity of RNA hairpin folding arises from the interplay between the loop formation, the disruption of the slow-breaking misfolded states, and the formation of the slow-forming native base stacks. We investigate the general physical mechanism for the dependence of the RNA hairpin folding kinetics on the sequence and the length of the hairpin loop and the helix stem. For example, 1), the folding would slow down when a stable GC basepair moves to the middle of the stem; 2), hairpin with GC basepair near the loop would fold/unfold faster than the one with GC near the tail of the stem; 3), within a certain range of the stem length, a longer stem can cause faster folding; and 4), certain misfolded states can assist folding through the formation of scaffold structures to lower the entropic barrier for the folding. All our findings are directly applicable and quantitatively testable in experiments. In addition, our results can be useful for molecular design to achieve desirable fast/slow-folding hairpins, hairpins with/without specific misfolded intermediates, and hairpins that fold along designed pathways.

A simple motif for protein recognition in DNA secondary structures

DNA in a single-stranded form (ssDNA) exists transiently within the cell and comprises the telomeres of linear chromosomes and the genomes of some DNA viruses. As with RNA, in the single-stranded state, some DNA sequences are able to fold into complex secondary and tertiary structures that may be recognized by proteins and participate in gene regulation. To better understand how such DNA elements might fold and interact with proteins, and to compare recognition features to those of a structured RNA, we used in vitro selection to identify ssDNAs that bind an RNA-binding peptide from the HIV Rev protein with high affinity and specificity. The large majority of selected binders contain a non-Watson-Crick G.T base-pair and an adjacent C:G base-pair and both are essential for binding. This GT motif can be presented in different DNA contexts, including a nearly perfect duplex and a branched three-helix structure, and appears to be recognized in large part by arginine residues separated by one turn of an alpha-helix. Interestingly, a very similar GT motif is necessary also for protein binding and function of a well-characterized model ssDNA regulatory element from the proenkephalin promoter.

Structural insights into SRP RNA: an induced fit mechanism for SRP assembly

Proper assembly of large protein-RNA complexes requires sequential binding of the proteins to the RNA. The signal recognition particle (SRP) is a multiprotein-RNA complex responsible for the cotranslational targeting of proteins to biological membranes. Here we describe the crystal structure at 2.6-A resolution of the S-domain of SRP RNA from the archeon Methanococcus jannaschii. Comparison of this structure with the SRP19-bound form reveals the nature of the SRP19-induced conformational changes, which promote subsequent SRP54 attachment. These structural changes are initiated at the SRP19 binding site and transmitted through helix 6 to looped-out adenosines, which form tertiary RNA interaction with helix 8. Displacement of these adenosines enforces a conformational change of the asymmetric loop structure in helix 8. In free RNA, the three unpaired bases A195, C196, and C197 are directed toward the helical axis, whereas upon SRP19 binding the loop backbone inverts and the bases are splayed out in a conformation that resembles the SRP54-bound form. Nucleotides adjacent to the bulged nucleotides seem to be particularly important in the regulation of this loop transition. Binding of SRP19 to 7S RNA reveals an elegant mechanism of how protein-induced changes are directed through an RNA molecule and may relate to those regulating the assembly of other RNPs.

Nucleic acid crystallography: current progress

Fifty years after the publication of the DNA double helix model by Watson and Crick, new nucleic acid structures keep emerging at an ever-increasing rate. The past three years have brought a flurry of new oligonucleotide structures, including those of a Hoogsteen-paired DNA duplex, Holliday junctions, DNA-drug complexes, quadruplexes, a host of RNA motifs and various nucleic acid analogues. Major advances were also made in terms of the structure and function of catalytic RNAs. These range from improved models of the phosphodiester cleavage reactions catalyzed by the hairpin and hepatitis delta virus ribozymes to the visualization of a complete active site of a group I self-splicing intron with bound 5'- and 3'-exons. These triumphs are complemented by a refined understanding of cation-nucleic-acid interactions and new routes to the generation of derivatives for phasing of DNA and RNA structures.

A novel cGUUAg tetraloop structure with a conserved yYNMGg-type backbone conformation from cloverleaf 1 of bovine enterovirus 1 RNA

The 5'-terminal cloverleaf (CL)-like RNA structures are essential for the initiation of positive- and negative-strand RNA synthesis of entero- and rhinoviruses. SLD is the cognate RNA ligand of the viral proteinase 3C (3C(pro)), which is an indispensable component of the viral replication initiation complex. The structure of an 18mer RNA representing the apical stem and the cGUUAg D-loop of SLD from the first 5'-CL of BEV1 was determined in solution to a root-mean-square deviation (r.m.s.d.) (all heavy atoms) of 0.59 A (PDB 1Z30). The first (antiG) and last (synA) nucleotide of the D-loop forms a novel 'pseudo base pair' without direct hydrogen bonds. The backbone conformation and the base-stacking pattern of the cGUUAg-loop, however, are highly similar to that of the coxsackieviral uCACGg D-loop (PDB 1RFR) and of the stable cUUCGg tetraloop (PDB 1F7Y) but surprisingly dissimilar to the structure of a cGUAAg stable tetraloop (PDB 1MSY), even though the cGUUAg BEV D-loop and the cGUAAg tetraloop differ by 1 nt only. Together with the presented binding data, these findings provide independent experimental evidence for our model [O. Ohlenschlager, J. Wohnert, E. Bucci, S. Seitz, S. Hafner, R. Ramachandran, R. Zell and M. Gorlach (2004) Structure, 12, 237-248] that the proteinase 3C(pro) recognizes structure rather than sequence.

Binding of an aminoacridine derivative to a GAAA RNA tetraloop

RNA tetraloops are common secondary structural motifs in many RNAs, especially ribosomal RNAs. There are few studies of small molecule recognition of RNA tetraloops although tetraloops are known to interact with RNA receptors and proteins, and to form nucleation sites for RNA folding. In this paper, we investigate the binding of neomycin, kanamycin, 2,4-diaminoquinazoline, quinacrine, and an aminoacridine derivative (AD1) to a GAAA tetraloop using fluorescence spectroscopy. We have found that AD1 and quinacrine bind to the GAAA tetraloop with the highest affinity of the molecules examined. The equilibrium dissociation constant of the AD1-GAAA tetraloop complex was determined to be 1.6 microM. RNase I and lead acetate footprinting experiments suggested that AD1 binds to the junction between the loop and stem of the GAAA tetraloop.

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.

The structure of a rigorously conserved RNA element within the SARS virus genome

We have solved the three-dimensional crystal structure of the stem-loop II motif (s2m) RNA element of the SARS virus genome to 2.7-Å resolution. SARS and related coronaviruses and astroviruses all possess a motif at the 3′ end of their RNA genomes, called the s2m, whose pathogenic importance is inferred from its rigorous sequence conservation in an otherwise rapidly mutable RNA genome. We find that this extreme conservation is clearly explained by the requirement to form a highly structured RNA whose unique tertiary structure includes a sharp 90° kink of the helix axis and several novel longer-range tertiary interactions. The tertiary base interactions create a tunnel that runs perpendicular to the main helical axis whose interior is negatively charged and binds two magnesium ions. These unusual features likely form interaction surfaces with conserved host cell components or other reactive sites required for virus function. Based on its conservation in viral pathogen genomes and its absence in the human genome, we suggest that these unusual structural features in the s2m RNA element are attractive targets for the design of anti-viral therapeutic agents. Structural genomics has sought to deduce protein function based on three-dimensional homology. Here we have extended this approach to RNA by proposing potential functions for a rigorously conserved set of RNA tertiary structural interactions that occur within the SARS RNA genome itself. Based on tertiary structural comparisons, we propose the s2m RNA binds one or more proteins possessing an oligomer-binding-like fold, and we suggest a possible mechanism for SARS viral RNA hijacking of host protein synthesis, both based upon observed s2m RNA macromolecular mimicry of a relevant ribosomal RNA fold.

The SARS RNA genome contains a unique structure that resembles a portion of ribosomal RNA; this may allow the virus to hijack its hosts protein synthesis machinery

“Deoxyribo Nanonucleic Acid”

The crystal structure of a single-stranded DNA oligonucleotide has revealed formation of a unique three-dimensional array by continuous antiparallel and parallel pairing between monomers. The array is based on tertiary interactions and represents a second-generation nanotechnological system.

Study of the Functional Interaction of the 900 Tetraloop of 16S Ribosomal RNA with Helix 24 within the Bacterial Ribosome

The 900 tetraloop that caps helix 27 of 16S ribosomal RNA (rRNA) is amongst the most conserved regions of rRNA. This tetraloop forms a GNRA motif that docks into the minor groove of three base-pairs at the bottom of helix 24 of 16S rRNA in the 30S subunit. Both the tetraloop and its receptor in helix 24 contact the 23S rRNA, forming the intersubunit bridge B2c. Here, we investigated the interaction between the 900 tetraloop and its receptor by genetic complementation. We used a specialized ribosome system in combination with an in vivo instant evolution approach to select mutations in helix 24 compensating for a mutation in the 900 tetraloop (A900G) that severely decreases ribosomal activity, impairing subunit association and translational fidelity. We selected two mutants where the G769-C810 base-pair of helix 24 was substituted with either U-A or C x A. When these mutations in helix 24 were investigated in the context of a wild-type 900 tetraloop, the C x A but not the U-A mutation severely impaired ribosome activity, interfering with subunit association and decreasing translational fidelity. In the presence of the A900G mutation, both mutations in helix 24 increased the ribosome activity to the same extent. Subunit association and translational fidelity were increased to the same level. Computer modeling was used to analyze the effect of the mutations in helix 24 on the interaction between the tetraloop and its receptor. This study demonstrates the functional importance of the interaction between the 900 tetraloop and helix 24.

The Neurospora crassa CYT-18 protein C-terminal RNA-binding domain helps stabilize interdomain tertiary interactions in group I introns

The Neurospora crassa mitochondrial tyrosyl-tRNA synthetase (CYT-18 protein) promotes the splicing of group I introns by stabilizing the catalytically active RNA structure. To accomplish this, CYT-18 recognizes conserved structural features of group I intron RNAs using regions of the N-terminal nucleotide-binding fold, intermediate alpha-helical, and C-terminal RNA-binding domains that also function in binding tRNA(Tyr). Curiously, whereas the splicing of the N. crassa mitochondrial large subunit rRNA intron is completely dependent on CYT-18's C-terminal RNA-binding domain, all other group I introns tested thus far are spliced efficiently by a truncated protein lacking this domain. To investigate the function of the C-terminal domain, we used an Escherichia coli genetic assay to isolate mutants of the Saccharomyces cerevisiae mitochondrial large subunit rRNA and phage T4 td introns that can be spliced in vivo by the wild-type CYT-18 protein, but not by the C-terminally truncated protein. Mutations that result in dependence on CYT-18's C-terminal domain include those disrupting two long-range GNRA tetraloop/receptor interactions: L2-P8, which helps position the P1 helix containing the 5'-splice site, and L9-P5, which helps establish the correct relative orientation of the P4-P6 and P3-P9 domains of the group I intron catalytic core. Our results indicate that different structural mutations in group I intron RNAs can result in dependence on different regions of CYT-18 for RNA splicing.

The common and the distinctive features of the bulged-G motif based on a 1.04 A resolution RNA structure

Bulged-G motifs are ubiquitous internal RNA loops that provide specific recognition sites for proteins and RNAs. To establish the common and distinctive features of the motif we determined the structures of three variants and compared them with related structures. The variants are 27-nt mimics of the sarcin/ricin loop (SRL) from Escherichia coli 23S ribosomal RNA that is an essential part of the binding site for elongation factors (EFs). The wild-type SRL has now been determined at 1.04 A resolution, supplementing data obtained before at 1.11 A and allowing the first calculation of coordinate error for an RNA motif. The other two structures, having a viable (C2658U*G2663A) or a lethal mutation (C2658G*G2663C), were determined at 1.75 and 2.25 A resolution, respectively. Comparisons reveal that bulged-G motifs have a common hydration and geometry, with flexible junctions at flanking structural elements. Six conserved nucleotides preserve the fold of the motif; the remaining seven to nine vary in sequence and alter contacts in both grooves. Differences between accessible functional groups of the lethal mutation and those of the viable mutation and wild-type SRL may account for the impaired elongation factor binding to ribosomes with the C2658G*G2663C mutation and may underlie the lethal phenotype.