Exploring RNA structure by integrative molecular modelling

Problem of Domain/Building Block Preservation in the Evolution of Biological Macromolecules and Evolutionary Computation

Structurally and functionally isolated domains in biological macromolecular evolution, both natural and artificial, are largely similar to "schemata", building blocks (BBs), in evolutionary computation (EC). The problem of preserving in subsequent evolutionary searches the already found domains / BBs is well known and quite relevant in biology as well as in EC. Both biology and EC are seeing parallel and independent development of several approaches to identifying and preserving previously identified domains / BBs. First, we notice the similarity of DNA shuffling methods in synthetic biology and multi-parent recombination algorithms in EC. Furthermore, approaches to computer identification of domains in proteins that are being developed in biology can be aligned with BB identification methods in EC. Finally, approaches to chimeric protein libraries optimization in biology can be compared to evolutionary search methods based on probabilistic models in EC. We propose to validate the prospects of mutual exchange of ideas and transfer of algorithms and approaches between evolutionary systems biology and EC in these three principal directions. A crucial aim of this transfer is the design of new advanced experimental techniques capable of solving more complex problems of in vitro evolution.

Conformational adaptation of UNCG loops upon crowding

If the A-form helix is the major structural motif found in RNA, the loops that cap them constitute the second most important family of motifs. Among those, two are overrepresented, GNRA and UNCG tetraloops. Recent surveys of RNA structures deposited in the PDB show that GNRA and UNCG tetraloops can adopt tertiary folds that are very different from their canonical conformations, characterized by the presence of a U-turn of a Z-turn, respectively. Crystallographic data from both a lariat-capping (LC) ribozyme and a group II intron ribozyme reveal that a given UUCG tetraloop can adopt a distinct fold depending on its structural environment. Specifically, when the crystal packing applies relaxed constraints on the loop, the canonical Z-turn conformation is observed. In contrast, a highly packed environment induces "squashing" of the tetraloop by distorting its sugar-phosphate backbone in a specific way that expels the first and fourth nucleobases out of the loop, and falls in van der Waals distance of the last base pair of the helix, taking the place of the pair formed between the first and fourth residues in Z-turn loops. The biological relevance of our observations is supported by the presence of similarly deformed loops in the highly packed environment of the ribosome and in a complex between a dsRNA and a RNase III. The finding that Z-turn loops change conformation under higher molecular packing suggests that, in addition to their demonstrated role in stabilizing RNA folding, they may contribute to the three-dimensional structure of RNA by mediating tertiary interactions with distal residues.

A T7 RNA Polymerase Mutant Enhances the Yield of 5'-Thienoguanosine-Initiated RNAs

Spectroscopic methods, which are used to establish RNA structure-function relationships, require strategies for post-synthetic, site-specific incorporation of chemical probes into target RNAs. For RNAs larger than 50 nt, the enzymatic incorporation of a nucleoside or nucleotide monophosphate guanosine analogue (G analogue) at their 5'-end is routinely achieved by T7 RNA polymerase (T7RNAP)-mediated in vitro transcription (IVT) of the appropriate DNA template containing a GTP-initiating class III Φ6.5 promoter. However, when high G analogue:GTP ratios are used to bias G analogue incorporation at the 5'-end, RNA yield is compromised. Here, we show that the use of a T7RNAP P266L mutant in IVT with 10:1 thienoguanosine (^th G):GTP increased the percent incorporation and yield of 5'-^th G-initiated precursor tRNA for a net ≈threefold gain compared to IVT with wild-type T7RNAP. We also demonstrated that a one-pot multienzyme approach, consisting of transcription by T7RNAP P266L and post-transcriptional cleanup by polyphosphatase and an exonuclease, led to essentially near-homogeneous 5'-^th G-modified transcripts. This approach should be of broad utility in preparing 5'-modified RNAs.

