RNA folding kinetics regulates translation of phage MS2 maturation gene

A roadmap for rRNA folding and assembly during transcription

Ribonucleoprotein (RNP) assembly typically begins during transcription when folding of the newly synthesized RNA is coupled with the recruitment of RNA-binding proteins (RBPs). Upon binding, the proteins induce structural rearrangements in the RNA that are crucial for the next steps of assembly. Focusing primarily on bacterial ribosome assembly, we discuss recent work showing that early RNA-protein interactions are more dynamic than previously supposed, and remain so, until sufficient proteins are recruited to each transcript to consolidate an entire domain of the RNP. We also review studies showing that stable assembly of an RNP competes against modification and processing of the RNA. Finally, we discuss how transcription sets the timeline for competing and cooperative RNA-RBP interactions that determine the fate of the nascent RNA. How this dance is coordinated is the focus of this review.

Regulation of the ThiM riboswitch is facilitated by the trapped structure formed during transcription of the wild-type sequence

The ThiM riboswitch from Escherichia coli is a typical mRNA device that modulates downstream gene expression by sensing TPP. The helix-based RNA folding theory is used to investigate its detailed regulatory behaviors in cells. This RNA molecule is transcriptionally trapped in a state with the unstructured SD sequence in the absence of TPP, which induces downstream gene expression. As a key step to turn on gene expression, formation of this trapped state (the genetic ON state) highly depends on the co-transcriptional folding of its wild-type sequence. Instead of stabilities of the genetic ON and OFF states, the transcription rate, pause, and ligand levels are combined to affect the ThiM riboswitch-mediated gene regulation, which is consistent with a kinetic control model.

Ribosome as a Translocase and Helicase

During protein synthesis, ribosome moves along mRNA to decode one codon after the other. Ribosome translocation is induced by a universally conserved protein, elongation factor G (EF-G) in bacteria and elongation factor 2 (EF-2) in eukaryotes. EF-G-induced translocation results in unwinding of the intramolecular secondary structures of mRNA by three base pairs at a time that renders the translating ribosome a processive helicase. Professor Alexander Sergeevich Spirin has made numerous seminal contributions to understanding the molecular mechanism of translocation. Here, we review Spirin's insights into the ribosomal translocation and recent advances in the field that stemmed from Spirin's pioneering work. We also discuss key remaining challenges in studies of translocase and helicase activities of the ribosome.

Uniqueness of RNA Coliphage Qβ Display System in Directed Evolutionary Biotechnology

Phage display technology involves the surface genetic engineering of phages to expose desirable proteins or peptides whose gene sequences are packaged within phage genomes, thereby rendering direct linkage between genotype with phenotype feasible. This has resulted in phage display systems becoming invaluable components of directed evolutionary biotechnology. The M13 is a DNA phage display system which dominates this technology and usually involves selected proteins or peptides being displayed through surface engineering of its minor coat proteins. The displayed protein or peptide’s functionality is often highly reduced due to harsh treatment of M13 variants. Recently, we developed a novel phage display system using the coliphage Qβ as a nano-biotechnology platform. The coliphage Qβ is an RNA phage belonging to the family of Leviviridae, a long investigated virus. Qβ phages exist as a quasispecies and possess features making them comparatively more suitable and unique for directed evolutionary biotechnology. As a quasispecies, Qβ benefits from the promiscuity of its RNA dependent RNA polymerase replicase, which lacks proofreading activity, and thereby permits rapid variant generation, mutation, and adaptation. The minor coat protein of Qβ is the readthrough protein, A₁. It shares the same initiation codon with the major coat protein and is produced each time the ribosome translates the UGA stop codon of the major coat protein with the of misincorporation of tryptophan. This misincorporation occurs at a low level (1/15). Per convention and definition, A₁ is the target for display technology, as this minor coat protein does not play a role in initiating the life cycle of Qβ phage like the pIII of M13. The maturation protein A₂ of Qβ initiates the life cycle by binding to the pilus of the F⁺ host bacteria. The extension of the A₁ protein with a foreign peptide probe recognizes and binds to the target freely, while the A₂ initiates the infection. This avoids any disturbance of the complex and the necessity for acidic elution and neutralization prior to infection. The combined use of both the A₁ and A₂ proteins of Qβ in this display system allows for novel bio-panning, in vitro maturation, and evolution. Additionally, methods for large library size construction have been improved with our directed evolutionary phage display system. This novel phage display technology allows 12 copies of a specific desired peptide to be displayed on the exterior surface of Qβ in uniform distribution at the corners of the phage icosahedron. Through the recently optimized subtractive bio-panning strategy, fusion probes containing up to 80 amino acids altogether with linkers, can be displayed for target selection. Thus, combined uniqueness of its genome, structure, and proteins make the Qβ phage a desirable suitable innovation applicable in affinity maturation and directed evolutionary biotechnology. The evolutionary adaptability of the Qβ phage display strategy is still in its infancy. However, it has the potential to evolve functional domains of the desirable proteins, glycoproteins, and lipoproteins, rendering them superior to their natural counterparts.

Self-powered RNA nanomachine driven by metastable structure

Many non-coding and regulatory RNA elements have evolved to exploit transient or metastable structures that emerge during transcription to control complex folding pathways or to encode dynamic functions. However, efforts to engineer synthetic RNA devices have mostly focused on the thermodynamically stable structures. Consequently, significant challenges and opportunities exist in engineering functional RNAs that explicitly take advantage of cotranscriptionally generated transient or metastable structures. In this work, we designed a short RNA sequence that adopts a robust metastable structure when transcribed by an RNA polymerase. Although the metastable structure persists for hours at low temperature, it refolds almost completely into the thermodynamically stable structure upon heat denaturation followed by cooling. The synthetic RNA was also equipped with the Broccoli aptamer so that it can bind its ligand and become fluorescent only in the thermodynamically stable structure. We further demonstrated that the relaxation to the thermodynamically stable and fluorescent structure can be catalyzed by a short trigger RNA in a sequence-specific manner. Finally, the RNA architecture was redesigned to sense and respond to microRNA sequences. In summary, we designed RNA nanomachines that can detect an RNA sequence, amplify signal and produce an optical output, all encoded in a single RNA transcript, self-powered by a metastable structure.

Cell-free expression of RNA encoded genes using MS2 replicase

