Self-induced structural switches in RNA

The Bacillaceae-1 RNA motif comprises two distinct classes

Non-coding RNAs are key regulatory players in bacteria. Many computationally predicted non-coding RNAs, however, lack functional associations. An example is the Bacillaceae-1 RNA motif, whose Rfam model consists of two hairpin loops. We find the motif conserved in nine of 13 non-pathogenic strains of the genus Bacillus but only in one pathogenic strain. To elucidate functional characteristics, we studied 118 hits of the Rfam model in 11 Bacillus spp. and found two distinct classes based on the ensemble diversity of their RNA secondary structure and the genomic context concerning the ribosomal RNA (rRNA) cluster. Forty hits are associated with the rRNA cluster, of which all 19 hits upstream flanking of 16S rRNA have a reverse complementary structure of low structural diversity. Fifty-two hits have large ensemble diversity, of which 38 are located between two coding genes. For eight hits in Bacillus subtilis, we investigated public expression data under various conditions and observed either the forward or the reverse complementary motif expressed. Five hits are associated with the rRNA cluster. Four of them are located upstream of the 16S rRNA and are not transcriptionally active, but instead, their reverse complements with low structural diversity are expressed together with the rRNA cluster. The three other hits are located between two coding genes in non-conserved genomic loci. Two of them are independently expressed from their surrounding genes and are structurally diverse. In summary, we found that Bacillaceae-1 RNA motifs upstream flanking of ribosomal RNA clusters tend to have one stable structure with the reverse complementary motif expressed in B. subtilis. In contrast, a subgroup of intergenic motifs has the thermodynamic potential for structural switches.

The 23S Ribosomal RNA From Pyrococcus furiosus Is Circularly Permuted

Synthesis and assembly of ribosomal components are fundamental cellular processes and generally well-conserved within the main groups of organisms. Yet, provocative variations to the general schemes exist. We have discovered an unusual processing pathway of pre-rRNA in extreme thermophilic archaea exemplified by Pyrococcus furiosus. The large subunit (LSU) rRNA is produced as a circularly permuted form through circularization followed by excision of Helix 98. As a consequence, the terminal domain VII that comprise the binding site for the signal recognition particle is appended to the 5´ end of the LSU rRNA that instead terminates in Domain VI carrying the Sarcin-Ricin Loop, the primary interaction site with the translational GTPases. To our knowledge, this is the first example of a true post-transcriptional circular permutation of a main functional molecule and the first example of rRNA fragmentation in archaea.

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.

Characterization of G-Quadruplex/Hairpin Transitions of RNAs by 19 F NMR Spectroscopy

2'-O-[(4-Trifluoromethyl-triazol-1-yl)methyl] reporter groups have been incorporated into guanosine-rich RNA models (including a known bistable Qd/Hp RNA and two G-rich regions of mRNA of human prion protein, PrP) and applied for the ¹⁹ F NMR spectroscopic characterization of plausible G-quadruplex/hairpin (Qd/Hp) transitions in these RNA structures. For the synthesis of the CF₃ -labeled RNAs, phosphoramidite building blocks of 2'-O-[(4-CF₃ -triazol-1-yl)methyl] nucleosides (cytidine, adenosine, and guanosine) were prepared and used as an integral part of the standard solid-phase RNA synthesis. The obtained ¹⁹ F NMR spectra supported the usual characterization data (obtained by UV- and CD-melting profiles and by ¹ H NMR spectra of the imino regions) and additionally gave more detailed information on the Qd/Hp transitions. The molar fractions of the secondary structural species (Qd, Hp) upon thermal denaturation and under varying ionic conditions could be determined from the intensities and shifts of the ¹⁹ F NMR signals. For a well-behaved Qd/Hp transition, thermodynamic parameters could be extracted.

Evolution of helix formation in the ribosomal Internal Transcribed Spacer 2 (ITS2) and its significance for RNA secondary structures

