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

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.

Ligand-dependent tRNA processing by a rationally designed RNase P riboswitch

We describe a synthetic riboswitch element that implements a regulatory principle which directly addresses an essential tRNA maturation step. Constructed using a rational in silico design approach, this riboswitch regulates RNase P-catalyzed tRNA 5′-processing by either sequestering or exposing the single-stranded 5′-leader region of the tRNA precursor in response to a ligand. A single base pair in the 5′-leader defines the regulatory potential of the riboswitch both in vitro and in vivo. Our data provide proof for prior postulates on the importance of the structure of the leader region for tRNA maturation. We demonstrate that computational predictions of ligand-dependent structural rearrangements can address individual maturation steps of stable non-coding RNAs, thus making them amenable as promising target for regulatory devices that can be used as functional building blocks in synthetic biology.

Regulation of the thiamine pyrophosphate (TPP)-sensing riboswitch in NMT1 mRNA from Neurospora crassa

The expression of Neurospora crassa NMT1 involved in thiamine pyrophosphate (TPP) metabolism is regulated at the level of mRNA splicing by a TPP-sensing riboswitch within the precursor NMT1 mRNA. Here, using the systematic helix-based computational method, we investigated the regulation of this riboswitch. We find that the function of the riboswitch does not depend on the transcription process. Whether TPP is present or not, the riboswitch predominately folds into the ON state, while the OFF state aptamer structure does not appear during transcription. Since the transition from the ON state to the aptamer structure is extremely slow, TPP may interact with the RNA before full formation of the aptamer structure, promoting the switch flipping. The potential to fully form helix P0 of the ON state is necessary to restore ligand-dependent gene control by the riboswitch.

Modulation and Visualization of EF-G Power Stroke During Ribosomal Translocation

During ribosome translocation, the elongation factor EF‐G undergoes large conformational change while maintaining its contact with the moving tRNA. We previously measured a power stroke accompanying EF‐G catalysis, which was consistent with structural studies. However, the role of power stroke in translocation fidelity remains unclear. Here, we report quantitative measurements of the power strokes of structurally modified EF‐Gs by using two different techniques and reveal the correlation between power stroke and translocation efficiency and fidelity. We discovered that the reduced power stroke only lowered the percentage of translocation but did not introduce translocation error. The established force ‐structure–function correlation for EF‐G indicates that power stroke drives ribosomal translocation, but the mRNA reading frame is probably maintained by ribosome itself. Furthermore, the microscope detection method reported here can be simply implemented for other biochemical applications.

Predicting Cotranscriptional Folding Kinetics For Riboswitch

On the basis of a helix-based transition rate model, we developed a new method for sampling cotranscriptional RNA conformational ensemble and the prediction of cotranscriptional folding kinetics. Applications to E. coli. SRP RNA and pbuE riboswitch indicate that the model may provide reliable predictions for the cotranscriptional folding pathways and population kinetics. For E. coli. SRP RNA, the predicted population kinetics and the folding pathway are consistent with the SHAPE profiles in the recent cotranscriptional SHAPE-seq experiments. For the pbuE riboswitch, the model predicts the transcriptional termination efficiency as a function of the force. The theoretical results show (a) a force-induced transition from the aptamer (antiterminator) to the terminator structure and (b) the different folding pathways for the riboswitch with and without the ligand (adenine). More specifically, without adenine, the aptamer structure emerges as a short-lived kinetic transient state instead of a thermodynamically stable intermediate state. Furthermore, from the predicted extension-time curves, the model identifies a series of conformational switches in the pulling process, where the predicted relative residence times for the different structures are in accordance with the experimental data. The model may provide a new tool for quantitative predictions of cotranscriptional folding kinetics, and results can offer useful insights into cotranscriptional folding-related RNA functions such as regulation of gene expression with riboswitches.

Insights into Stability and Folding of GNRA and UNCG Tetraloops Revealed by Microsecond Molecular Dynamics and Well-Tempered Metadynamics