RNA replicases catalyse transcription and replication of viral RNA genomes. Of particular interest for in vitro studies are phage replicases due to their small number of host factors required for activity and their ability to initiate replication in the absence of any primers. However, the requirements for template recognition by most phage replicases are still only poorly understood. Here, we show that the active replicase of the archetypical RNA phage MS2 can be produced in a recombinant cell-free expression system. We find that the 3' terminal fusion of antisense RNAs with a domain derived from the reverse complement of the wild type MS2 genome generates efficient templates for transcription by the MS2 replicase. The new system enables DNA-independent gene expression both in batch reactions and in microcompartments. Finally, we demonstrate that MS2-based RNA-dependent transcription-translation reactions can be used to control DNA-dependent gene expression by encoding a viral DNA-dependent RNA polymerase on a MS2 RNA template. Our study sheds light on the template requirements of the MS2 replicase and paves the way for new in vitro applications including the design of genetic circuits combining both DNA- and RNA-encoded systems.

A Polymer Physics Framework for the Entropy of Arbitrary Pseudoknots

The accurate prediction of RNA secondary structure from primary sequence has had enormous impact on research from the past 40 years. Although many algorithms are available to make these predictions, the inclusion of non-nested loops, termed pseudoknots, still poses challenges arising from two main factors: 1) no physical model exists to estimate the loop entropies of complex intramolecular pseudoknots, and 2) their NP-complete enumeration has impeded their study. Here, we address both challenges. First, we develop a polymer physics model that can address arbitrarily complex pseudoknots using only two parameters corresponding to concrete physical quantities-over an order of magnitude fewer than the sparsest state-of-the-art phenomenological methods. Second, by coupling this model to exhaustive enumeration of the set of possible structures, we compute the entire free energy landscape of secondary structures resulting from a primary RNA sequence. We demonstrate that for RNA structures of ∼80 nucleotides, with minimal heuristics, the complete enumeration of possible secondary structures can be accomplished quickly despite the NP-complete nature of the problem. We further show that despite our loop entropy model's parametric sparsity, it performs better than or on par with previously published methods in predicting both pseudoknotted and non-pseudoknotted structures on a benchmark data set of RNA structures of ≤80 nucleotides. We suggest ways in which the accuracy of the model can be further improved.

Conformational Dynamics of the RNA G-Quadruplex and its Effect on Translation Efficiency

During translation, intracellular mRNA folds co-transcriptionally and must refold following the passage of ribosome. The mRNAs can be entrapped in metastable structures during these folding events. In the present study, we evaluated the conformational dynamics of the kinetically favored, metastable, and hairpin-like structure, which disturbs the thermodynamically favored G-quadruplex structure, and its effect on co-transcriptional translation in prokaryotic cells. We found that nascent mRNA forms a metastable hairpin-like structure during co-transcriptional folding instead of the G-quadruplex structure. When the translation progressed co-transcriptionally before the metastable hairpin-like structure transition to the G-quadruplex, function of the G-quadruplex as a roadblock of the ribosome was sequestered. This suggested that kinetically formed RNA structures had a dominant effect on gene expression in prokaryotes. The results of this study indicate that it is critical to consider the conformational dynamics of RNA-folding to understand the contributions of the mRNA structures in controlling gene expression.

Mimicking Co-Transcriptional RNA Folding Using a Superhelicase

Vectorial folding of RNA during transcription can produce intermediates with distinct biochemical activities. Here, we design an artificial minimal system to mimic cotranscriptional RNA folding in vitro. In this system, a presynthesized RNA molecule begins to fold from its 5'-end, as it is released from a heteroduplex by an engineered helicase that translocates on the complementary DNA strand in the 3'-to-5' direction. This chemically stabilized "superhelicase" Rep-X processively unwinds thousands of base pairs of DNA. The presynthesized RNA enables us to flexibly position fluorescent labels on the RNA for single-molecule fluorescence resonance energy transfer analysis and allows us to study real-time conformational dynamics during the vectorial folding process. We observed distinct signatures of the maiden secondary and tertiary folding of the Oryza sativa twister ribozyme. The maiden vectorial tertiary folding transitions occurred faster than Mg²⁺-induced refolding, but were also more prone to misfolding, likely due to sequential formation of alternative secondary structures. This novel assay can be applied to studying other kinetically controlled processes, such as riboswitch control and RNA-protein assembly.

RNA Phage Biology in a Metagenomic Era

The number of novel bacteriophage sequences has expanded significantly as a result of many metagenomic studies of phage populations in diverse environments. Most of these novel sequences bear little or no homology to existing databases (referred to as the "viral dark matter"). Also, these sequences are primarily derived from DNA-encoded bacteriophages (phages) with few RNA phages included. Despite the rapid advancements in high-throughput sequencing, few studies enrich for RNA viruses, i.e., target viral rather than cellular fraction and/or RNA rather than DNA via a reverse transcriptase step, in an attempt to capture the RNA viruses present in a microbial communities. It is timely to compile existing and relevant information about RNA phages to provide an insight into many of their important biological features, which should aid in sequence-based discovery and in their subsequent annotation. Without comprehensive studies, the biological significance of RNA phages has been largely ignored. Future bacteriophage studies should be adapted to ensure they are properly represented in phageomic studies.

Function of the RNA Coliphage Qβ Proteins in Medical In Vitro Evolution

Qβ is a positive (+) single-stranded RNA bacteriophage covered by a 25 nm icosahedral shell. Qβ belongs to the family of Leviviridae and is found throughout the world (bacterial isolates and sewage). The genome of Qβ is about 4.2 kb, coding for four proteins. This genome is surrounded by 180 copies of coat proteins (capsomers) each comprised of 132 residues of amino acids. The other proteins, the subunit II (β) of a replicase, the maturation protein (A₂) and the read-through or minor coat protein (A₁), play a key role in phage infection. With the replicase protein, which lacks proofreading activity, as well as its short replication time, and high population size, Qβ phage has attractive features for in vitro evolution. The A₁ protein gene shares the same initiation codon with the coat protein gene and is produced during translation when the coat protein’s UGA stop codon triplet (about 400 nucleotides from the initiation) is suppressed by a low level of ribosome misincorporation of tryptophan. Thus, A₁ is termed the read-through protein. This RNA phage platform technology not only serves to display foreign peptides but is also exceptionally suited to address questions about in vitro evolution. The C-terminus of A₁ protein confers to this RNA phage platform an exceptional feature of not only a linker for foreign peptide to be displayed also a model for evolution. This platform was used to present a peptide library of the G-H loop of the capsid region P1 of the foot-and-mouth disease virus (FMDV) called VP1 protein. The library was exposed on the exterior surface of Qβ phages, evolved and selected with the monoclonal antibodies (mAbs) SD6 of the FMDV. These hybrid phages could principally be good candidates for FMDV vaccine development. Separately, the membrane proximal external region (MPER) of human immunodeficiency virus type 1 (HIV-1) epitopes was fused with the A₁ proteins and exposed on the Qβ phage exterior surface. The engineered phages with MPER epitopes were recognized by anti-MPER specific antibodies. This system could be used to overcome the challenge of effective presentation of MPER to the immune system. A key portion of this linear epitope could be randomized and evolved with the Qβ system. Overall, antigens and epitopes of RNA viruses relevant to public health can be randomized, evolved and selected in pools using the proposed Qβ model to overcome their plasticity and the challenge of vaccine development. Major epitopes of a particular virus can be engineered or displayed on the Qβ phage surface and used for vaccine efficacy evaluation, thus avoiding the use of live viruses.