Helices are the most common elements of RNA secondary structure. Despite intensive investigations of various types of RNAs, the evolutionary history of the formation of new helices (novel helical structures) remains largely elusive. Here, by studying the nuclear ribosomal Internal Transcribed Spacer 2 (ITS2), a fast-evolving part of the eukaryotic nuclear ribosomal operon, we identify two possible types of helix formation: one type is "dichotomous helix formation"--transition from one large helix to two smaller helices by invagination of the apical part of a helix, which significantly changes the shape of the original secondary structure but does not increase its complexity (i.e., the total length of the RNA). An alternative type is "lateral helix formation"--origin of an extra helical region by the extension of a bulge loop or a spacer in a multi-helix loop of the original helix, which does not disrupt the pre-existing structure but increases RNA size. Moreover, we present examples from the RNA sequence literature indicating that both types of helix formation may have implications for RNA evolution beyond ITS2.

An RNA hairpin to G-quadruplex conformational transition

RNA molecules can fold into noncanonical structures such as the four-stranded structures known as G-quadruplexes. G-quadruplexes in the transcriptome have recently emerged as relevant regulatory elements of gene expression. Conformational transitions in RNA molecules offer an important way to regulate their biological functions. Here we report on the competition between a canonical hairpin structure and a G-quadruplex structure within an RNA molecule. We show that the conformational preference strongly depends on the relative amounts of mono- and divalent metal ions present in solution. In our system, the G-quadruplex, whose formation is not predicted by available predictive RNA folding programs, is the major conformer at physiologically relevant K(+) and Mg(2+) concentrations. Furthermore, we show that a synthetic small molecule can displace the structural dynamic equilibrium in favor of the hairpin conformer. This work highlights a new and important level of complexity in RNA folding that could be relevant to the biological functions and targeting of RNAs comprising G-quadruplex motifs.

Computational analysis of noncoding RNAs

Noncoding RNAs have emerged as important key players in the cell. Understanding their surprisingly diverse range of functions is challenging for experimental and computational biology. Here, we review computational methods to analyze noncoding RNAs. The topics covered include basic and advanced techniques to predict RNA structures, annotation of noncoding RNAs in genomic data, mining RNA-seq data for novel transcripts and prediction of transcript structures, computational aspects of microRNAs, and database resources.

Kinetic mechanism of conformational switch between bistable RNA hairpins

Transitions between the different conformational states play a critical role in many RNA catalytic and regulatory functions. In this study, we use the Kinetic Monte Carlo method to investigate the kinetic mechanism for the conformational switches between bistable RNA hairpins. We find three types of conformational switch pathways for RNA hairpins: refolding after complete unfolding, folding through basepair-exchange pathways and through pseudoknot-assisted pathways, respectively. The result of the competition between the three types of pathways depends mainly on the location of the rate-limiting base stacks (such as the GC base stacks) in the structures. Depending on the structural relationships between the two bistable hairpins, the conformational switch can follow single or multiple dominant pathways. The predicted folding pathways are supported by the activation energy results derived from the Arrhenius plot as well as the NMR spectroscopy data.

The importance of peripheral sequences in determining the metal selectivity of an in vitro-selected Co(2+) -dependent DNAzyme

DNAzymes are catalytically active DNA molecules that use metal cofactors for their enzymatic functions. While a growing number of DNAzymes with diverse functions and metal selectivities have been reported, the relationships between metal ion selectivity, conserved sequences and structures responsible for selectivity remain to be elucidated. To address this issue, we report biochemical assays of a family of previously reported in vitro selected DNAzymes. This family includes the clone 11 DNAzyme, which was isolated by positive and negative selection, and the clone 18 DNAzyme, which was isolated by positive selection alone. The clone 11 DNAzyme has a higher selectivity for Co(2+) over Pb(2+) compared with clone 18. The reasons for this difference are explored here through phylogenetic comparison, mutational analysis and stepwise truncation. A novel DNAzyme truncation method incorporated a nick in the middle of the DNAzyme to allow for truncation close to the nicked site while preserving peripheral sequences at both ends of the DNAzyme. The results demonstrate that peripheral sequences within the substrate binding arms, most notably the stem loop, loop II, are sufficient to restore its selectivity for Co(2+) over Pb(2+) to levels observed in clone 11. A comparison of these sequences' secondary structures and Co(2+) selectivities suggested that metastable structures affect metal ion selectivity. The Co(2+) selectivity of the clone 11 DNAzyme showed that the metal ion binding and selectivities of small, in vitro selected DNAzymes may be more complex than previously appreciated, and that clone 11 may be more similar to larger ribozymes than to other small DNAzymes in its structural complexity and behavior. These factors should be taken into account when metal-ion selectivity is required in rationally designed DNAzymes and DNAzyme-based biosensors.

