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

pH-dependent complex formation with TAR RNA and DNA: application towards logic gates

Nucleic acid-based logic gates have shown great potential in biotechnology, medicine as well as diagnostics. Herein, we have constructed pH-responsive logic devices by utilizing HIV-1 TAR hairpins in combination with a thiazole peptide that exhibits turn-on fluorescence upon interacting with TAR RNA or DNA. Based on this, INHIBIT-AND and YES-INHIBIT-AND logic gates were constructed in parallel. The pH alteration leads to conformational changes of the hairpin structure, enabling the construction of a multi-reset reusable logic system which could be developed for in vitro sensing of the HIV-1 viral RNA.

Kinetic Mechanism of RNA Helix-Terminal Basepairing-A Kinetic Minima Network Analysis

RNA functions are often kinetically controlled. The folding kinetics of RNAs involves global structural changes and local nucleotide movement, such as base flipping. The most elementary step in RNA folding is the closing and opening of a basepair. By integrating molecular dynamics simulation, master equation, and kinetic Monte Carlo simulation, we investigate the kinetics mechanism of RNA helix-terminal basepairing. The study reveals a six-state folding scheme with three dominant folding pathways of tens, hundreds, and thousands of nanoseconds of folding timescales, respectively. The overall kinetics is rate limited by the detrapping of a misfolded state with the overall folding time of 10^-5 s. Moreover, the analysis examines the different roles of the various driving forces, such as the basepairing and stacking interactions and the ion binding/dissociation effects on structural changes. The results may provide useful insights for developing a basepair opening/closing rate model and further kinetics models of large RNAs.

Understanding the kinetic mechanism of RNA single base pair formation

RNA functions are intrinsically tied to folding kinetics. The most elementary step in RNA folding is the closing and opening of a base pair. Understanding this elementary rate process is the basis for RNA folding kinetics studies. Previous studies mostly focused on the unfolding of base pairs. Here, based on a hybrid approach, we investigate the folding process at level of single base pairing/stacking. The study, which integrates molecular dynamics simulation, kinetic Monte Carlo simulation, and master equation methods, uncovers two alternative dominant pathways: Starting from the unfolded state, the nucleotide backbone first folds to the native conformation, followed by subsequent adjustment of the base conformation. During the base conformational rearrangement, the backbone either retains the native conformation or switches to nonnative conformations in order to lower the kinetic barrier for base rearrangement. The method enables quantification of kinetic partitioning among the different pathways. Moreover, the simulation reveals several intriguing ion binding/dissociation signatures for the conformational changes. Our approach may be useful for developing a base pair opening/closing rate model.

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.

TRANSAT-- method for detecting the conserved helices of functional RNA structures, including transient, pseudo-knotted and alternative structures

The prediction of functional RNA structures has attracted increased interest, as it allows us to study the potential functional roles of many genes. RNA structure prediction methods, however, assume that there is a unique functional RNA structure and also do not predict functional features required for in vivo folding. In order to understand how functional RNA structures form in vivo, we require sophisticated experiments or reliable prediction methods. So far, there exist only a few, experimentally validated transient RNA structures. On the computational side, there exist several computer programs which aim to predict the co-transcriptional folding pathway in vivo, but these make a range of simplifying assumptions and do not capture all features known to influence RNA folding in vivo. We want to investigate if evolutionarily related RNA genes fold in a similar way in vivo. To this end, we have developed a new computational method, Transat, which detects conserved helices of high statistical significance. We introduce the method, present a comprehensive performance evaluation and show that Transat is able to predict the structural features of known reference structures including pseudo-knotted ones as well as those of known alternative structural configurations. Transat can also identify unstructured sub-sequences bound by other molecules and provides evidence for new helices which may define folding pathways, supporting the notion that homologous RNA sequence not only assume a similar reference RNA structure, but also fold similarly. Finally, we show that the structural features predicted by Transat differ from those assuming thermodynamic equilibrium. Unlike the existing methods for predicting folding pathways, our method works in a comparative way. This has the disadvantage of not being able to predict features as function of time, but has the considerable advantage of highlighting conserved features and of not requiring a detailed knowledge of the cellular environment.

Many non-coding genes exert their function via an RNA structure which starts emerging while the RNA sequence is being transcribed from the genome. The resulting folding pathway is known to depend on a variety of features such as the transcription speed, the concentration of various ions and the binding of proteins and other molecules. Not all of these influences can be adequately captured by the existing computational methods which try to replicate what happens in vivo. So far, it has been challenging to experimentally investigate co-transcriptional folding pathways in vivo and only little data from in vitro experiments exists. In order to investigate if functionally similar RNA sequences from different organisms fold in a similar way, we have developed a new computational method, called Transat, which does not require the detailed computational modeling of the cellular environment. We show in a comprehensive analysis that our method is capable of detecting known structural features and provide evidence that structural features of the in vivo folding pathways have been conserved for several biologically interesting classes of RNA sequences.

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.

Structures, kinetics, thermodynamics, and biological functions of RNA hairpins