Dual function of a cis-acting RNA element that acts as a replication enhancer and a translation repressor in a plant positive-stranded RNA virus

The genome of red clover necrotic mosaic virus is divided into two positive-stranded RNA molecules of RNA1 and RNA2, which have no 5' cap structure and no 3' poly(A) tail. Previously, we showed that any mutations in the cis-acting RNA replication elements of RNA2 abolished its cap-independent translational activity, suggesting a strong link between RNA replication and translation. Here, we investigated the functions of the 5' untranslated region (UTR) of RNA2 and revealed that the basal stem-structure (5'BS) predicted in the 5' UTR is essential for robust RNA replication. Interestingly, RNA2 mutants with substitution or deletion in the right side of the 5'BS showed strong translational activity, despite their impaired replication competency. Furthermore, nucleotide sequences other than the 5'BS of the 5' UTR were essential to facilitate the replication-associated translation. Overall, these cis-acting RNA elements seem to coordinately regulate the balance between RNA replication and replication-associated translation.

Structures of Qβ virions, virus-like particles, and the Qβ-MurA complex reveal internal coat proteins and the mechanism of host lysis

In single-stranded RNA bacteriophages (ssRNA phages) a single copy of the maturation protein binds the genomic RNA (gRNA) and is required for attachment of the phage to the host pilus. For the canonical Allolevivirus Qβ the maturation protein, A₂, has an additional role as the lysis protein, by its ability to bind and inhibit MurA, which is involved in peptidoglycan biosynthesis. Here, we determined structures of Qβ virions, virus-like particles, and the Qβ-MurA complex using single-particle cryoelectron microscopy, at 4.7-Å, 3.3-Å, and 6.1-Å resolutions, respectively. We identified the outer surface of the β-region in A₂ as the MurA-binding interface. Moreover, the pattern of MurA mutations that block Qβ lysis and the conformational changes of MurA that facilitate A₂ binding were found to be due to the intimate fit between A₂ and the region encompassing the closed catalytic cleft of substrate-liganded MurA. Additionally, by comparing the Qβ virion with Qβ virus-like particles that lack a maturation protein, we observed a structural rearrangement in the capsid coat proteins that is required to package the viral gRNA in its dominant conformation. Unexpectedly, we found a coat protein dimer sequestered in the interior of the virion. This coat protein dimer binds to the gRNA and interacts with the buried α-region of A₂, suggesting that it is sequestered during the early stage of capsid formation to promote the gRNA condensation required for genome packaging. These internalized coat proteins are the most asymmetrically arranged major capsid proteins yet observed in virus structures.

In situ structures of the genome and genome-delivery apparatus in a single-stranded RNA virus

Genome packaging into a protein capsid and its subsequent delivery into a host cell are two fundamental processes in the life cycle of a virus. Unlike dsDNA viruses which pump their genome into a preformed capsid^1-3, ssRNA viruses, such as bacteriophage MS2, co-assemble their capsid with genome^4-7; however, the structural basis of this co-assembly is poorly understood. MS2 infects Escherichia coli via host “sex” pilus (F-pilus)⁸ and is the first fully-sequenced organism⁹ and a model system for studies of gene translational regulations^10,11, RNA-protein interactions^12-14, and RNA virus assembly^15-17. Its positive-sense ssRNA genome of 3569 bases is enclosed in a capsid with one maturation protein (MP) monomer and 89 coat protein (CP) dimers arranged in a T=3 icosahedral lattice^18,19. MP is responsible for attaching the virus to an F-pilus and delivering the viral genome into the host during infection⁸, but how the genome is organized and delivered are not known. Here we show the MS2 structure at 3.6Å resolution determined by electron-counting cryo electron microscopy (cryoEM) and asymmetric reconstruction. We traced ~80% backbone of the viral genome, built atomic models for 16 RNA stem-loops, and identified three conserved motifs of RNA-CP interactions among 15 of these stem-loops with diverse sequences. The stem-loop at 3’ end of the genome interacts extensively with the MP, which, with just a six-helix bundle and a six-stranded β-sheet, forms a genome-delivery apparatus, and joins 89 CP-dimers to form a capsid. This first atomic description of genome-capsid interactions in a spherical ssRNA virus provides insights into genome delivery via host “sex” pilus and mechanisms underlying ssRNA-capsid co-assembly, and inspires imaginations about links between nucleoprotein complexes and the origin of viruses.

Asymmetric cryo-EM reconstruction of phage MS2 reveals genome structure in situ

In single-stranded ribonucleic acid (RNA) viruses, virus capsid assembly and genome packaging are intertwined processes. Using cryo-electron microscopy and single particle analysis we determined the asymmetric virion structure of bacteriophage MS2, which includes 178 copies of the coat protein, a single copy of the A-protein and the RNA genome. This reveals that in situ, the viral RNA genome can adopt a defined conformation. The RNA forms a branched network of stem-loops that almost all allocate near the capsid inner surface, while predominantly binding to coat protein dimers that are located in one-half of the capsid. This suggests that genomic RNA is highly involved in genome packaging and virion assembly.

MS2 is a single-stranded RNA bacteriophage that infects its host via adsorption to bacterial pili. Here the authors visualize the MS2 virion with asymmetric cryo-EM reconstruction, revealing that the genome of MS2 adopts a specific structure of asymmetrically distributed stem-loops connected to the capsid.

Coordination of Genomic RNA Packaging with Viral Assembly in HIV-1

The tremendous progress made in unraveling the complexities of human immunodeficiency virus (HIV) replication has resulted in a library of drugs to target key aspects of the replication cycle of the virus. Yet, despite this accumulated wealth of knowledge, we still have much to learn about certain viral processes. One of these is virus assembly, where the viral genome and proteins come together to form infectious progeny. Here we review this topic from the perspective of how the route to production of an infectious virion is orchestrated by the viral genome, and we compare and contrast aspects of the assembly mechanisms employed by HIV-1 with those of other RNA viruses.

Translation Initiation is Controlled by RNA Folding Kinetics via a Ribosome Drafting Mechanism