RNA hairpins capped by 5'-GNRA-3' or 5'-UNCG-3' tetraloops (TLs) are prominent RNA structural motifs. Despite their small size, a wealth of experimental data, and recent progress in theoretical simulations of their structural dynamics and folding, our understanding of the folding and unfolding processes of these small RNA elements is still limited. Theoretical description of the folding and unfolding processes requires robust sampling, which can be achieved by either an exhaustive time scale in standard molecular dynamics simulations or sophisticated enhanced sampling methods, using temperature acceleration or biasing potentials. Here, we study structural dynamics of 5'-GNRA-3' and 5'-UNCG-3' TLs by 15-μs-long standard simulations and a series of well-tempered metadynamics, attempting to accelerate sampling by bias in a few chosen collective variables (CVs). Both methods provide useful insights. The unfolding and refolding mechanisms of the GNRA TL observed by well-tempered metadynamics agree with the (reverse) folding mechanism suggested by recent replica exchange molecular dynamics simulations. The orientation of the glycosidic bond of the GL4 nucleobase is critical for the UUCG TL folding pathway, and our data strongly support the hypothesis that GL4-anti forms a kinetic trap along the folding pathway. Along with giving useful insight, our study also demonstrates that using only a few CVs apparently does not capture the full folding landscape of the RNA TLs. Despite using several sophisticated selections of the CVs, formation of the loop appears to remain a hidden variable, preventing a full convergence of the metadynamics. Finally, our data suggest that the unfolded state might be overstabilized by the force fields used.

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.

Design criteria for synthetic riboswitches acting on transcription

Riboswitches are RNA-based regulators of gene expression composed of a ligand-sensing aptamer domain followed by an overlapping expression platform. The regulation occurs at either the level of transcription (by formation of terminator or antiterminator structures) or translation (by presentation or sequestering of the ribosomal binding site). Due to a modular composition, these elements can be manipulated by combining different aptamers and expression platforms and therefore represent useful tools to regulate gene expression in synthetic biology. Using computationally designed theophylline-dependent riboswitches we show that 2 parameters, terminator hairpin stability and folding traps, have a major impact on the functionality of the designed constructs. These have to be considered very carefully during design phase. Furthermore, a combination of several copies of individual riboswitches leads to a much improved activation ratio between induced and uninduced gene activity and to a linear dose-dependent increase in reporter gene expression. Such serial arrangements of synthetic riboswitches closely resemble their natural counterparts and may form the basis for simple quantitative read out systems for the detection of specific target molecules in the cell.

Effect of loop composition on the stability and folding kinetics of RNA hairpins with large loops

RNA hairpins are ubiquitous structural elements in biological RNAs, where they have the potential to regulate RNA folding and interactions with other molecules. There are established methods for predicting the thermodynamic stability of an RNA hairpin, but there are still relatively few detailed examinations of the kinetics of folding. Nonetheless, several recent studies indicate that hairpin folding does not proceed via a simple two-state model. Here, we monitor fluorescence from hairpins constructed as molecular beacons in ensemble, fluorescence correlation spectroscopy, and stopped-flow experiments to describe the folding of RNA hairpins with long (15 nucleotide) loops. Our results show that folding of these hairpins occurs through more than two states and that the mechanism of folding includes a fast intermediate phase observed on the tens of microseconds time scale and a slow phase, attributed to formation of the native folded hairpin loop and stem, observed on the milliseconds time scale. The composition of the RNA loop determines the time scale of intermediate and native folded states. Hairpins with a polyuracil loop sequence exhibit slower relaxation of the intermediate state and faster relaxation of the native folded state when compared to that of hairpins with cytosine or adenine in the loop. We hypothesize this composition dependence could be attributed to nucleobase stacking in cytosine and adenine containing regions of the loop, which would be absent in hairpins containing polyuracil loops. Such base stacking could destabilize the intermediate folds, thereby speeding the relaxation of the intermediate relative to similar sized hairpins with no base stacking in the loop. Likewise, the lower intermediate stability could prolong the relaxation of the native folded state.

RNA folding: structure prediction, folding kinetics and ion electrostatics

Beyond the "traditional" functions such as gene storage, transport and protein synthesis, recent discoveries reveal that RNAs have important "new" biological functions including the RNA silence and gene regulation of riboswitch. Such functions of noncoding RNAs are strongly coupled to the RNA structures and proper structure change, which naturally leads to the RNA folding problem including structure prediction and folding kinetics. Due to the polyanionic nature of RNAs, RNA folding structure, stability and kinetics are strongly coupled to the ion condition of solution. The main focus of this chapter is to review the recent progress in the three major aspects in RNA folding problem: structure prediction, folding kinetics and ion electrostatics. This chapter will introduce both the recent experimental and theoretical progress, while emphasize the theoretical modelling on the three aspects in RNA folding.

Counterion and polythymidine loop-length-dependent folding and thermodynamic stability of DNA hairpins reveal the unusual counterion-dependent stability of tetraloop hairpins