Structure and stability of RNA/RNA kissing complex: with application to HIV dimerization initiation signal

We develop a statistical mechanical model to predict the structure and folding stability of the RNA/RNA kissing-loop complex. One of the key ingredients of the theory is the conformational entropy for the RNA/RNA kissing complex. We employ the recently developed virtual bond-based RNA folding model (Vfold model) to evaluate the entropy parameters for the different types of kissing loops. A benchmark test against experiments suggests that the entropy calculation is reliable. As an application of the model, we apply the model to investigate the structure and folding thermodynamics for the kissing complex of the HIV-1 dimerization initiation signal. With the physics-based energetic parameters, we compute the free energy landscape for the HIV-1 dimer. From the energy landscape, we identify two minimal free energy structures, which correspond to the kissing-loop dimer and the extended-duplex dimer, respectively. The results support the two-step dimerization process for the HIV-1 replication cycle. Furthermore, based on the Vfold model and energy minimization, the theory can predict the native structure as well as the local minima in the free energy landscape. The root-mean-square deviations (RMSDs) for the predicted kissing-loop dimer and extended-duplex dimer are ~3.0 Å. The method developed here provides a new method to study the RNA/RNA kissing complex.

Evolution of a helper virus-derived, ribosome binding translational enhancer in an untranslated satellite RNA of Turnip crinkle virus

SatC is a noncoding subviral RNA associated with Turnip crinkle virus (TCV). A 100-nt stretch in the 3' UTR of TCV contains three hairpins and two pseudoknots that fold into a tRNA-shaped structure (TSS) that binds 80S ribosomes. The 3' half of satC is derived from TCV and contains 6-nt differences in the TSS-analogous region. SatC binds poorly to 80S ribosomes, and molecular modeling that predicted the 3D structure of the TSS did not predict a similar structure for satC. When the satC TSS region was step-wise converted to the original TCV TSS bases, ribosome binding increased to TCV TSS levels without significantly affecting satC replication. However, mutant satC was less fit when accumulating in plants and gave rise to numerous second site changes that weakened one of two satC conformations. These results suggest that minor changes from the original TCV sequence in satC reflect requirements other than elimination of ribosome binding.

Thermodynamics of RNA melting, one base pair at a time

The melting of base pairs is a ubiquitous feature of RNA structural transitions, which are widely used to sense and respond to cellular stimuli. A recent study employing solution nuclear magnetic resonance (NMR) imino proton exchange spectroscopy provides a rare base-pair-specific view of duplex melting in the Salmonella FourU RNA thermosensor, which regulates gene expression in response to changes in temperature at the translational level by undergoing a melting transition. The authors observe "microscopic" enthalpy-entropy compensation--often seen "macroscopically" across a series of related molecular species--across base pairs within the same RNA. This yields variations in base-pair stabilities that are an order of magnitude smaller than corresponding variations in enthalpy and entropy. A surprising yet convincing link is established between the slopes of enthalpy-entropy correlations and RNA melting points determined by circular dichroism (CD), which argues that unfolding occurs when base-pair stabilities are equalized. A single AG-to-CG mutation, which enhances the macroscopic hairpin thermostability and folding cooperativity and renders the RNA thermometer inactive in vivo, spreads its effect microscopically throughout all base pairs in the RNA, including ones far removed from the site of mutation. The authors suggest that an extended network of hydration underlies this long-range communication. This study suggests that the deconstruction of macroscopic RNA unfolding in terms of microscopic unfolding events will require careful consideration of water interactions.

mRNA secondary structures fold sequentially but exchange rapidly in vivo