RNA folding plays an important role in controlling protein synthesis as well as other cellular processes. Existing models have focused on how RNA folding energetics control translation initiation rate under equilibrium conditions but have largely ignored the effects of nonequilibrium RNA folding. We introduce a new mechanism, called "ribosome drafting", that explains how a mRNA's folding kinetics and the ribosome's binding rate collectively control its translation initiation rate. During cycles of translation, ribosome drafting emerges whenever successive ribosomes bind to a mRNA faster than the mRNA can refold, maintaining it in a nonequilibrium state with an acceleration of protein synthesis. Using computational design, time-correlated single photon counting, and expression measurements, we demonstrate that slow-folding and fast-folding RNA structures with equivalent folding energetics can vary protein synthesis rates by 1000-fold. We determine the necessary conditions for ribosome drafting by characterizing mRNAs with rationally designed ribosome binding rates, folding kinetics, and folding energetics, confirming the predictions of a nonequilibrium Markov model of translation. Our results have widespread implications, illustrating how competitive folding and assembly kinetics can shape the gene expression machinery's sequence-structure-function relationship inside cells.

The Shine-Dalgarno sequence of riboswitch-regulated single mRNAs shows ligand-dependent accessibility bursts

In response to intracellular signals in Gram-negative bacteria, translational riboswitches—commonly embedded in messenger RNAs (mRNAs)—regulate gene expression through inhibition of translation initiation. It is generally thought that this regulation originates from occlusion of the Shine-Dalgarno (SD) sequence upon ligand binding; however, little direct evidence exists. Here we develop Single Molecule Kinetic Analysis of RNA Transient Structure (SiM-KARTS) to investigate the ligand-dependent accessibility of the SD sequence of an mRNA hosting the 7-aminomethyl-7-deazaguanine (preQ₁)-sensing riboswitch. Spike train analysis reveals that individual mRNA molecules alternate between two conformational states, distinguished by ‘bursts' of probe binding associated with increased SD sequence accessibility. Addition of preQ₁ decreases the lifetime of the SD's high-accessibility (bursting) state and prolongs the time between bursts. In addition, ligand-jump experiments reveal imperfect riboswitching of single mRNA molecules. Such complex ligand sensing by individual mRNA molecules rationalizes the nuanced ligand response observed during bulk mRNA translation.

In response to intracellular signals, bacterial translational riboswitches embedded in mRNAs can regulate gene expression through inhibition of translation initiation. Here, the authors describe SiM-KARTS, a novel approach for detecting changes in the structure of single RNA molecules in response to a ligand.

Four RNA families with functional transient structures

Protein-coding and non-coding RNA transcripts perform a wide variety of cellular functions in diverse organisms. Several of their functional roles are expressed and modulated via RNA structure. A given transcript, however, can have more than a single functional RNA structure throughout its life, a fact which has been previously overlooked. Transient RNA structures, for example, are only present during specific time intervals and cellular conditions. We here introduce four RNA families with transient RNA structures that play distinct and diverse functional roles. Moreover, we show that these transient RNA structures are structurally well-defined and evolutionarily conserved. Since Rfam annotates one structure for each family, there is either no annotation for these transient structures or no such family. Thus, our alignments either significantly update and extend the existing Rfam families or introduce a new RNA family to Rfam. For each of the four RNA families, we compile a multiple-sequence alignment based on experimentally verified transient and dominant (dominant in terms of either the thermodynamic stability and/or attention received so far) RNA secondary structures using a combination of automated search via covariance model and manual curation. The first alignment is the Trp operon leader which regulates the operon transcription in response to tryptophan abundance through alternative structures. The second alignment is the HDV ribozyme which we extend to the 5' flanking sequence. This flanking sequence is involved in the regulation of the transcript's self-cleavage activity. The third alignment is the 5' UTR of the maturation protein from Levivirus which contains a transient structure that temporarily postpones the formation of the final inhibitory structure to allow translation of maturation protein. The fourth and last alignment is the SAM riboswitch which regulates the downstream gene expression by assuming alternative structures upon binding of SAM. All transient and dominant structures are mapped to our new alignments introduced here.

Highly efficient in vitro and in vivo delivery of functional RNAs using new versatile MS2-chimeric retrovirus-like particles

RNA delivery is an attractive strategy to achieve transient gene expression in research projects and in cell- or gene-based therapies. Despite significant efforts investigating vector-directed RNA transfer, there is still a requirement for better efficiency of delivery to primary cells and in vivo. Retroviral platforms drive RNA delivery, yet retrovirus RNA-packaging constraints limit gene transfer to two genome-molecules per viral particle. To improve retroviral transfer, we designed a dimerization-independent MS2-driven RNA packaging system using MS2-Coat-retrovirus chimeras. The engineered chimeric particles promoted effective packaging of several types of RNAs and enabled efficient transfer of biologically active RNAs in various cell types, including human CD34(+) and iPS cells. Systemic injection of high-titer particles led to gene expression in mouse liver and transferring Cre-recombinase mRNA in muscle permitted widespread editing at the ROSA26 locus. We could further show that the VLPs were able to activate an osteoblast differentiation pathway by delivering RUNX2- or DLX5-mRNA into primary human bone-marrow mesenchymal-stem cells. Thus, the novel chimeric MS2-lentiviral particles are a versatile tool for a wide range of applications including cellular-programming or genome-editing.

Interconversion between parallel and antiparallel conformations of a 4H RNA junction in domain 3 of foot-and-mouth disease virus IRES captured by dynamics simulations

RNA junctions are common secondary structural elements present in a wide range of RNA species. They play crucial roles in directing the overall folding of RNA molecules as well as in a variety of biological functions. In particular, there has been great interest in the dynamics of RNA junctions, including conformational pathways of fully base-paired 4-way (4H) RNA junctions. In such constructs, all nucleotides participate in one of the four double-stranded stem regions, with no connecting loops. Dynamical aspects of these 4H RNAs are interesting because frequent interchanges between parallel and antiparallel conformations are thought to occur without binding of other factors. Gel electrophoresis and single-molecule fluorescence resonance energy transfer experiments have suggested two possible pathways: one involves a helical rearrangement via disruption of coaxial stacking, and the other occurs by a rotation between the helical axes of coaxially stacked conformers. Employing molecular dynamics simulations, we explore this conformational variability in a 4H junction derived from domain 3 of the foot-and-mouth disease virus internal ribosome entry site (IRES); this junction contains highly conserved motifs for RNA-RNA and RNA-protein interactions, important for IRES activity. Our simulations capture transitions of the 4H junction between parallel and antiparallel conformations. The interconversion is virtually barrier-free and occurs via a rotation between the axes of coaxially stacked helices with a transient perpendicular intermediate. We characterize this transition, with various interhelical orientations, by pseudodihedral angle and interhelical distance measures. The high flexibility of the junction, as also demonstrated experimentally, is suitable for IRES activity. Because foot-and-mouth disease virus IRES structure depends on long-range interactions involving domain 3, the perpendicular intermediate, which maintains coaxial stacking of helices and thereby consensus primary and secondary structure information, may be beneficial for guiding the overall organization of the RNA system in domain 3.

