Structural parameters affecting the kinetics of RNA hairpin formation

Utilizing RNA-seq Data to Infer Bacterial Transcription Termination Sites and Validate Predictions

The transcription termination process is an important part of the gene expression process in the cell. It has been studied extensively, but many aspects of the mechanism are not well understood. The widespread availability of experimental RNA-seq data from high-throughput experiments provides a unique opportunity to infer the end of the transcription units genome wide. This data is available for both Rho-dependent and Rho-independent termination pathways that drive transcription termination in bacteria. Our book chapter gives an overview of the current knowledge of Rho-independent transcription termination mechanisms and the prediction approaches currently deployed to infer the termination sites. Thereafter, we describe our method that uses cluster hairpins to detect Rho-independent transcription termination sites. These clusters are a group of hairpins that lies at <15 bp from each other and are together capable of enforcing the termination process. The idea of a group of hairpins being extensively used for transcription termination is new, and results show that at least 52% of the total cases are of this type, while in the remaining cases, a single strong hairpin is capable of driving transcription termination. The reads derived from the RNA-seq data for corresponding bacteria have been used to validate the predicted sites. The predictions that match these RNA-seq derived sites have higher confidence, and we find almost 98% of the predicted sites, including alternate termination sites, to match the RNA-seq data. We discuss the features of predicted hairpins in detail for a better understanding of the Rho-independent transcription termination mechanism in bacteria. We also explain how users can use the tools developed by us to do transcription terminator predictions and design their experiments through genome-level visualization of the transcription termination sites from the precomputed INTERPIN database.

INTERPIN: A repository for intrinsic transcription termination hairpins in bacteria

The large-scale detection of putative intrinsic transcription terminators is limited to only a few bacteria currently. We discovered a group of hairpins, called cluster hairpins, present within 15 nucleotides from each other. These are expected to work in tandem to cause intrinsic transcription termination (ITT), while the single hairpin can do the same alone. Therefore, exploring these ITT sites and the hairpins across bacterial genomes becomes highly desirable. INTERPIN is the largest archived collection of in silico inferred ITT hairpins in bacteria, covering 12745 bacterial genomes and encompassing ten bacterial phyla for ∼25 million hairpins. Users can obtain details on operons, individual cluster, and single ITT hairpins that were screened therein. Integrated Genome Viewer (IGV) software interactively visualizes hairpin secondary and tertiary structures in the genomic context. We also discuss statistics for the occurrence of cluster or single hairpins and other termination alternatives while showing the validation of predicted hairpins against in vivo detected hairpins. The database is freely available at http://pallab.cds.iisc.ac.in/INTERPIN/. INTERPIN (database and software) can make predictions for both AT and GC-rich genomes, which has not been achieved by any other program so far. It can also be used to improve genome annotation as well as to get predictions to improve the understanding of the ITT pathway by further analysis.

Clusters of hairpins induce intrinsic transcription termination in bacteria

Intrinsic transcription termination (ITT) sites are currently identified by locating single and double-adjacent RNA hairpins downstream of the stop codon. ITTs for a limited number of genes/operons in only a few bacterial genomes are currently known. This lack of coverage is a lacuna in the existing ITT inference methods. We have studied the inter-operon regions of 13 genomes covering all major phyla in bacteria, for which good quality public RNA-seq data exist. We identify ITT sites in 87% of cases by predicting hairpin(s) and validate against 81% of cases for which the RNA-seq derived sites could be calculated. We identify 72% of these sites correctly, with 98% of them located ≤ 80 bases downstream of the stop codon. The predicted hairpins form a cluster (when present < 15 bases) in two-thirds of the cases, the remaining being single hairpins. The largest number of clusters is formed by two hairpins, and the occurrence decreases exponentially with an increasing number of hairpins in the cluster. Our study reveals that hairpins form an effective ITT unit when they act in concert in a cluster. Their pervasiveness along with single hairpin terminators corroborates a wider utilization of ITT mechanisms for transcription control across bacteria.

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.

Structure-guided mutational analysis of gene regulation by the Bacillus subtilis pbuE adenine-responsive riboswitch in a cellular context