RNAs adopt defined structures to perform biological activities, and conformational transitions among alternative structures are critical to virtually all RNA-mediated processes ranging from metabolite-activation of bacterial riboswitches to pre-mRNA splicing and viral replication in eukaryotes. Mechanistic analysis of an RNA folding reaction in a biological context is challenging because many steps usually intervene between assembly of a functional RNA structure and execution of a biological function. We developed a system to probe mechanisms of secondary structure folding and exchange directly in vivo using self-cleavage to monitor competition between mutually exclusive structures that promote or inhibit ribozyme assembly. In previous work, upstream structures were more effective than downstream structures in blocking ribozyme assembly during transcription in vitro, consistent with a sequential folding mechanism. However, upstream and downstream structures blocked ribozyme assembly equally well in vivo, suggesting that intracellular folding outcomes reflect thermodynamic equilibration or that annealing of contiguous sequences is favored kinetically. We have extended these studies to learn when, if ever, thermodynamic stability becomes an impediment to rapid equilibration among alternative RNA structures in vivo. We find that a narrow thermodynamic threshold determines whether kinetics or thermodynamics govern RNA folding outcomes in vivo. mRNA secondary structures fold sequentially in vivo, but exchange between adjacent secondary structures is much faster in vivo than it is in vitro. Previous work showed that simple base-paired RNA helices dissociate at similar rates in vivo and in vitro so exchange between adjacent structures must occur through a different mechanism, one that likely involves facilitation of branch migration by proteins associated with nascent transcripts.

A conformational switch in the DiGIR1 ribozyme involved in release and folding of the downstream I-DirI mRNA

DiGIR1 is a group I-like cleavage ribozyme found as a structural domain within a nuclear twin-ribozyme group I intron. DiGIR1 catalyzes cleavage by branching at an Internal Processing Site (IPS) leading to formation of a lariat cap at the 5'-end of the 3'-cleavage product. The 3'-cleavage product is subsequently processed into an mRNA encoding a homing endonuclease. By analysis of combinations of 5'- and 3'-deletions, we identify a hairpin in the 5'-UTR of the mRNA (HEG P1) that is formed by conformational switching following cleavage. The formation of HEG P1 inhibits the reversal of the branching reaction, thus giving it directionality. Furthermore, the release of the mRNA is a consequence of branching rather than hydrolytic cleavage. A model is put forward that explains the release of the I-DirI mRNA with a lariat cap and a structured 5'-UTR as a direct consequence of the DiGIR1 branching reaction. The role of HEG P1 in GIR1 branching is reminiscent of that of hairpin P-1 in splicing of the Tetrahymena rRNA group I intron and illustrates a general principle in RNA-directed RNA processing.

RNA dynamics: it is about time

Many recently discovered RNA functions rely on highly complex multistep conformational transitions that occur in response to an array of cellular signals. These dynamics accompany and guide, for example, RNA cotranscriptional folding, ligand sensing and signaling, site-specific catalysis in ribozymes, and the hierarchically ordered assembly of ribonucleoproteins. RNA dynamics are encoded by both the inherent properties of RNA structure, spanning many motional modes with a large range of amplitudes and timescales, and external trigger factors, ranging from proteins, nucleic acids, metal ions, metabolites, and vitamins to temperature and even directional RNA biosynthesis itself. Here, we review recent advances in our understanding of RNA dynamics as highlighted by biophysical tools.

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.

Simulating RNA folding kinetics on approximated energy landscapes

We present a general computational approach to simulate RNA folding kinetics that can be used to extract population kinetics, folding rates and the formation of particular substructures that might be intermediates in the folding process. Simulating RNA folding kinetics can provide unique insight into RNA whose functions are dictated by folding kinetics and not always by nucleotide sequence or the structure of the lowest free-energy state. The method first builds an approximate map (or model) of the folding energy landscape from which the population kinetics are analyzed by solving the master equation on the map. We present results obtained using an analysis technique, map-based Monte Carlo simulation, which stochastically extracts folding pathways from the map. Our method compares favorably with other computational methods that begin with a comprehensive free-energy landscape, illustrating that the smaller, approximate map captures the major features of the complete energy landscape. As a result, our method scales to larger RNAs. For example, here we validate kinetics of RNA of more than 200 nucleotides. Our method accurately computes the kinetics-based functional rates of wild-type and mutant ColE1 RNAII and MS2 phage RNAs showing excellent agreement with experiment.

In silico design of small RNA switches