Dynamic energy landscapes of riboswitches help interpret conformational rearrangements and function

Riboswitches are RNAs that modulate gene expression by ligand-induced conformational changes. However, the way in which sequence dictates alternative folding pathways of gene regulation remains unclear. In this study, we compute energy landscapes, which describe the accessible secondary structures for a range of sequence lengths, to analyze the transcriptional process as a given sequence elongates to full length. In line with experimental evidence, we find that most riboswitch landscapes can be characterized by three broad classes as a function of sequence length in terms of the distribution and barrier type of the conformational clusters: low-barrier landscape with an ensemble of different conformations in equilibrium before encountering a substrate; barrier-free landscape in which a direct, dominant “downhill” pathway to the minimum free energy structure is apparent; and a barrier-dominated landscape with two isolated conformational states, each associated with a different biological function. Sharing concepts with the “new view” of protein folding energy landscapes, we term the three sequence ranges above as the sensing, downhill folding, and functional windows, respectively. We find that these energy landscape patterns are conserved in various riboswitch classes, though the order of the windows may vary. In fact, the order of the three windows suggests either kinetic or thermodynamic control of ligand binding. These findings help understand riboswitch structure/function relationships and open new avenues to riboswitch design.

Riboswitches are RNAs that modulate gene expression by ligand-induced conformational changes. However, the way that sequence dictates alternative folding pathways of gene regulation remains unclear. In this study, we mimic transcription by computing energy landscapes which describe accessible secondary structures for a range of sequence lengths. Consistent with experimental evidence, we find that most riboswitch landscapes can be characterized by three broad classes as a function of sequence length in terms of the distribution and barrier type of the conformational clusters: Low-barrier landscape with an ensemble of conformations in equilibrium before encountering a substrate; barrier-free landscape with a dominant “downhill” pathway to the minimum free energy structure; and barrier-dominated landscape with two isolated conformational states with different functions. Sharing concepts with the “new view” of protein folding energy landscapes, we term the three sequence ranges above as the sensing, downhill folding, and functional windows, respectively. We find that these energy landscape patterns are conserved between riboswitch classes, though the order of the windows may vary. In fact, the order of the three windows suggests either kinetic or thermodynamic control of ligand binding. These findings help understand riboswitch structure/function relationships and open new avenues to riboswitch design.

The jerky and knotty dynamics of RNA

RNA is known to exhibit a jerky dynamics, as intramolecular thermal motion, on <0.1 micros time scales, is punctuated by infrequent structural rearrangements on much longer time scales, i.e. from >10 micros up to a few minutes or even hours. These rare stochastic events correspond to the formation or dissociation of entire stems through cooperative base pairing/unpairing transitions. Such a clear separation of time scales in RNA dynamics has made it possible to implement coarse grained RNA simulations, which predict RNA folding and unfolding pathways including kinetically trapped structures on biologically relevant time scales of seconds to minutes. RNA folding simulations also enable to predict the formation of pseudoknots, that is, helices interior to loops, which mechanically restrain the relative orientations of other non-nested helices. But beyond static structural constraints, pseudoknots can also strongly affect the folding and unfolding dynamics of RNA, as the order by which successive helices are formed and dissociated can lead to topologically blocked transition intermediates. The resulting knotty dynamics can enhance the stability of RNA switches, improve the efficacy of co-transcriptional folding pathways and lead to unusual self-assembly properties of RNA.

Killer and protective ribosomes

In prokaryotes, translation influences mRNA decay. The breakdown of most Escherichia coli mRNAs is initiated by RNase E, a 5'-dependent endonuclease. Some mRNAs are protected by ribosomes even if these are located far upstream of cleavage sites ("protection at a distance"), whereas others require direct shielding of these sites. I argue that these situations reflect different modes of interaction of RNase E with mRNAs. Protection at a distance is most impressive in Bacilli, where ribosomes can protect kilobases of unstable downstream sequences. I propose that this protection reflects the role in mRNA decay of RNase J1, a 5'-->3' exonuclease with no E. coli equivalent. Finally, recent years have shown that besides their protective role, ribosomes can also cleave their mRNA under circumstances that cause ribosome stalling. The endonuclease associated with this "killing" activity, which has a eukaryotic counterpart ("no-go decay"), is not characterized; it may be borne by the distressed ribosome itself.

Design of simple synthetic RNA thermometers for temperature-controlled gene expression in Escherichia coli

RNA thermometers are thermosensors that regulate gene expression by temperature-induced changes in RNA conformation. Naturally occurring RNA thermometers exhibit complex secondary structures which are believed to undergo a series of gradual structural changes in response to temperature shifts. Here, we report the de novo design of considerably simpler RNA thermometers that provide useful RNA-only tools to regulate bacterial gene expression by a shift in the growth temperature. We show that a single small stem-loop structure containing the ribosome binding site is sufficient to construct synthetic RNA thermometers that work efficiently at physiological temperatures. Our data suggest that the thermometers function by a simple melting mechanism and thus provide minimum size on/off switches to experimentally induce or repress gene expression by temperature.

Folding kinetics of large RNAs

We introduce here a heuristic approach to kinetic RNA folding that constructs secondary structures by stepwise combination of building blocks. These blocks correspond to subsequences and their thermodynamically optimal structures. These are determined by the standard dynamic programming approach to RNA folding. Folding trajectories are modeled at base-pair resolution using the Morgan-Higgs heuristic and a barrier tree-based heuristic to connect combinations of the local building blocks. Implemented in the program Kinwalker, the algorithm allows co-transcriptional folding and can be used to fold sequences of up to about 1500 nucleotides in length. A detailed comparison with several well-studied examples from the literature, including the delayed folding of bacteriophage cloverleaf structures, the adenine sensing riboswitch, and the hok RNA, shows an excellent agreement of predicted trajectories and experimental evidence. The software is available as part of the ViennaRNA Package.

A model for the origin of protein synthesis as coreplicational scanning of nascent RNA