Most RNA comprises one strand and therefore can fold back on itself to form complex structures. At the heart of these structures is the hairpin, which is composed of a stem having Watson-Crick base pairing and a loop wherein the backbone changes directionality. First, we review the structure of hairpins including diversity in the stem, loop, and closing base pair. The function of RNA hairpins in biology is discussed next, including roles for isolated hairpins, as well as hairpins in the context of complex tertiary structures. We describe the kinetics and thermodynamics of hairpin folding including models for hairpin folding, folding transition states, and the cooperativity of folding. Lastly, we discuss some ways in which hairpins can influence the folding and function of tertiary structures, both directly and indirectly. RNA hairpins provide a simple means of controlling gene expression that can be understood in the language of physical chemistry.

Ligand-induced folding of the guanine-sensing riboswitch is controlled by a combined predetermined induced fit mechanism

All known guanine-sensing riboswitches regulate gene expression by specifically binding to guanine (G) or related analogs with high affinity to switch off transcription. The aptamers of this class of riboswitches are characterized by three helices (P1-P3), surrounding a central core of phylogenetically conserved nucleotides and a long-range loop-loop interaction. To gain more insight into the switching mechanism, we present here a comparison between the solution-state structures of the G-free and G-bound forms of the guanine aptamer from the xpt-pbuX operon of Bacillus subtilis, as derived from NMR chemical shifts and magnetic-field-induced residual dipolar couplings. The high-resolution NMR analysis shows the G-free aptamer is highly structured with parallel P2 and P3 helices and the long-range loop-loop interaction already present, implying that the structure is largely preformed to bind the ligand. Structural changes upon guanine binding are found to be localized to the central core. In the free state, the G-quadruple interaction and two base pairs of the P1 stem flanking the central core appear to be largely disordered. The ligand thus binds via a combined predetermined-induced fit mechanism, involving a previously unstructured five-residue loop of the J2-3 junction that folds over the ligand. These limited additional interactions within a preorganized setting possibly explain how the aptamer rapidly responds to ligand binding, which is necessary to switch the structural state of the expression platform within a narrow time frame before the RNA polymerase escapes the 5'-UTR.

Kinetics of RNA refolding in dynamic equilibrium by 1H-detected 15N exchange NMR spectroscopy

By implementing new NMR methods that were designed to map very slow exchange processes we have investigated and characterized the refolding kinetics of a thermodynamically stable 34mer RNA sequence in dynamic equilibrium. The RNA sequence was designed to undergo a topologically favored conformational exchange between different hairpin folds, serving as a model to estimate the minimal time required for more complex RNA folding processes. Chemically prepared RNA sequences with sequence-selective (15)N labels provided the required signal separation and allowed a straightforward signal assignment of the imino protons by HNN correlation experiments. The 2D version of the new (1)H-detected (15)N exchange spectroscopy (EXSY) pulse sequence provided cross-peaks for resonances belonging to different folds that interchange on the time scale of longitudinal relaxation of (15)N nuclei bound to imino protons. The 34mer RNA sequence exhibits two folds which exchange on the observable time scale (tau(obs) approximately T(1){(15)Nu} < 5 s) and a third fold which is static on this time scale. A 1D version of the (15)N exchange experiment allowed the measurement of the exchange rates between the two exchanging folds as a function of temperature and the determination of the corresponding activation energies E(a) and frequency factors A. We found that the refolding rates are strongly affected by an entropically favorable preorientation of the replacing strand. The activation energies are comparable to values obtained for the slow refolding of RNA sequences of similar thermodynamic stability but less favorable topology.

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.

On A Six-Dimensional Representation of RNA Secondary Structures

In this paper, we proposed a 6-D representation of RNA secondary structures. The use of the 6-D representation is illustrated by constructing structure invariants. Comparisons with the similarity/dissimilarity results based on 6-D representation for a set of RNA secondary structures, are considered to illustrate the use of our structure invariants based on the entries in derived sequence matrices restricted to a selected width of a band along the main diagonal.

Monomer and dimer of Chandipura virus unphosphorylated P-protein binds leader RNA differently: implications for viral RNA synthesis

Interaction of the leader RNA with the unphosphorylated P-protein has been proposed to play a key role in the transcription-replication transition of Chandipura virus, a model rhabdovirus. Electrophoretic mobility shift assay with the leader RNA and the unphosphorylated P-protein demonstrated existence of two distinct complexes in vitro. Measurements of stoichiometry indicate the protein monomer/RNA ratio to be 1:1 and 2:1 for faster and slower migrating bands, respectively. We have also observed a concentration-dependent oligomerization of the unphosphorylated P-protein, in sub-micromolar to low micromolar range. Sedimentation velocity, dynamic light scattering and large zone gel filtration experiments suggest a monomer-dimer-tetramer model of association. RNA binding experiments suggest that the two complexes assembled from one molecule of the leader RNA binding to either a protein monomer or a dimer. A truncated RNA consisting of a 3' region of the leader transcript exclusively formed the 1:1 complex, whereas a RNA consisting of only the 5' region forms the 2:1 complex exclusively. RNA binding experiments at different protein concentrations suggest that binding of the RNA comprising the 3' region weakens significantly at higher P(0) concentrations, whereas in contrast the binding of the RNA comprising the 5' region becomes modestly tighter. Implications of two different types of leader RNA-P-protein complexes in viral RNA synthesis are discussed.

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.