The discovery of natural RNA sensors that respond to a change in the environment by a conformational switch can be utilized for various biotechnological and nanobiotechnological advances. One class of RNA sensors is the riboswitch: an RNA genetic control element that is capable of sensing small molecules, responding to a deviation in ligand concentration with a structural change. Riboswitches are modularly built from smaller components. Computational methods can potentially be utilized in assembling these building block components and offering improvements in the biochemical design process. We describe a computational procedure to design RNA switches from building blocks with favorable properties. To achieve maximal throughput for genetic control purposes, future designer RNA switches can be assembled based on a computerized preprocessing buildup of the constituent domains, namely the aptamer and the expression platform in the case of a synthetic riboswitch. Conformational switching is enabled by the RNA versatility to possess two highly stable states that are energetically close to each other but topologically distinct, separated by an energy barrier between them. Initially, computer simulations can produce a list of short sequences that switch between two conformers when trigerred by point mutations or temperature. The short sequences should possess an additional desirable property; when these selected small RNA switch segments are attached to various aptamers, the ligand binding mechanism should replace the aforementioned event triggers, which will no longer be effective for crossing the energy barrier. In the assembled RNA sequence, energy minimization folding predictions should then show no difference between the folded structure of the entire sequence relative to the folded structure of each of its constituents. Moreover, energy minimization methods applied on the entire sequence could aid at this preprocessing stage by exhibiting high mutational robustness to capture the stability of the formed hairpin in the expression platform. The above computer-assisted assembly procedure together with application specific considerations may further be tailored for therapeutic gene regulation. Index Terms-Design of RNA switches, energy minimization methods, RNA folding predictions.

Shape similarity measures for the design of small RNA switches

Conformational switching in the secondary structure of RNAs has recently attracted considerable attention, fostered by the discovery of 'riboswitches' in living organisms. These are genetic control elements that were found in bacteria and offer a unique regulation mechanism based on switching between two highly stable states, separated by an energy barrier between them. In riboswitches, the energy barrier is crossed by direct metabolite binding, which facilitates regulation by allosteric means. However, other event triggers can cause switching to occur, such as single-point mutations and slight variations in temperature. Examples of switches with these event triggers have already been reported experimentally in the past. Here, the goal is to computationally design small RNA switches that rely on these triggers. Towards this end, our computer simulations utilize a variety of different similarity measures to assess the distances between an initial state and triggered states, based on the topology of the secondary structure itself. We describe these combined similarity measures that rely on both coarse-grained and fine-grained graph representations of the RNA secondary structure. As a result of our simulations, we provide some candidate sequences of approximately 30-50 nt, along with the exact triggers that drive the switching. The event triggers under consideration can be modelled by Zuker's mfold or the Vienna package. The proposed methodology that rely on shape measures can further be used to computationally generate more candidates by simulating various event triggers and calculating their effect on the shape.

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.

A pseudoknot in a preactive form of a viral RNA is part of a structural switch activating minus-strand synthesis

RNA can adopt different conformations in response to changes in the metabolic status of cells, which can regulate processes such as transcription, translation, and RNA cleavage. We previously proposed that an RNA conformational switch in an untranslated satellite RNA (satC) of Turnip crinkle virus (TCV) regulates initiation of minus-strand synthesis (G. Zhang, J. Zhang, A. T. George, T. Baumstark, and A. E. Simon, RNA 12:147-162, 2006). This model was based on the lack of phylogenetically inferred hairpins or a known pseudoknot in the "preactive" structure assumed by satC transcripts in vitro. We now provide evidence that a second pseudoknot (Psi(2)), whose disruption reduces satC accumulation in vivo and enhances transcription by the TCV RNA-dependent RNA polymerase in vitro, stabilizes the preactive satC structure. Alteration of either Psi(2) partner caused nearly identical structural changes, including single-stranded-specific cleavages in the pseudoknot sequences and strong cleavages in a distal element previously proposed to mediate the conformational switch. These results indicate that the preactive structure identified in vitro has biological relevance in vivo and support a requirement for this alternative structure and a conformational switch in high-level accumulation of satC in vivo.

The role of a metastable RNA secondary structure in hepatitis delta virus genotype III RNA editing

RNA editing plays a critical role in the life cycle of hepatitis delta virus (HDV). The host editing enzyme ADAR1 recognizes specific RNA secondary structure features around the amber/W site in the HDV antigenome and deaminates the amber/W adenosine. A previous report suggested that a branched secondary structure is necessary for editing in HDV genotype III. This branched structure, which is distinct from the characteristic unbranched rod structure required for HDV replication, was only partially characterized, and knowledge concerning its formation and stability was limited. Here, we examine the secondary structures, conformational dynamics, and amber/W site editing of HDV genotype III RNA using a miniaturized HDV genotype III RNA in vitro. Computational analysis of this RNA using the MPGAfold algorithm indicated that the RNA has a tendency to form both metastable and stable unbranched secondary structures. Moreover, native polyacrylamide gel electrophoresis demonstrated that this RNA forms both branched and unbranched rod structures when transcribed in vitro. As predicted, the branched structure is a metastable structure that converts readily to the unbranched rod structure. Only branched RNA was edited at the amber/W site by ADAR1 in vitro. The structural heterogeneity of HDV genotype III RNA is significant because not only are both conformations of the RNA functionally important for viral replication, but the ratio of the two forms could modulate editing by determining the amount of substrate RNA available for modification.