The origin of protein synthesis is one of the major riddles of molecular biology. It was proposed a decade ago that the ribosomal RNA evolved from an earlier RNA-replisome (a ribozyme fulfilling RNA replication) while transfer RNA (tRNA) evolved from a genomic replication origin. Applying these hypotheses, I suggest that protein synthesis arose for the purpose of segregating copy and template RNA during replication through the conventional formation of a complementary strand. Nascent RNA was scanned in 5' to 3' direction following the progress of replication. The base pairing of several tRNA-like molecules with nascent RNA released the replication intermediates trapped in duplex. Synthesis of random peptides evolved to fuel the turnover of tRNAs. Then the combination of replication-coupled peptide formation and the independent development of amino acid-specific tRNA aminoacylation resulted in template-based protein synthesis. Therefore, the positioning of tRNAs adjacent to each other developed for the purpose of replication rather than peptide synthesis. This hypothesis does not include either selection for useful peptides or specific recognition of amino acids at the initial evolution of translation. It does, however, explain a number of features of modern translation apparatus, such as the relative flexibility of genetic code, the number of proteins shared by the transcription and translation machines, the universal participation of an RNA subunit in co-translational protein secretion, 'unscheduled translation', and factor-independent translocation. Assistance of original ribosomes in keeping apart the nascent transcript from its template is still widely explored by modern bacteria and perhaps by other domains of life.

Encoding folding paths of RNA switches

RNA co-transcriptional folding has long been suspected to play an active role in helping proper native folding of ribozymes and structured regulatory motifs in mRNA untranslated regions (UTRs). Yet, the underlying mechanisms and coding requirements for efficient co-transcriptional folding remain unclear. Traditional approaches have intrinsic limitations to dissect RNA folding paths, as they rely on sequence mutations or circular permutations that typically perturb both RNA folding paths and equilibrium structures. Here, we show that exploiting sequence symmetries instead of mutations can circumvent this problem by essentially decoupling folding paths from equilibrium structures of designed RNA sequences. Using bistable RNA switches with symmetrical helices conserved under sequence reversal, we demonstrate experimentally that native and transiently formed helices can guide efficient co-transcriptional folding into either long-lived structure of these RNA switches. Their folding path is controlled by the order of helix nucleations and subsequent exchanges during transcription, and may also be redirected by transient antisense interactions. Hence, transient intra- and inter-molecular base pair interactions can effectively regulate the folding of nascent RNA molecules into different native structures, provided limited coding requirements, as discussed from an information theory perspective. This constitutive coupling between RNA synthesis and RNA folding regulation may have enabled the early emergence of autonomous RNA-based regulation networks.

Self-splicing of a group I intron reveals partitioning of native and misfolded RNA populations in yeast

Stable RNAs must form specific three-dimensional structures, yet many RNAs become kinetically trapped in misfolded conformations. To understand the factors that control the accuracy of RNA folding in the cell, the self-splicing activity of the Tetrahymena group I intron was compared in different genetic contexts in budding yeast. The extent of splicing was 98% when the intron was placed in its natural rDNA context, but only 3% when the intron was expressed in an exogenous pre-mRNA. Further experiments showed that the probability of forming the active intron structure depends on local sequence context and transcription by Pol I. Pre-rRNAs decayed at similar rates, whether the intron was wild type or inactivated by an internal deletion, suggesting that most of the unreacted pre-rRNA is incompetent to splice. Northern blots and complementation assays showed that mutations that destabilize the intron tertiary structure inhibited self-splicing and processing of internal transcribed spacer 2. The data are consistent with partitioning of pre-rRNAs into active and inactive populations. The misfolded RNAs are sequestered and degraded without refolding to a significant extent. Thus, the initial fidelity of folding can dictate the intracellular fate of transcripts containing this group I intron.

Computer simulation of chaperone effects of Archaeal C/D box sRNA binding on rRNA folding

Archaeal C/D box small RNAs (sRNAs) are homologues of eukaryotic C/D box small nucleolar RNAs (snoRNAs). Their main function is guiding 2'-O-ribose methylation of nucleotides in rRNAs. The methylation requires the pairing of an sRNA antisense element to an rRNA target site with formation of an RNA-RNA duplex. The temporary formation of such a duplex during rRNA maturation is expected to influence rRNA folding in a chaperone-like way, in particular in thermophilic Archaea, where multiple sRNAs with two binding sites are found. Here we investigate possible mechanisms of chaperone function of Archaeoglobus fulgidus and Pyrococcus abyssi C/D box sRNAs using computer simulations of rRNA secondary structure formation by genetic algorithm. The effects of sRNA binding on rRNA structure are introduced as temporary structural constraints during co-transcriptional folding. Comparisons of the final predictions with simulations without sRNA binding and with phylogenetic structures show that sRNAs with two antisense elements may significantly facilitate the correct formation of long-range interactions in rRNAs, in particular at elevated temperatures. The simulations suggest that the main mechanism of this effect is a transient restriction of folding in rRNA domains where the termini are brought together by binding to double-guide sRNAs.

Dynamics of intramolecular recognition: base-pairing in DNA/RNA near and far from equilibrium

The physics of the base-pairing interaction in DNA and RNA molecules plays a fundamental role in biology. Past experimental and theoretical research has led to a fairly complete and quantitative understanding of the equilibrium properties such as the different phases, the melting behavior, and the response to slow stretching. The non-equilibrium behavior is even richer than might be expected on the basis of thermodynamics. However, the non-equilibrium behavior is also far less understood. Here, we review different theoretical approaches to the study of base-pairing thermodynamics and kinetics, and illustrate the rich phenomenology with several examples that use these approaches.

Dynamics-based amplification of RNA function and its characterization by using NMR spectroscopy

The ever-increasing cellular roles ascribed to RNA raise fundamental questions regarding how a biopolymer composed of only four chemically similar building-block nucleotides achieves such functional diversity. Here, I discuss how RNA achieves added mechanistic and chemical complexity by undergoing highly controlled conformational changes in response to a variety of cellular signals. I examine pathways for achieving selectivity in these conformational changes that rely to different extents on the structure and dynamics of RNA. Finally, I review solution-state NMR techniques that can be used to characterize RNA structural dynamics and its relationship to function.

Regulation of translation via mRNA structure in prokaryotes and eukaryotes

The mechanism of initiation of translation differs between prokaryotes and eukaryotes, and the strategies used for regulation differ accordingly. Translation in prokaryotes is usually regulated by blocking access to the initiation site. This is accomplished via base-paired structures (within the mRNA itself, or between the mRNA and a small trans-acting RNA) or via mRNA-binding proteins. Classic examples of each mechanism are described. The polycistronic structure of mRNAs is an important aspect of translational control in prokaryotes, but polycistronic mRNAs are not usable (and usually not produced) in eukaryotes. Four structural elements in eukaryotic mRNAs are important for regulating translation: (i) the m7G cap; (ii) sequences flanking the AUG start codon; (iii) the position of the AUG codon relative to the 5' end of the mRNA; and (iv) secondary structure within the mRNA leader sequence. The scanning model provides a framework for understanding these effects. The scanning mechanism also explains how small open reading frames near the 5' end of the mRNA can down-regulate translation. This constraint is sometimes abrogated by changing the structure of the mRNA, sometimes with clinical consequences. Examples are described. Some mistaken ideas about regulation of translation that have found their way into textbooks are pointed out and corrected.