Stem-loop DNA hairpins containing a 5-base-pair (bp) stem and single-stranded polythymidine loop were investigated using thermodynamic melting analysis and stopped-flow kinetics. These studies revealed the thermodynamic stability and folding kinetics as a function of loop length and counterion concentration. Our results show the unusually high thermodynamic stability for tetraloop or 4 poly(dT) loop hairpin as compared with longer loop length hairpins. Furthermore, this exceptional stability is highly counterion-dependent. For example, in the higher counterion concentration regime of 50 mM NaCl and above, the tetraloop hairpin displays enhanced stability as compared with longer loop length hairpins. However, at lower counterion concentration of 25 mM NaCl and below, the thermal stability of tetraloop hairpin is consistent with the longer loop hairpins. The enhanced stability of tetraloop hairpins at higher counterion concentration can be explained on the basis of the combined entropic effect of loop closure as well as base stacking in the loop regions. The stability of longer loop length hairpins at all counterion concentrations as well as tetraloop hairpin at lower counterion concentration can be explained on the basis of entropic effect of loop closure alone. The thermodynamic parameters at lower and higher counterion concentrations were determined to quantify the enhanced stability of base-stacking effects occurring at higher counterion concentrations. For example, for 100 mM NaCl, excess Gibbs energy and enthalpy due to base stacking within the tetraloops were measured to be -1.2 ± 0.14 and -3.28 ± 0.32 kcal/mol, respectively, whereas, no excess of Gibbs energy and enthalpy was observed for 0, 5, 10, and 25 mM NaCl. These findings suggest significant base-stacking interactions occurring in the loop region of the tetraloop hairpins at higher counterion concentration and less significant base-stacking interactions in the lower counterion concentration regime. We suggest that at higher counterion concentrations, hydrophobic collapse of the nucleotides in the loop may be enhanced due to the increased polarity of the solvent, thereby enhancing base-stacking interactions that contribute to unusually high stability.

Investigating the thermodynamics of UNCG tetraloops using infrared spectroscopy

Using infrared (IR) absorption spectroscopy, we have explored the folding thermodynamics of the UNCG class of RNA hairpin tetraloops (N = U, A, C, or G). Without the need to introduce non-native probes, IR spectroscopy makes it possible to distinguish specific structural elements such as base pairing versus base stacking or loop versus stem motions. Our results show that different structural components exhibit different thermodynamics. Specifically, we have found that tetraloop melting proceeds in a thermally sequential fashion, where base pairing in the stem is disrupted before (i.e., at a lower temperature) base stacking along the entire chain. In addition, for N = A, our data argue that the structure immediately surrounding adenine is particularly stable and melts at a higher temperature than either base-pairing or base-stacking interactions. Taken together, these results suggest that hairpin loop formation is not a simple two-state process, even in the equilibrium limit.

Panorama of DNA hairpin folding observed via diffusion-decelerated fluorescence correlation spectroscopy

A wide kinetic range from ∼0.1 μs to 1 s is offered by diffusion-decelerated FCS, allowing simultaneously monitoring multi-kinetic components.

Importance of diffuse metal ion binding to RNA

RNAs are highly charged polyanionic molecules. RNA structure and function are strongly correlated with the ionic condition of the solution. The primary focus of this article is on the role of diffusive ions in RNA folding. Due to the long-range nature of electrostatic interactions, the diffuse ions can contribute significantly to RNA structural stability and folding kinetics. We present an overview of the experimental findings as well as the theoretical developments on the diffuse ion effects in RNA folding. This review places heavy emphasis on the effect of magnesium ions. Magnesium ions play a highly efficient role in stabilizing RNA tertiary structures and promoting tertiary structural folding. The highly efficient role goes beyond the mean-field effect such as the ionic strength. In addition to the effects of specific ion binding and ion dehydration, ion-ion correlation for the diffuse ions can contribute to the efficient role of the multivalent ions such as the magnesium ions in RNA folding.

Perturbation-based Markovian transmission model for probing allosteric dynamics of large macromolecular assembling: a study of GroEL-GroES

Large macromolecular assemblies are often important for biological processes in cells. Allosteric communications between different parts of these molecular machines play critical roles in cellular signaling. Although studies of the topology and fluctuation dynamics of coarse-grained residue networks can yield important insights, they do not provide characterization of the time-dependent dynamic behavior of these macromolecular assemblies. Here we develop a novel approach called Perturbation-based Markovian Transmission (PMT) model to study globally the dynamic responses of the macromolecular assemblies. By monitoring simultaneous responses of all residues (>8,000) across many (>6) decades of time spanning from the initial perturbation until reaching equilibrium using a Krylov subspace projection method, we show that this approach can yield rich information. With criteria based on quantitative measurements of relaxation half-time, flow amplitude change, and oscillation dynamics, this approach can identify pivot residues that are important for macromolecular movement, messenger residues that are key to signal mediating, and anchor residues important for binding interactions. Based on a detailed analysis of the GroEL-GroES chaperone system, we found that our predictions have an accuracy of 71-84% judged by independent experimental studies reported in the literature. This approach is general and can be applied to other large macromolecular machineries such as the virus capsid and ribosomal complex.