Structural parameters affecting the kinetics of RNA hairpin formation

There is little experimental knowledge on the sequence dependent rate of hairpin formation in RNA. We have therefore designed RNA sequences that can fold into either of two mutually exclusive hairpins and have determined the ratio of folding of the two conformations, using structure probing. This folding ratio reflects their respective folding rates. Changing one of the two loop sequences from a purine- to a pyrimidine-rich loop did increase its folding rate, which corresponds well with similar observations in DNA hairpins. However, neither changing one of the loops from a regular non-GNRA tetra-loop into a stable GNRA tetra-loop, nor increasing the loop size from 4 to 6 nt did affect the folding rate. The folding kinetics of these RNAs have also been simulated with the program ‘Kinfold’. These simulations were in agreement with the experimental results if the additional stabilization energies for stable tetra-loops were not taken into account. Despite the high stability of the stable tetra-loops, they apparently do not affect folding kinetics of these RNA hairpins. These results show that it is possible to experimentally determine relative folding rates of hairpins and to use these data to improve the computer-assisted simulation of the folding kinetics of stem–loop structures.

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.

A mutational analysis of the 8-17 deoxyribozyme core

The 8-17 deoxyribozyme is a small RNA-cleaving DNA enzyme of significant applicative interest. We measured the kinetics of over 60 variants of 8-17, mutated within the "core" region. The data were analyzed according to a conceptual framework in which deleterious substitutions can either decrease the stability of the reaction's transition state, or favor unreactive ground-state conformations. In agreement with earlier in vitro evolution studies, the most severe functional effects were observed upon mutating four conserved residues, whose role was further explored by replacing them with non-standard nucleotides. Removal or modification of individual functional groups on the A6 and G7 bases suggested that these residues are involved in a close-contact interaction and form a network of functionally important hydrogen bonds. Mutagenesis of residues C13 and G14 was less revealing, but argued strongly against a role of C13 as a general acid/base catalyst. The use of non-standard nucleotides also led to the identification of one deoxyribozyme variant that, under some ionic conditions, is substantially more active than the wild-type construct. Finally, the effects of mutations in the intramolecular "core stem" correlated only in part with changes in helical stability, suggesting that a stable stem is required but not sufficient for optimal activity.

Discontinuous subgenomic RNA synthesis in arteriviruses is guided by an RNA hairpin structure located in the genomic leader region

Nidoviruses produce an extensive 3'-coterminal nested set of subgenomic (sg) mRNAs, which are used to express structural proteins and sometimes accessory proteins. In arteriviruses and coronaviruses, these mRNAs contain a common 5' leader sequence, derived from the genomic 5' end. The joining of the leader sequence to different segments derived from the 3'-proximal part of the genome (mRNA bodies) presumably involves a unique mechanism of discontinuous minus-strand RNA synthesis in which base pairing between sense and antisense transcription-regulating sequences (TRSs) plays an essential role. The leader TRS is present in the loop of a hairpin structure that functions in sg mRNA synthesis. In this study, the minimal sequences in the 5'-proximal region of the Equine arteritis virus genome that are required for sg RNA synthesis were delimited through mutagenesis. A full-length cDNA clone was engineered in which this domain was duplicated, allowing us to make mutations and monitor their effects on sg RNA synthesis without seriously affecting genome replication and translation. The leader TRS present in the duplicated sequence was used and yielded novel sg mRNAs with significantly extended leaders. Our combined findings suggest that the leader TRS hairpin (LTH) and its immediate flanking sequences are essential for efficient sg RNA synthesis and form an independent functional entity that could be moved 300 nucleotides downstream of its original position in the genome. We hypothesize that a conformational switch in the LTH region regulates the role of the 5'-proximal region of the arterivirus genome in subgenomic RNA synthesis.