Initiation of protein synthesis in bacteria

Valuable information on translation initiation is available from biochemical data and recently solved structures. We present a detailed description of current knowledge about the structure, function, and interactions of the individual components involved in bacterial translation initiation. The first section describes the ribosomal features relevant to the initiation process. Subsequent sections describe the structure, function, and interactions of the mRNA, the initiator tRNA, and the initiation factors IF1, IF2, and IF3. Finally, we provide an overview of mechanisms of regulation of the translation initiation event. Translation occurs on ribonucleoprotein complexes called ribosomes. The ribosome is composed of a large subunit and a small subunit that hold the activities of peptidyltransfer and decode the triplet code of the mRNA, respectively. Translation initiation is promoted by IF1, IF2, and IF3, which mediate base pairing of the initiator tRNA anticodon to the mRNA initiation codon located in the ribosomal P-site. The mechanism of translation initiation differs for canonical and leaderless mRNAs, since the latter is dependent on the relative level of the initiation factors. Regulation of translation occurs primarily in the initiation phase. Secondary structures at the mRNA ribosomal binding site (RBS) inhibit translation initiation. The accessibility of the RBS is regulated by temperature and binding of small metabolites, proteins, or antisense RNAs. The future challenge is to obtain atomic-resolution structures of complete initiation complexes in order to understand the mechanism of translation initiation in molecular detail.

cis-acting RNA signals in the NS5B C-terminal coding sequence of the hepatitis C virus genome

The cis-replicating RNA elements in the 5' and 3' nontranslated regions (NTRs) of the hepatitis C virus (HCV) genome have been thoroughly studied before. However, no cis-replicating elements have been identified in the coding sequences of the HCV polyprotein until very recently. The existence of highly conserved and stable stem-loop structures in the RNA polymerase NS5B coding sequence, however, has been previously predicted (A. Tuplin, J. Wood, D. J. Evans, A. H. Patel, and P. Simmonds, RNA 8:824-841, 2002). We have selected for our studies a 249-nt-long RNA segment in the C-terminal NS5B coding region (NS5BCR), which is predicted to form four stable stem-loop structures (SL-IV to SL-VII). By deletion and mutational analyses of the RNA structures, we have determined that two of the stem-loops (SL-V and SL-VI) are essential for replication of the HCV subgenomic replicon in Huh-7 cells. Mutations in the loop and the top of the stem of these RNA elements abolished replicon RNA synthesis but had no effect on translation. In vitro gel shift and filter-binding assays revealed that purified NS5B specifically binds to SL-V. The NS5B-RNA complexes were specifically competed away by unlabeled homologous RNA, to a small extent by 3' NTR RNA, and only poorly by 5' NTR RNA. The other two stem-loops (SL-IV and SL-VII) of the NS5BCR domain were found to be important but not essential for colony formation by the subgenomic replicon. The precise function(s) of these cis-acting RNA elements is not known.

Translational repression mechanisms in prokaryotes

Translational repression results from a complex choreography of macromolecular interactions interfering with the formation of translational initiation complexes. The relationship between the rate and extent of formation of these interactions to form repressed mRNA complexes determines the extent of repression. A novel analysis of repression mechanisms is presented here and it indicates that the reversibility of repressed complex formation influences the steady state balance of the distribution of translationally active and inactive complexes and therefore has an impact on the efficiency of repression. Reviewed here is evidence for three distinct translational repression mechanisms, regulating expression of the transcription factor sigma32, threonine tRNA synthetase and ribosomal proteins on the alpha operon in Escherichia coli. Efficient regulation of expression in these systems makes use of specific mRNA structures in quite different ways.

Effect of transcription on folding of the Tetrahymena ribozyme

Sequential formation of RNA interactions during transcription can bias the folding pathway and ultimately determine the functional state of a transcript. The kinetics of cotranscriptional folding of the Tetrahymena L-21 ribozyme was compared with refolding of full-length transcripts under the same conditions. Sequential folding after transcription by phage T7 or Escherichia coli polymerase is only twice as fast as refolding, and the yield of native RNA is the same. By contrast, a greater fraction of circularly permuted variants folded correctly at early times during transcription than during refolding. Hybridization of complementary oligonucleotides suggests that cotranscriptional folding enables a permuted RNA beginning at G303 to escape non-native interactions in P3 and P9. We propose that base pairing of upstream sequences during transcription elongation favors branched secondary structures that increase the probability of forming the native ribozyme structure.

Bistable secondary structures of small RNAs and their structural probing by comparative imino proton NMR spectroscopy

We investigate 25-34 nucleotide RNA sequences, that have been rationally designed to adopt two different secondary structures that are in thermodynamic equilibrium. Experimental evidence for the co-existence of the two conformers results from the NH...N 1H NMR spectra. When compared to the NH...N 1H NMR spectra of appropriate reference sequences the equilibrium position is easily quantifiable even without the assignment of the individual NH resonances. The reference sequences represent several Watson-Crick base-paired double helical segments, each encountered in either of the two conformers of the bistable target sequence. In addition, we rationalize the influence of nucleotide mutations on the equilibrium position of one of the bistable RNA sequences. The approach further allows a detailed thermodynamic analysis and the evaluation of secondary structure predictions for multistable RNAs obtained by computational methods.

Self-induced structural switches in RNA

Many biologically active RNAs show a switch in their secondary structure, which is accompanied by changes in their function. Such changes in secondary structure often require trans-acting factors, e.g. RNA chaperones. However, several biologically active RNAs do not require trans-acting factors for this structural switch, which is therefore indicated here as a "self-induced switch". These self-induced structural switches have several characteristics in common. They all start from a metastable structure, which is maintained for some time allowing or blocking a particular function of the RNA. Hereafter, a structural element becomes available, e.g. during transcription, triggering a rapid transition into a stable conformation, which again is accompanied by either a gain or loss of function. A further common element of this type of switches is the involvement of a branch migration or strand displacement reaction, which lowers the energy barrier of the reaction sufficiently to allow rapid refolding. Here, we review a number of these self-induced switches in RNA secondary structure as proposed for several systems. A general model for this type of switches is presented, showing its importance in the biology of functionally active RNAs.

A pH-jump approach for investigating secondary structure refolding kinetics in RNA