A coarse-grained model for assisting the investigation of structure and dynamics of large nucleic acids by ion mobility spectrometry-mass spectrometry

Ion Mobility Spectrometry-Mass Spectrometry (IMS-MS) is a rapidly emerging tool for the investigation of nucleic acid structure and dynamics. IMS-MS determinations can provide valuable information regarding alternative topologies, folding intermediates, and conformational heterogeneities, which are not readily accessible to other analytical techniques. The leading strategies for data interpretation rely on computational and experimental approaches to correctly assign experimental observations to putative structures. A very effective strategy involves the application of molecular dynamics (MD) simulations to predict the structure of the analyte molecule, calculate its collision cross section (CCS), and then compare this computational value with the corresponding experimental data. While this approach works well for small nucleic acid species, analyzing larger nucleic acids of biological interest is hampered by the computational cost associated with capturing their extensive structure and dynamics in all-atom detail. In this report, we describe the implementation of a coarse graining (CG) approach to reduce the cost of the computational methods employed in the data interpretation workflow. Our framework employs a five-bead model to accurately represent each nucleotide in the nucleic acid structure. The beads are appropriately parameterized to enable the direct calculation of CCS values from CG models, thus affording the ability to pursue the analysis of larger, highly dynamic constructs. The validity of this approach was successfully confirmed by the excellent correlation between the CCS values obtained in parallel by all-atom and CG workflows.

Transfer RNA: From pioneering crystallographic studies to contemporary tRNA biology

Transfer RNAs (tRNAs) play a key role in protein synthesis as adaptor molecules between messenger RNA and protein sequences on the ribosome. Their discovery in the early sixties provoked a worldwide infatuation with the study of their architecture and their function in the decoding of genetic information. tRNAs are also emblematic molecules in crystallography: the determination of the first tRNA crystal structures represented a milestone in structural biology and tRNAs were for a long period the sole source of information on RNA folding, architecture, and post-transcriptional modifications. Crystallographic data on tRNAs in complex with aminoacyl-tRNA synthetases (aaRSs) also provided the first insight into protein:RNA interactions. Beyond the translation process and the history of structural investigations on tRNA, this review also illustrates the renewal of tRNA biology with the discovery of a growing number of tRNA partners in the cell, the involvement of tRNAs in a variety of regulatory and metabolic pathways, and emerging applications in biotechnology and synthetic biology.

Characterizing 3D RNA structure by single molecule FRET

The importance of elucidating the three dimensional structures of RNA molecules is becoming increasingly clear. However, traditional protein structural techniques such as NMR and X-ray crystallography have several important drawbacks when probing long RNA molecules. Single molecule Förster resonance energy transfer (smFRET) has emerged as a useful alternative as it allows native sequences to be probed in physiological conditions and allows multiple conformations to be probed simultaneously. This review serves to describe the method of generating a three dimensional RNA structure from smFRET data from the biochemical probing of the secondary structure to the computational refinement of the final model.

Speciation of a group I intron into a lariat capping ribozyme

The lariat-capping (LC) ribozyme is a natural ribozyme isolated from eukaryotic microorganisms. Despite apparent structural similarity to group I introns, the LC ribozyme catalyzes cleavage by a 2',5' branching reaction, leaving the 3' product with a 3-nt lariat cap that functionally substitutes for a conventional mRNA cap in the downstream pre-mRNA encoding a homing endonuclease. We describe the crystal structures of the precleavage and postcleavage LC ribozymes, which suggest that structural features inherited from group I ribozymes have undergone speciation due to profound changes in molecular selection pressure, ultimately giving rise to an original branching ribozyme family. The structures elucidate the role of key elements that regulate the activity of the LC ribozyme by conformational switching and suggest a mechanism by which the signal for branching is transmitted to the catalytic core. The structures also show how conserved interactions twist residues, forming the lariat to join chemical groups involved in branching.

cis-Acting 5' hammerhead ribozyme optimization for in vitro transcription of highly structured RNAs