Riboswitches are a broadly distributed form of RNA-based gene regulation in Bacteria and, more rarely, Archaea and Eukarya. Most often found in the 5'-leader sequence of bacterial mRNAs, they are generally composed of two functional domains: a receptor (aptamer) domain that binds an effector molecule and a regulatory domain (or expression platform) that instructs the expression machinery. One of the most studied riboswitches is the Bacillus subtilis adenine-responsive pbuE riboswitch, which regulates gene expression at the transcriptional level, up-regulating expression in response to increased intracellular effector concentrations. In this work, we analyzed sequence and structural elements that contribute to efficient ligand-dependent regulatory activity in a co-transcriptional and cellular context. Unexpectedly, we found that the P1 helix, which acts as the antitermination element of the switch in this RNA, supported ligand-dependent activation of a reporter gene over a broad spectrum of lengths from 3 to 10 bp. This same trend was also observed using a minimal in vitro single-turnover transcription assay, revealing that this behavior is intrinsic to the RNA sequence. We also found that the sequences at the distal tip of the terminator not directly involved in alternative secondary structure formation are highly important for efficient regulation. These data strongly support a model in which the switch is highly localized to the P1 helix adjacent to the ligand-binding pocket that likely presents a local kinetic block to invasion of the aptamer by the terminator.

Tospovirus ambisense genomic RNA segments use almost complete repertoire of stable tetraloops in the intergenic region

The intergenic regions of the ambisense RNA segments of viruses from the Tospovirus genus form large extended RNA structures that regulate virus replication. Using comparative structure analysis, we show the presence of conserved alternative conformations at the apical parts of these structures. In one conformation, a branched Y-shape, the 5'-proximal hairpin arms are mostly capped by exceptionally stable tetraloop motifs. The tetraloop hairpins are folded in both virus and virus-complementary sense RNAs, and different tetraloops can functionally replace each other. Folding simulations show that the branched Y-shape structures can undergo a conformational transition to alternative extended rod-like conformations. Functional importance of both alternatives is supported by nucleotide covariations. The balanced equilibrium between alternative structures is evidenced by native gel electrophoresis of mutant RNA transcripts with shifted equilibria. The tetraloops play a role in the stability and dynamics of structures but may also be recognized by proteins involved in translation and/or replication.

Stem-loop structures can effectively substitute for an RNA pseudoknot in -1 ribosomal frameshifting

-1 Programmed ribosomal frameshifting (PRF) in synthesizing the gag-pro precursor polyprotein of Simian retrovirus type-1 (SRV-1) is stimulated by a classical H-type pseudoknot which forms an extended triple helix involving base-base and base-sugar interactions between loop and stem nucleotides. Recently, we showed that mutation of bases involved in triple helix formation affected frameshifting, again emphasizing the role of the triple helix in -1 PRF. Here, we investigated the efficiency of hairpins of similar base pair composition as the SRV-1 gag-pro pseudoknot. Although not capable of triple helix formation they proved worthy stimulators of frameshifting. Subsequent investigation of ∼30 different hairpin constructs revealed that next to thermodynamic stability, loop size and composition and stem irregularities can influence frameshifting. Interestingly, hairpins carrying the stable GAAA tetraloop were significantly less shifty than other hairpins, including those with a UUCG motif. The data are discussed in relation to natural shifty hairpins.

Predicting secondary structural folding kinetics for nucleic acids

We report a new computational approach to the prediction of RNA secondary structure folding kinetics. In this approach, each elementary kinetic step is represented as the transformation between two secondary structures that differ by a helix. Based on the free energy landscape analysis, we identify three types of dominant pathways and the rate constants for the kinetic steps: 1), formation; 2), disruption of a helix stem; and 3), helix formation with concomitant partial melting of a competing (incompatible) helix. The third pathway, termed the tunneling pathway, is the low-barrier dominant pathway for the conversion between two incompatible helices. Comparisons with experimental data indicate that this new method is quite reliable in predicting the kinetics for RNA secondary structural folding and structural rearrangements. The approach presented here may provide a robust first step for further systematic development of a predictive theory for the folding kinetics for large RNAs.

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.

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.

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

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

Influence of C-5 halogenation of uridines on hairpin versus duplex RNA folding

Halogenation of bases is a widespread method used for solving crystal structures of nucleic acids. However, this modification may have important consequences on RNA folding and thus on the success of crystallization. We have used a combination of UV thermal melting, steady-state fluorescence, X-ray crystallography, and gel electrophoresis techniques to study the influence of uridine halogenation (bromination or iodination) on the RNA folding. The HIV-1 Dimerization Initiation Site is an RNA hairpin that can adopt an alternative duplex conformation and was used as a model. We have shown that, unexpectedly, the RNA hairpin/duplex ratio is strongly dependent not only on the presence but also on the position of halogenation.

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.