It has been shown that premature translation of the plasmid-mediated toxin in hok/sok of plasmid R1 and pnd/pndB of plasmid R483 is prevented during transcription of the hok and pnd mRNAs by the formation of metastable hairpins at the 5'-end of the mRNA. Here, an experimental approach is presented, which allows the accurate measurement of the refolding kinetics of the 5'-end RNA fragments in vitro without chemically modifying the RNA. The method is based on acid denaturation followed by a pH-jump to neutral pH as a novel way to trap kinetically favoured RNA secondary structures, allowing the measurement of a wide range of biologically relevant refolding rates, with or without the use of standard stopped-flow equipment. The refolding rates from the metastable to the stable conformation in both the hok74 and pnd58 5'-end RNA fragments were determined by using UV absorbance changes corresponding to the structural rearrangements. The measured energy barriers showed that the refolding path does not need complete unfolding of the metastable structures before the formation of the final structures. Two alternative models of such a pathway are discussed.

In vivo generation of hybrids between different species of RNA phages

Hybrids between different species or genera of the single-stranded RNA coliphages have not been found in nature. Here, it has been shown that viable hybrids between different phage species can easily be generated in the laboratory by in vivo recombination. cDNA of species I phage MS2 located on a plasmid and lacking part of its 5' untranslated leader (5' UTR) was complemented with another plasmid carrying the 5' half of the genome of fr, a species I phage, or of KU1, a species II representative with low sequence similarity. When the two plasmids were present in the same cell there was spontaneous production of hybrid phages. Interestingly, these hybrids did not arise by a double or single crossover that would replace the missing MS2 sequences with those of fr or KU1. Rather, hybrids arose by attaching the complete 5' UTR of fr or KU1 to the 5' terminus of the defective MS2 phage. Several elements of the 5' UTR then occurred twice, one from KU1 (or fr) and the other from MS2. These redundant elements are in most cases deleted upon evolution of the hybrids. As a result, the 5' UTR of KU1 (or fr) then replaced that of MS2. It was earlier shown that this 5' UTR could assume two alternating structures that facilitated transient translation of the proximal maturation gene. Apparently, this timer function of the 5' UTR was exchangeable and could function independently of the rest of the genome. When hybrids were competed against wild-type, they were quickly outgrown, probably explaining their absence from natural isolates.

Hairpin ribozymes with four-way helical junctions mediate intracellular RNA ligation

Virtually all RNA-mediated reactions require transitions among alternative RNA conformations. The complexity of biological reactions can obscure specific conformational changes in vivo and important features of the intracellular environment are difficult to reproduce in vitro. However, simple RNA self-cleavage and ligation reactions offer a unique opportunity to measure the kinetics and equilibria of specific RNA conformational transitions directly in living cells. Hairpin ribozymes that incorporate the natural four-way helical junction self-cleave rapidly in vivo, but only when cleavage products dissociate rapidly. Cleavage rates fall when cleavage products remain bound in stable base-paired helices, providing evidence that bound products undergo re-ligation. These results provide the first detailed kinetic description of an intracellular ribozyme reaction that includes cleavage, ligation and product dissociation rates. Kinetic and equilibrium parameters measured in vivo correspond well, but not perfectly, with values measured for the same reactions in vitro under conditions that approximate an intracellular ionic environment.

Translational control by delayed RNA folding: identification of the kinetic trap

The maturation or A-protein gene of single-stranded RNA phage MS2 is preceded by a 130-nt long untranslated leader. When MS2 RNA folding is at equilibrium, the gene is untranslatable because the leader adopts a well-defined cloverleaf structure in which the Shine-Dalgarno (SD) sequence of the maturation gene is taken up in long-distance base pairing with an upstream complementary sequence (UCS). Synthesis of the A-protein takes place transiently while the RNA is synthesized from the minus strand. This requires that formation of the inhibitory cloverleaf is slow. In vitro, the folding delay was on the order of minutes. Here, we present evidence that this postponed folding is caused by the formation of a metastable intermediate. This intermediate is a small local hairpin that contains the UCS in its loop, thereby preventing or slowing down its pairing with the SD sequence. Mutants in which the small hairpin could not be formed made no detectable amounts of A-protein and were barely viable. Apparently, here the cloverleaf formed quicker than ribosomes could bind. On the other hand, mutants in which the small intermediary hairpin was stabilized produced more A-protein than wild type and were viable. One hardly growing mutant that could not form the metastable hairpin and did not make detectable amounts of A-protein was evolved. The emerging pseudo-revertant had selected two second site repressor mutations that allowed reconstruction of a variant of the metastable intermediate. The pseudo-revertant had also regained the capacity to produce the A-protein.

Folding of branched RNA species

The global structures of branched RNA species are important to their function. Branched RNA species are defined as molecules in which double-helical segments are interrupted by abrupt discontinuities. These include helical junctions of different orders, and base bulges and loops. Common helical junctions are three- and four-way junctions, often interrupted by mispairs or additional nucleotides. There are many interesting examples of functional RNA junctions, including the hammerhead and hairpin ribozymes, and junctions that serve as binding sites for proteins. The junctions display some common structural properties. These include a tendency to undergo pairwise helical stacking and ion-induced conformational transitions. Helical branchpoints can act as key architectural components and as important sites for interactions with proteins. Copyright 1999 John Wiley & Sons, Inc.

Coupled nucleotide covariations reveal dynamic RNA interaction patterns

Evolutionarily conserved structures in related RNA molecules contain coordinated variations (covariations) of paired nucleotides. Analysis of covariations is a very powerful approach to deduce phylogenetically conserved (i.e., functional) conformations, including tertiary interactions. Here we discuss conserved RNA folding pathways that are revealed by covariation patterns. In such pathways, structural requirements for alternative pairings cause some nucleotides to covary with two different partners. Such "coupled" covariations between three or more nucleotides were found in various types of RNAs. The analysis of coupled covariations can unravel important features of RNA folding dynamics and improve phylogeny reconstruction in some cases. Importantly, it is necessary to distinguish between multiple covariations determined by mutually exclusive structures and those determined by tertiary contacts.

Alternative conformations of a nucleic acid four-way junction

A crystal structure of a 108 nucleotide RNA-DNA complex containing a four-way junction was solved at 3.1 A resolution. The structure of the junction differs substantially from the "stacked-X" conformation observed previously, due to a 135 degrees rotation of the branches. Comparison of the two conformers provides insight into the factors contributing to the flexibility of four-way junctions. The stacked-X conformation maximizes base-stacking but causes unfavorable repulsion between phosphate groups, whereas the 135 degrees -rotated "crossed" conformation minimizes electrostatic clashes at the expense of reduced base-stacking. Despite the large rotation of the branches, both junction structures exhibit an antiparallel arrangement of the continuous strands and opposite polarity of the crossover strands.