Single-molecule fluorescence resonance energy transfer assays reveal heterogeneous folding ensembles in a simple RNA stem-loop

We have examined the folding ensembles present in solution for a series of RNA oligonucleotides that encompass the replicase translational operator stem-loop of the RNA bacteriophage MS2. Single-molecule (SM) fluorescence assays suggest that these RNAs exist in solution as ensembles of differentially base-paired/base-stacked states at equilibrium. There are two distinct ensembles for the wild-type sequence, implying the existence of a significant free energy barrier between "folded" and "unfolded" ensembles. Experiments with sequence variants are consistent with an unfolding mechanism in which interruptions to base-paired duplexes, in this example by the single-stranded loop and a single-base bulge in the base-paired stem, as well as the free ends, act as nucleation points for unfolding. The switch between folded and unfolded ensembles is consistent with a transition that occurs when all base-pairing and/or base-stacking interactions that would orientate the legs of the RNA stem are broken. Strikingly, a U-to-C replacement of a residue in the loop, which creates a high-affinity form of the operator for coat protein binding, results in dramatically different (un)folding behaviour, revealing distinct subpopulations that are either stabilised or destabilised with respect to the wild-type sequence. This result suggests additional reasons for selection against the C-variant stem-loop in vivo and provides an explanation for the increased affinity.

RNA folding: conformational statistics, folding kinetics, and ion electrostatics

RNA folding is a remarkably complex problem that involves ion-mediated electrostatic interaction, conformational entropy, base pairing and stacking, and noncanonical interactions. During the past decade, results from a variety of experimental and theoretical studies pointed to (a) the potential ion correlation effect in Mg2+-RNA interactions, (b) the rugged energy landscapes and multistate RNA folding kinetics even for small RNA systems such as hairpins and pseudoknots, (c) the intraloop interactions and sequence-dependent loop free energy, and (d) the strong nonadditivity of chain entropy in RNA pseudoknot and other tertiary folds. Several related issues, which have not been thoroughly resolved, require combined approaches with thermodynamic and kinetic experiments, statistical mechanical modeling, and all-atom computer simulations.

Perturbation-based Markovian Transmission Model for macromolecular machinery in cell

The study of the dynamics of a complex system is an important problem that includes large macromolecular complexes, molecular interaction networks, and cell functional modules. Large macromolecular complexes in cellular machinery can be modeled as a connected network, as in the elastic or Gaussian network models as demonstrated by Bahar and colleagues. Here we propose the Perturbation-based Markovian Transmission Model for studying the dynamics of signal transmission in macromolecular machinery. The initial perturbation is transmitted by a Markovian processes, and the dynamics of the probability flow is analytically solved using the master equation. Due to the large size of macromolecular complexes, it is very difficult to obtain analytical time-dependent Markovian dynamics of all atoms from the first perturbation until stationary state. To overcome it, we decrease the level of complexity of the transition matrix using a Krylov subspace method. This method is equivalent to integrating all eigen modes, and we show it can provide a globally accurate solution to the dynamics problem of signal transmission for very large macromolecular complexes with reasonable computational time. We give results of the dynamics of the GroEL-GroES chaperone system by applying uniform perturbation to all residues. We are able to identify experimentally found important residues and provide a set of predicted pivot, messenger, and effector residues, each with distinct dynamic behavior. Further results of selective perturbation on the surface of ATP binding pocket identifies the path of maximal probability flow of signal.

Salt dependence of nucleic acid hairpin stability

Single-stranded junctions/loops are frequently occurring structural motifs in nucleic acid structures. Due to the polyanionic nature of the nucleic acid backbone, metal ions play a crucial role in the loop stability. Here we use the tightly bound ion theory, which can account for the possible ion correlation and ensemble (fluctuation) effects, to predict the ion-dependence of loop and stem-loop (hairpin) free energies. The predicted loop free energy is a function of the loop length, the loop end-to-end distance, and the ion (Na(+) and Mg(2+) in this study) concentrations. Based on the statistical mechanical calculations, we derive a set of empirical formulas for the loop thermodynamic parameters as functions of Na(+) and Mg(2+) concentrations. For three specific types of loops, namely, hairpin, bulge, and internal loops, the predicted free energies agree with the experimental data. Further applications of these empirical formulas to RNA and DNA hairpin stability lead to good agreements with the available experimental data. Our results indicate that the ion-dependent loop stability makes significant contribution to the overall ion-dependence of the hairpin stability.