RNA-mediated biological processes usually require precise definition of 5' and 3' ends. RNA ends obtained by in vitro transcription using T7 RNA polymerase are often heterogeneous in length and sequence. An efficient strategy to overcome these drawbacks consists of inserting an RNA with known boundaries in between two ribozymes, usually a 5' hammerhead and a 3' hepatitis delta virus ribozymes, that cleave off the desired RNA. In practice, folding of the three RNAs challenges each other, potentially preventing thorough processing. Folding and cleavage of the 5' hammerhead ribozyme relies on a sequence of nucleotides belonging to the central RNA making it more sensitive than the usual 3' hepatitis delta virus ribozyme. The intrinsic stability of the central RNA may thus prevent correct processing of the full transcript. Here, we present a method in which incorporation of a full-length hammerhead ribozyme with a specific tertiary interaction prevents alternative folding with the lariat capping GIR1 ribozyme and enables complete cleavage in the course of the transcription. This strategy may be transposable for in vitro transcription of any highly structured RNA.

Predicting helical topologies in RNA junctions as tree graphs

RNA molecules are important cellular components involved in many fundamental biological processes. Understanding the mechanisms behind their functions requires knowledge of their tertiary structures. Though computational RNA folding approaches exist, they often require manual manipulation and expert intuition; predicting global long-range tertiary contacts remains challenging. Here we develop a computational approach and associated program module (RNAJAG) to predict helical arrangements/topologies in RNA junctions. Our method has two components: junction topology prediction and graph modeling. First, junction topologies are determined by a data mining approach from a given secondary structure of the target RNAs; second, the predicted topology is used to construct a tree graph consistent with geometric preferences analyzed from solved RNAs. The predicted graphs, which model the helical arrangements of RNA junctions for a large set of 200 junctions using a cross validation procedure, yield fairly good representations compared to the helical configurations in native RNAs, and can be further used to develop all-atom models as we show for two examples. Because junctions are among the most complex structural elements in RNA, this work advances folding structure prediction methods of large RNAs. The RNAJAG module is available to academic users upon request.

The functional exchangeability of pk- and k-turns in RNA structure

Ribonuclease P RNA requires a sharply kinked RNA helix to make a loop-receptor interaction that creates the binding site for the substrate. In some forms of the ribozyme, this is accomplished by a k-turn, while others have a different element called the pk-turn. The structure of the pk-turn in RNase P of Thermotoga maritima is globally very similar to a k-turn, but lacks all the standard features of that structure, including long-range hydrogen bonds between the two helical arms. We show here that in an isolated RNA duplex, the pk-turn fails to adopt a tightly kinked structure, but rather is a flexible element. This suggests that the tertiary contacts of RNase P assist its folding into the required kinked structure. We find that we can replace the k-turn of the SAM-I riboswitch with the pk-turn, such that the resulting RNA retains its ability to bind SAM, although with lower affinity. We also find that we can replace the pk-turn of T. maritima RNase P with a standard k-turn (in either orientation) with retention of ribozyme activity. Thus, although the pk-turn cannot intrinsically fold into the kinked structure, it can be induced to fold correctly in context. And the pk-turn and k-turns can substitute functionally for one another.

The use of a combination of computer-assisted structure prediction and SHAPE probing to elucidate the secondary structures of five viroids

The elucidation of the structures of viroids, noncoding infectious RNA species, is paramount to obtain an understanding of the various aspects of their life cycles (including replication, transport and pathogenesis). In general, the secondary structures of viroids have been predicted using computer software programs which have been shown to possess several important limitations. Clearly, the predicted structure of a viroid needs to receive physical support prior to its use in the accurate interpretation of any mechanistic studies. Here, SHAPE probing coupled to computer-assisted structure prediction using the RNAstructure software program was employed to determine the structures of five viroids. These species belong to four genera of the Pospiviroidae family, and none have had their structure characterized in solution. In addition, several interesting questions were addressed by either studying various sequence variants or varying the SHAPE conditions. More importantly, this approach is novel in the study of viroids, and should be of significant aid in the determination of the structures of other RNA species.