Review fluorescence correlation spectroscopy for probing the kinetics and mechanisms of DNA hairpin formation

This article reviews the application of fluorescence correlation spectroscopy (FCS) and related techniques to the study of nucleic acid hairpin conformational fluctuations in free aqueous solutions. Complimentary results obtained using laser-induced temperature jump spectroscopy, single-molecule fluorescence spectroscopy, optical trapping, and biophysical theory are also discussed. The studies cited reveal that DNA and RNA hairpin folding occurs by way of a complicated reaction mechanism involving long- and short-lived reaction intermediates. Reactions occurring on the subnanoseconds to seconds time scale have been observed, pointing out the need for experimental techniques capable of probing a broad range of reaction times in the study of such complex, multistate reactions.

Loop dependence of the stability and dynamics of nucleic acid hairpins

Hairpin loops are critical to the formation of nucleic acid secondary structure, and to their function. Previous studies revealed a steep dependence of single-stranded DNA (ssDNA) hairpin stability with length of the loop (L) as approximately L(8.5 +/- 0.5), in 100 mM NaCl, which was attributed to intraloop stacking interactions. In this article, the loop-size dependence of RNA hairpin stabilities and their folding/unfolding kinetics were monitored with laser temperature-jump spectroscopy. Our results suggest that similar mechanisms stabilize small ssDNA and RNA loops, and show that salt contributes significantly to the dependence of hairpin stability on loop size. In 2.5 mM MgCl2, the stabilities of both ssDNA and RNA hairpins scale as approximately L(4 +/- 0.5), indicating that the intraloop interactions are weaker in the presence of Mg2+. Interestingly, the folding times for ssDNA hairpins (in 100 mM NaCl) and RNA hairpins (in 2.5 mM MgCl2) are similar despite differences in the salt conditions and the stem sequence, and increase similarly with loop size, approximately L(2.2 +/- 0.5) and approximately L(2.6 +/- 0.5), respectively. These results suggest that hairpins with small loops may be specifically stabilized by interactions of the Na+ ions with the loops. The results also reinforce the idea that folding times are dominated by an entropic search for the correct nucleating conformation.

Elucidation of the RNA target of linezolid by using a linezolid-neomycin B heteroconjugate and genomic SELEX

A covalently modified heteroconjugate between linezolid and neomycin B leads to an enhanced and more specific binding affinity to hairpin RNA targets in comparison to neomycin B itself. This heteroconjugate was used as a lure to select linezolid-specific hairpin RNA from an Escherichia coli genome RNA. The selected RNA obtained after eight cycles not only has typical stem-loop structures but also includes known sequences of the linezolid binding site. The results of RNA footprinting show that the binding site of the heteroconjugate encompasses both stem and loop regions, suggesting that the possible binding site for linezolid is in the terminal loop. In addition, findings from application of a surface plasmon resonance assay clearly demonstrate that linezolid binds to selected hairpin RNA in a highly specific manner with a low millimolar affinity. The results suggest that heteroconjugates might represent a generally useful approach in studies aimed at uncovering loop-specific RNA binding ligands that would be otherwise difficult to identify owing to their weak affinities.

RNA helix stability in mixed Na+/Mg2+ solution

A recently developed tightly bound ion model can account for the correlation and fluctuation (i.e., different binding modes) of bound ions. However, the model cannot treat mixed ion solutions, which are physiologically relevant and biologically significant, and the model was based on B-DNA helices and thus cannot directly treat RNA helices. In the present study, we investigate the effects of ion correlation and fluctuation on the thermodynamic stability of finite length RNA helices immersed in a mixed solution of monovalent and divalent ions. Experimental comparisons demonstrate that the model gives improved predictions over the Poisson-Boltzmann theory, which has been found to underestimate the roles of multivalent ions such as Mg2+ in stabilizing DNA and RNA helices. The tightly bound ion model makes quantitative predictions on how the Na+-Mg2+ competition determines helix stability and its helix length-dependence. In addition, the model gives empirical formulas for the thermodynamic parameters as functions of Na+/Mg2+ concentrations and helix length. Such formulas can be quite useful for practical applications.