A structural module in RNase P expands the variety of RNA kinks

RNA structures are built from recurrent modules that can be identified by structural and comparative sequence analysis. In order to assemble sets of helices in compact architectures, modules that introduce bends and kinks are necessary. Among such modules, kink-turns form an important family that presents sequence and structural characteristics. Here, we describe an internal loop in the bacterial type A RNase P RNA that sets helices bound at the junctions exactly in the same relative positions as in kink-turns but without the structural signatures typical of kink-turns. Our work suggests that identifying a structural module in a subset of RNA sequences constitutes a strategy to identify distinct sequential motifs sharing common structural characteristics.

Selection of the most potent specific on/off adaptor-hepatitis delta virus ribozymes for use in gene targeting

The Hepatitis Delta Virus (HDV) ribozyme, which is well adapted to the environment of the human cell, is an excellent candidate for the future development of gene-inactivation systems. On top of this, a new generation of HDV ribozymes now exists that benefits from the addition of a specific on/off adaptor (specifically the SOFA-HDV ribozymes) which greatly increases both the ribozyme's specificity and its cleavage activity. Unlike RNAi and hammerhead ribozymes, the designing of SOFA-HDV ribozymes to cleave, in trans, given RNA species has never been the object of a systematic optimization study, even with their recent use for the gene knockdown of various targets. This report aims at both improving and clarifying the design process of SOFA-HDV ribozymes. Both the ribozyme and the targeted RNA substrate were analyzed in order to provide new criteria that are useful in the selection of the most potent SOFA-HDV ribozymes. The crucial features present in both the ribozyme's biosensor and blocker, as well as at the target site, were identified and characterized. Simple rules were derived and tested using hepatitis C virus NS5B RNA as a model target. Overall, this method should promote the use of the SOFA-HDV ribozymes in a plethora of applications in both functional genomics and gene therapy.

Intermolecular interaction between a branching ribozyme and associated homing endonuclease mRNA

RNA tertiary interactions involving docking of GNRA (N; any base; R; purine) hairpin loops into helical stem structures on other regions of the same RNA are one of the most common RNA tertiary interactions. In this study, we investigated a tertiary association between a GAAA hairpin tetraloop in a small branching ribozyme (DiGIR1) and a receptor motif (HEG P1 motif) present in a hairpin structure on a separate mRNA molecule. DiGIR1 generates a 2', 5' lariat cap at the 5' end of its downstream homing endonuclease mRNA by catalysing a self-cleavage branching reaction at an internal processing site. Upon release, the 5' end of the mRNA forms a distinct hairpin structure termed HEG P1. Our biochemical data, in concert with molecular 3D modelling, provide experimental support for an intermolecular tetraloop receptor interaction between the L9 GAAA in DiGIR1 and a GNRA tetraloop receptor-like motif (UCUAAG-CAAGA) found within the HEG P1. The biological role of this interaction appears to be linked to the homing endonuclease expression by promoting post-cleavage release of the lariat capped mRNA. These findings add to our understanding of how protein-coding genes embedded in nuclear ribosomal DNA are expressed in eukaryotes and controlled by ribozymes.

Binding-interaction prediction of RNA-binding ligands

RNA molecules are involved in a wide range of biological processes and have been recognized as very important therapeutic targets. Mainly owing to the scarcity of information and experimental studies, the application of computational approaches and, in particular, of docking methodologies in the RNA field has developed slowly. However, in recent years the docking of RNA-binding ligands has experienced significant expansion. This article focuses attention on the docking of RNA-binding ligands, analyzing the development of RNA-docking approaches, the reliability of the docking methods and, finally, evaluating the results of docking-based virtual screening studies reported in the literature.