The role of specific 2'-hydroxyl groups in the stabilization of the folded conformation of kink-turn RNA

RIP-PEN-seq identifies a class of kink-turn RNAs as splicing regulators

A kink-turn (K-turn) is a three-dimensional RNA structure that exists in all three primary phylogenetic domains. In this study, we developed the RIP-PEN-seq method to identify the full-length sequences of RNAs bound by the K-turn binding protein 15.5K and discovered a previously uncharacterized class of RNAs with backward K-turn motifs (bktRNAs) in humans and mice. All bktRNAs share two consensus sequence motifs at their fixed terminal position and have complex folding properties, expression and evolution patterns. We found that a highly conserved bktRNA1 guides the methyltransferase fibrillarin to install RNA methylation of U12 small nuclear RNA in humans. Depletion of bktRNA1 causes global splicing dysregulation of U12-type introns by impairing the recruitment of ZCRB1 to the minor spliceosome. Most bktRNAs regulate the splicing of local introns by interacting with the 15.5K protein. Taken together, our findings characterize a class of small RNAs and uncover another layer of gene expression regulation that involves crosstalk among bktRNAs, RNA splicing and RNA methylation.

Structure and ion-dependent folding of k-junctions

k-Junctions are elaborated forms of kink turns with an additional helix on the nonbulged strand, thus forming a three-way helical junction. Two were originally identified in the structures of Arabidopsis and Escherichia coli thiamine pyrophosphate (TPP) riboswitches, and another called DUF-3268 was tentatively identified from sequence information. In this work we show that the Arabidopsis and E. coli riboswitch k-junctions fold in response to the addition of magnesium or sodium ions, and that atomic mutations that should disrupt key hydrogen bonding interactions greatly impair folding. Using X-ray crystallography, we have determined the structure of the DUF-3268 RNA and thus confirmed that it is a k-junction. It also folds upon the addition of metal ions, though requiring a 40-fold lower concentration of either divalent or monovalent ions. The key difference between the DUF-3268 and riboswitch k-junctions is the lack of nucleotides inserted between G1b and A2b in the former. We show that this insertion is primarily responsible for the difference in folding properties. Finally, we show that the DUF-3268 can functionally substitute for the k-junction in the E. coli TPP riboswitch such that the chimera can bind the TPP ligand, although less avidly.

Natural deep eutectic solvents protect RNA from thermal-induced degradation

RNA is a fundamental nucleic acid for life and it plays important roles in the regulation of gene transcription, post-transcriptional regulation, and epigenetic regulation. Recently, the focus on this nucleic acid has significantly increased due to the development of mRNA vaccines and RNA-based gene therapy protocols. Unfortunately, RNA based products show constrains mainly owing to instability and easy degradability of the RNA molecules. Indeed, unlike the DNA molecule which has a great intrinsic stability, RNA is more prone to degradation and this process is accelerated under thermal treatment. Here we describe a method that involves the use of Natural Deep Eutectic Solvents (NaDES) capable of slowing down RNA degradation process. Our results show that this technology seems suitable for improving the stability of specific RNA molecules particularly susceptible to thermal-induced degradation. Therefore, this technique represents a valuable tool to stabilize RNA molecules used in gene therapy and mRNA vaccines.

RNAMotifComp: a comprehensive method to analyze and identify structurally similar RNA motif families

MOTIVATION

The 3D structures of RNA play a critical role in understanding their functionalities. There exist several computational methods to study RNA 3D structures by identifying structural motifs and categorizing them into several motif families based on their structures. Although the number of such motif families is not limited, a few of them are well-studied. Out of these structural motif families, there exist several families that are visually similar or very close in structure, even with different base interactions. Alternatively, some motif families share a set of base interactions but maintain variation in their 3D formations. These similarities among different motif families, if known, can provide a better insight into the RNA 3D structural motifs as well as their characteristic functions in cell biology.

RESULTS

In this work, we proposed a method, RNAMotifComp, that analyzes the instances of well-known structural motif families and establishes a relational graph among them. We also have designed a method to visualize the relational graph where the families are shown as nodes and their similarity information is represented as edges. We validated our discovered correlations of the motif families using RNAMotifContrast. Additionally, we used a basic Naïve Bayes classifier to show the importance of RNAMotifComp. The relational analysis explains the functional analogies of divergent motif families and illustrates the situations where the motifs of disparate families are predicted to be of the same family.

AVAILABILITY AND IMPLEMENTATION

Source code publicly available at https://github.com/ucfcbb/RNAMotifFamilySimilarity.

Synthesis and properties of fully-modified 4'-selenoRNA, an endonuclease-resistant RNA analog

A large number of chemically modified oligonucleotides (ONs) have been developed for RNA-based technologies. In most modified RNAs, the characteristic 2'-hydroxyl (2'-OH) groups are removed to enhance both nuclease resistance and hybridization ability. However, the importance of the 2'-OH group in RNA structure and function is well known. Here, we report the synthesis and properties of 4'-selenoRNA in which all four nucleoside units retain the 2'-OH groups but contain a selenium atom instead of an oxygen atom at the 4'-position of the furanose ring. 4'-SelenoRNA has enhanced ability to form duplexes with RNA, and high endonuclease resistance despite the presence of the 2'-OH groups. X-ray crystallography analysis showed that the 4'-selenoRNA duplex adopts an A-conformation, similar to natural RNA, although one 4'-selenocytidine residue has unusual South-type sugar puckering. Furthermore, preliminary studies using 4'-seleno-modified siRNAs suggest that 4'-selenoRNA may be applicable to RNA interference technology. Collectively, our results raise the possibility of a new class of modified RNA in which 2'-OH groups do not need to be masked.

Structures of artificially designed discrete RNA nanoarchitectures at near-atomic resolution

Although advances in nanotechnology have enabled the construction of complex and functional synthetic nucleic acid–based nanoarchitectures, high-resolution discrete structures are lacking because of the difficulty in obtaining good diffracting crystals. Here, we report the design and construction of RNA nanostructures based on homooligomerizable one-stranded tiles for x-ray crystallographic determination. We solved three structures to near-atomic resolution: a 2D parallelogram, a 3D nanobracelet unexpectedly formed from an RNA designed for a nanocage, and, eventually, a bona fide 3D nanocage designed with the guidance of the two previous structures. Structural details of their constituent motifs, such as kissing loops, branched kissing loops, and T-junctions, that resemble natural RNA motifs and resisted x-ray determination are revealed, providing insights into those natural motifs. This work unveils the largely unexplored potential of crystallography in gaining high-resolution feedback for nanoarchitectural design and suggests a route to investigate RNA motif structures by configuring them into nanoarchitectures.

Directed evolution of orthogonal RNA-RBP pairs through library-vs-library in vitro selection

RNA-binding proteins (RBPs) and their RNA ligands play many critical roles in gene regulation and RNA processing in cells. They are also useful for various applications in cell biology and synthetic biology. However, re-engineering novel and orthogonal RNA-RBP pairs from natural components remains challenging while such synthetic RNA-RBP pairs could significantly expand the RNA-RBP toolbox for various applications. Here, we report a novel library-vs-library in vitro selection strategy based on Phage Display coupled with Systematic Evolution of Ligands by EXponential enrichment (PD-SELEX). Starting with pools of 1.1 × 1012 unique RNA sequences and 4.0 × 108 unique phage-displayed L7Ae-scaffold (LS) proteins, we selected RNA-RBP complexes through a two-step affinity purification process. After six rounds of library-vs-library selection, the selected RNAs and LS proteins were analyzed by next-generation sequencing (NGS). Further deconvolution of the enriched RNA and LS protein sequences revealed two synthetic and orthogonal RNA-RBP pairs that exhibit picomolar affinity and >4000-fold selectivity.

Structure and folding of four putative kink turns identified in structured RNA species in a test of structural prediction rules

k-Turns are widespread key architectural elements that occur in many classes of RNA molecules. We have shown previously that their folding properties (whether or not they fold into their tightly kinked structure on addition of metal ions) and conformation depend on their local sequence, and we have elucidated a series of rules for prediction of these properties from sequence. In this work, we have expanded the rules for prediction of folding properties, and then applied the full set to predict the folding and conformation of four probable k-turns we have identified amongst 224 structured RNA species found in bacterial intergenenic regions by the Breaker lab (1). We have analyzed the ion-dependence of folding of the four k-turns using fluorescence resonance energy transfer, and determined the conformation of two of them using X-ray crystallography. We find that the experimental data fully conform to both the predicted folding and conformational properties. We conclude that our folding rules are robust, and can be applied to new k-turns of unknown characteristics with confidence.

Using an L7Ae-Tethered, Hydroxyl Radical-Mediated Footprinting Strategy to Identify and Validate Kink-Turns in RNAs

Kink-turns are important RNA structural modules that facilitate long-range tertiary interactions and form binding sites for members of the L7Ae family of proteins. Present in a wide variety of functional RNAs, kink-turns play key organizational roles in many RNA-based cellular processes, including translation, modification, and tRNA biogenesis. It is important to determine the contribution of kink-turns to the overall architecture of resident RNAs, as these modules dictate ribonucleoprotein (RNP) assembly and function. This chapter describes a site-directed, hydroxyl radical-mediated footprinting strategy that utilizes L7Ae-tethered chemical nucleases to experimentally validate computationally identified kink-turns in any RNA and under a wide variety of conditions. The work plan described here uses the catalytic RNase P RNA as an example to provide a blueprint for using this footprinting method to map RNA-protein interactions in other RNP complexes.

Simulation Study of the Plasticity of k-Turn Motif in Different Environments

The k-turn is a widespread and important motif in RNA. According to the internal hydrogen bond network, it has two stable states, called N1 and N3. The relative stability between the states changes with the environment. It is able to accept different conformations in different environments. This is called the "plasticity" of a molecule. In this work, we study the plasticity of k-turn by the mixing REMD method in explicit solvent. The results are concluded as follows. First, N1 and N3 are almost equally stable when k-turn is in the solvent alone. The molecule is quite flexible as a hinge. However, after binding to different proteins, such as the proteins L7Ae and L24e, k-turn falls into one global minimum. The preferred state could be either N1 or N3. On the contrary, the other nonpreferred state becomes unstable with a weaker binding affinity to the protein. It reveals that RNA-binding protein is able to modulate the representative state of k-turn at equilibrium. This is in agreement with the findings in experiments. Moreover, free energy calculations show that the free energy barrier between the N1 and N3 states of k-turn increases in the complexes. The state-to-state transition is greatly impeded. We also give a deep discussion on the mechanism of the high plasticity of k-turn in different environments.

Effect of methylation of adenine N6 on kink turn structure depends on location

N6-methyladenine is the most common covalent modification in cellular RNA species, with demonstrated functional consequences. At the molecular level this methylation could alter local RNA structure, and/or modulate the binding of specific proteins. We have previously shown that trans-Hoogsteen-sugar (sheared) A:G base pairs can be completely disrupted by methylation, and that this occurs in a sub-set ofD/D k-turn structures. In this work we have investigated to what extent sequence context affects the severity with which inclusion of N6-methyladenine into different A:G base pairs of a standard k-turn affects RNA folding and L7Ae protein binding. We find that local sequence has a major influence, ranging from complete absence of folding and protein binding to a relatively mild effect. We have determined the crystal structure of one of these species both free and protein-bound, showing the environment of the methyl group and the way the modification is accommodated into the k-turn structure.

The kink-turn in the structural biology of RNA

The kink-turn (k-turn) is a widespread structural motif found in functional RNA species. It typically comprises a three-nucleotide bulge followed by tandem trans sugar edge-Hoogsteen G:A base pairs. It introduces a sharp kink into the axis of duplex RNA, juxtaposing the minor grooves. Cross-strand H-bonds form at the interface, accepted by the conserved adenine nucleobases of the G:A basepairs. Alternative acceptors for one of these divides the k-turns into two conformational classes N3 and N1. The base pair that follows the G:A pairs (3b:3n) determines which conformation is adopted by a given k-turn. k-turns often mediate tertiary contacts in folded RNA species and frequently bind proteins. Common k-turn binding proteins include members of the L7Ae family, such as the human 15·5k protein. A recognition helix within these proteins binds in the widened major groove on the outside of the k-turn, that makes specific H-bonds with the conserved guanine nucleobases of the G:A pairs. L7Ae binds with extremely high affinity, and single-molecule data are consistent with folding by conformational selection. The standard, simple k-turn can be elaborated in a variety of ways, that include the complex k-turns and the k-junctions. In free solution in the absence of added metal ions or protein k-turns do not adopt the tightly-kinked conformation. They undergo folding by the binding of proteins, by the formation of tertiary contacts, and some (but not all) will fold on the addition of metal ions. Whether or not folding occurs in the presence of metal ions depends on local sequence, including the 3b:3n position, and the -1b:-1n position (5' to the bulge). In most cases -1b:-1n = C:G, so that the 3b:3n position is critical since it determines both folding properties and conformation. In general, the selection of these sequence matches a given k-turn to its biological requirements. The k-turn structure is now very well understood, to the point at which they can be used as a building block for the formation of RNA nano-objects, including triangles and squares.

Methylation guide RNA evolution in archaea: structure, function and genomic organization of 110 C/D box sRNA families across six Pyrobaculum species

Archaeal homologs of eukaryotic C/D box small nucleolar RNAs (C/D box sRNAs) guide precise 2′-O-methyl modification of ribosomal and transfer RNAs. Although C/D box sRNA genes constitute one of the largest RNA gene families in archaeal thermophiles, most genomes have incomplete sRNA gene annotation because reliable, fully automated detection methods are not available. We expanded and curated a comprehensive gene set across six species of the crenarchaeal genus Pyrobaculum, particularly rich in C/D box sRNA genes. Using high-throughput small RNA sequencing, specialized computational searches and comparative genomics, we analyzed 526 Pyrobaculum C/D box sRNAs, organizing them into 110 families based on synteny and conservation of guide sequences which determine methylation targets. We examined gene duplications and rearrangements, including one family that has expanded in a pattern similar to retrotransposed repetitive elements in eukaryotes. New training data and inclusion of kink-turn secondary structural features enabled creation of an improved search model. Our analyses provide the most comprehensive, dynamic view of C/D box sRNA evolutionary history within a genus, in terms of modification function, feature plasticity, and gene mobility.

The role of RNA structure in translational regulation by L7Ae protein in archaea

A recent study has shown that archaeal L7Ae binds to a putative k-turn structure in the 5'-leader of the mRNA of its structural gene to regulate translation. To function as a regulator, the RNA should be unstructured in the absence of protein, but it should adopt a k-turn-containing stem-loop on binding L7Ae. Sequence analysis of UTR sequences indicates that their k-turn elements will be unable to fold in the absence of L7Ae, and we have demonstrated this experimentally in solution using FRET for the Archaeoglobus fulgidus sequence. We have solved the X-ray crystal structure of the complex of the A. fulgidus RNA bound to its cognate L7Ae protein. The RNA adopts a standard k-turn conformation that is specifically recognized by the L7Ae protein, so stabilizing the stem-loop. In-line probing of the natural-sequence UTR shows that the RNA is unstructured in the absence of L7Ae binding, but folds on binding the protein such that the ribosome binding site is occluded. Thus, L7Ae regulates its own translation by switching the conformation of the RNA to alter accessibility.

Crystal structures of the archaeal RNase P protein Rpp38 in complex with RNA fragments containing a K-turn motif

A characteristic feature of archaeal ribonuclease P (RNase P) RNAs is that they have extended helices P12.1 and P12.2 containing kink-turn (K-turn) motifs to which the archaeal RNase P protein Rpp38, a homologue of the human RNase P protein Rpp38, specifically binds. PhoRpp38 from the hyperthermophilic archaeon Pyrococcus horikoshii is involved in the elevation of the optimum temperature of the reconstituted RNase P by binding the K-turns in P12.1 and P12.2. Previously, the crystal structure of PhoRpp38 in complex with the K-turn in P12.2 was determined at 3.4 Å resolution. In this study, the crystal structure of PhoRpp38 in complex with the K-turn in P12.2 was improved to 2.1 Å resolution and the structure of PhoRpp38 in complex with the K-turn in P12.1 was also determined at a resolution of 3.1 Å. Both structures revealed that Lys35, Asn38 and Glu39 in PhoRpp38 interact with characteristic G·A and A·G pairs in the K-turn, while Thr37, Asp59, Lys84, Glu94, Ala96 and Ala98 in PhoRpp38 interact with the three-nucleotide bulge in the K-turn. Moreover, an extended stem-loop containing P10-P12.2 in complex with PhoRpp38, as well as PhoRpp21 and PhoRpp29, which are the archaeal homologues of the human proteins Rpp21 and Rpp29, respectively, was affinity-purified and crystallized. The crystals thus grown diffracted to a resolution of 6.35 Å. Structure determination of the crystals will demonstrate the previously proposed secondary structure of stem-loops including helices P12.1 and P12.2 and will also provide insight into the structural organization of the specificity domain in P. horikoshii RNase P RNA.

Sequence determinants of the folding properties of box C/D kink-turns in RNA

Folding properties differ markedly between kink-turns (k-turns) that have different biological functions. While ribosomal and riboswitch k-turns generally fold into their kinked conformation on addition of metal ions, box C/D snoRNP k-turns remain completely unfolded under these conditions, although they fold on addition of L7Ae protein. Sequence elements have been systematically exchanged between a standard ribosomal k-turn (Kt-7) that folds on addition of metal ions, and a box C/D k-turn. Folding was studied using fluorescence resonance energy transfer and gel electrophoresis. Three sequence elements each contribute in an approximately additive manner to the different folding properties of Kt-7 and box C/D k-turns from archaea. Bioinformatic analysis indicates that k-turn sequences evolve sequences that suit their folding properties to their biological function. The majority of ribosomal and riboswitch k-turns have sequences allowing unassisted folding in response to the presence of metal ions. In contrast, box C/D k-turns have sequences that require the binding of proteins to drive folding into the kinked conformation, consistent with their role in the assembly of the box C/D snoRNP apparatus. The rules governing the influence of sequence on folding properties can be applied to other standard k-turns to predict their folding characteristics.

Role of hydroxyl groups in the B-ring of flavonoids in stabilization of the Hoogsteen paired third strand of Poly(U).Poly(A)*Poly(U) triplex

We have reported the interaction of two flavonoids namely quercetin (Q) and morin (M) with double stranded poly(A).poly(U) (herein after A.U) and triple stranded poly(U).poly(A)*poly(U) (herein after U.A*U, dot represents the Watson-Crick and asterisk represents Hoogsteen base pairing respectively) in this article. It has been observed that relative positions of hydroxyl groups on the B-ring of the flavonoids affect the stabilization of RNA. The double strand as well as the triple strand of RNA-polymers become more stabilized in presence of Q, however both the duplex and triplex remain unaffected in presence of M. The presence of catechol moiety on the B-ring of Q is supposed to be responsible for the stabilization. Moreover, after exploiting a series of biophysical experiments, it has been found that, triple helical RNA becomes more stabilized over its parent duplex in presence of Q. Fluorescence quenching, viscosity measurement and helix melting results establish the fact that Q binds with both forms of RNA through the mode of intercalation while M does not bind at all to either forms of RNA.

A quasi-cyclic RNA nano-scale molecular object constructed using kink turns

k-Turns are widespread RNA architectural elements that mediate tertiary interactions. We describe a double-kink-turn motif comprising two inverted k-turns that forms a tight horse-shoe structure that can assemble into a variety of shapes by coaxial association of helical ends. Using X-ray crystallography we show that these assemble with two (dumbell), three (triangle) and four units (square), with or without bound protein, within the crystal lattice. In addition, exchange of a single basepair can almost double the pore radius or shape of a molecular assembly. On the basis of this analysis we synthesized a 114 nt self-complementary RNA containing six k-turns. The crystal structure of this species shows that it forms a quasi-cyclic triangular object. These are randomly disposed about the three-fold axis in the crystal lattice, generating a circular RNA of quasi D3 symmetry with a shape reminiscent of that of a cyclohexane molecule in its chair conformation. This work demonstrates that the k-turn is a powerful building block in the construction of nano-scale molecular objects, and illustrates why k-turns are widely used in natural RNA molecules to organize long-range architecture and mediate tertiary contacts.

A critical base pair in k-turns determines the conformational class adopted, and correlates with biological function

k-turns are commonly-occurring motifs that introduce sharp kinks into duplex RNA, thereby facilitating tertiary contacts. Both the folding and conformation of k-turns are determined by their local sequence. k-turns fall into two conformational classes, called N3 and N1, that differ in the pattern of hydrogen bonding in the core. We show here that this is determined by the basepair adjacent to the critical G•A pairs. We determined crystal structures of a series of Kt-7 variants in which this 3b,3n position has been systematically varied, showing that this leads to a switch in the conformation. We have previously shown that the 3b,3n position also determines the folding characteristics of the k-turn, i.e. whether or not the k-turn can fold in the presence of metal ions alone. We have analyzed the distribution of 3b,3n sequences from four classes of k-turns from ribosomes, riboswitches and U4 snRNA, finding a strong conservation of properties for a given k-turn type. We thus demonstrate a strong association between biological function, 3b,3n sequence and k-turn folding and conformation. This has strong predictive power, and can be applied to the modeling of large RNA architectures.

Structural basis for recognition of a kink-turn motif by an archaeal homologue of human RNase P protein Rpp38

PhoRpp38 in the hyperthermophilic archaeon Pyrococcus horikoshii, a homologue of human ribonuclease P (RNase P) protein Rpp38, belongs to the ribosomal protein L7Ae family that specifically recognizes a kink-turn (K-turn) motif. A previous biochemical study showed that PhoRpp38 specifically binds to two stem-loops, SL12 and SL16, containing helices P12.1/12.2 and P15/16 respectively, in P. horikoshii RNase P RNA (PhopRNA). In order to gain insight into the PhoRpp38 binding mode to PhopRNA, we determined the crystal structure of PhoRpp38 in complex with the SL12 mutant (SL12M) at a resolution of 3.4 Å. The structure revealed that Lys35 on the β-strand (β1) and Asn38, Glu39, and Lys42 on the α-helix (α2) in PhoRpp38 interact with characteristic G•A and A•G pairs in SL12M, where Ile93, Glu94, and Val95, on a loop between α4 and β4 in PhoRpp38, interact with the 3-nucleotide bulge (G-G-U) in the SL12M. The structure demonstrates the previously proposed secondary structure of SL12, including helix P12.2. Structure-based mutational analysis indicated that amino acid residues involved in the binding to SL12 are also responsible for the binding to SL16. This result suggested that each PhoRpp38 binds to the K-turns in SL12 and SL16 in PhopRNA. A pull-down assay further suggested the presence of a second K-turn in SL12. Based on the present results, together with available data, we discuss a structural basis for recognition of K-turn motifs in PhopRNA by PhoRpp38.

The Kink Turn, a Key Architectural Element in RNA Structure

Kink turns (k-turns) are widespread structural elements that introduce an axial bend into duplex RNA with an included angle of 50°. These mediate key tertiary interactions and bind specific proteins including members of the L7Ae family. The standard k-turn comprises a three-nucleotide bulge followed by G·A and A·G pairs. The RNA kinks by an association of the two minor grooves, stabilized by the formation of a number of key cross-strand hydrogen bonds mostly involving the adenine bases of the G·A and A·G pairs. The k-turns may be divided into two conformational classes, depending on the receptor for one of these hydrogen bonds. k-turns become folded by one of three different processes. Some, but not all, k-turns become folded in the presence of metal ions. Whether or not a given k-turn is folded under these conditions is determined by its sequence. We present a set of rules for the prediction of folding properties and the structure adopted on local sequence.

A Highly Coupled Network of Tertiary Interactions in the SAM-I Riboswitch and Their Role in Regulatory Tuning

RNA folding in vivo is significantly influenced by transcription, which is not necessarily recapitulated by Mg(2+)-induced folding of the corresponding full-length RNA in vitro. Riboswitches that regulate gene expression at the transcriptional level are an ideal system for investigating this aspect of RNA folding as ligand-dependent termination is obligatorily co-transcriptional, providing a clear readout of the folding outcome. The folding of representative members of the SAM-I family of riboswitches has been extensively analyzed using approaches focusing almost exclusively upon Mg(2+) and/or S-adenosylmethionine (SAM)-induced folding of full-length transcripts of the ligand binding domain. To relate these findings to co-transcriptional regulatory activity, we have investigated a set of structure-guided mutations of conserved tertiary architectural elements of the ligand binding domain using an in vitro single-turnover transcriptional termination assay, complemented with phylogenetic analysis and isothermal titration calorimetry data. This analysis revealed a conserved internal loop adjacent to the SAM binding site that significantly affects ligand binding and regulatory activity. Conversely, most single point mutations throughout key conserved features in peripheral tertiary architecture supporting the SAM binding pocket have relatively little impact on riboswitch activity. Instead, a secondary structural element in the peripheral subdomain appears to be the key determinant in observed differences in regulatory properties across the SAM-I family. These data reveal a highly coupled network of tertiary interactions that promote high-fidelity co-transcriptional folding of the riboswitch but are only indirectly linked to regulatory tuning.

A critical base pair in k-turns that confers folding characteristics and correlates with biological function

Kink turns (k-turns) are widespread elements in RNA that mediate tertiary contacts by kinking the helical axis. We have found that the ability of k-turns to undergo ion-induced folding is conferred by a single base pair that follows the conserved A·G pairs, that is, the 3b·3n position. A Watson–Crick pair leads to an inability to fold in metal ions alone, while 3n=G or 3b=C (but not both) permits folding. Crystallographic study reveals two hydrated metal ions coordinated to O6 of G3n and G2n of Kt-7. Removal of either atom impairs Mg²⁺-induced folding in solution. While SAM-I riboswitches have 3b·3n sequences that would predispose them to ion-induced folding, U4 snRNA are strongly biased to an inability to such folding. Thus riboswitch sequences allow folding to occur independently of protein binding, while U4 should remain unfolded until bound by protein. The empirical rules deduced for k-turn folding have strong predictive value.

The k-turn is a widespread RNA element that adopts a kinked structure that mediates tertiary contacts and frequently binds specific proteins. Here, McPhee et al. show that the ability of a given k-turn to fold in the presence of metal ions alone—or to otherwise require protein binding—is attributable to a specific base pair.

The K-turn motif in riboswitches and other RNA species

The kink turn is a widespread structure motif that introduces a tight bend into the axis of duplex RNA. This generally functions to mediate tertiary interactions, and to serve as a specific protein binding site. K-turns or closely related structures are found in at least seven different riboswitch structures, where they function as key architectural elements that help generate the ligand binding pocket. This article is part of a Special Issue entitled: Riboswitches.

The k-junction motif in RNA structure

The k-junction is a structural motif in RNA comprising a three-way helical junction based upon kink turn (k-turn) architecture. A computer program written to examine relative helical orientation identified the three-way junction of the Arabidopsis TPP riboswitch as an elaborated k-turn. The Escherichia coli TPP riboswitch contains a related k-junction, and analysis of >11 000 sequences shows that the structure is common to these riboswitches. The k-junction exhibits all the key features of an N1-class k-turn, including the standard cross-strand hydrogen bonds. The third helix of the junction is coaxially aligned with the C (canonical) helix, while the k-turn loop forms the turn into the NC (non-canonical) helix. Analysis of ligand binding by ITC and global folding by gel electrophoresis demonstrates the importance of the k-turn nucleotides. Clearly the basic elements of k-turn structure are structurally well suited to generate a three-way helical junction, retaining all the key features and interactions of the k-turn.

Structure of a rare non-standard sequence k-turn bound by L7Ae protein

Kt-23 from Thelohania solenopsae is a rare RNA kink turn (k-turn) where an adenine replaces the normal guanine at the 2n position. L7Ae is a member of a strongly conserved family of proteins that bind a range of k-turn structures in the ribosome, box C/D and H/ACA small nucleolar RNAs and U4 small nuclear RNA. We have solved the crystal structure of T. solenopsae Kt-23 RNA bound to Archeoglobus fulgidus L7Ae protein at a resolution of 2.95 Å. The protein binds in the major groove displayed on the outer face of the k-turn, in a manner similar to complexes with standard k-turn structures. The k-turn adopts a standard N3 class conformation, with a single hydrogen bond from A2b N6 to A2n N3. This contrasts with the structure of the same sequence located in the SAM-I riboswitch, where it adopts an N1 structure, showing the inherent plasticity of k-turn structure. This potentially can affect any tertiary interactions in which the RNA participates.

The molecular recognition of kink-turn structure by the L7Ae class of proteins

L7Ae is a member of a protein family that binds kink-turns (k-turns) in many functional RNA species. We have solved the X-ray crystal structure of the near-consensus sequence Kt-7 of Haloarcula marismortui bound by Archaeoglobus fulgidus L7Ae at 2.3-Å resolution. We also present a structure of Kt-7 in the absence of bound protein at 2.2-Å resolution. As a result, we can describe a general mode of recognition of k-turn structure by the L7Ae family proteins. The protein makes interactions in the widened major groove on the outer face of the k-turn. Two regions of the protein are involved. One is an α-helix that enters the major groove of the NC helix, making both nonspecific backbone interactions and specific interactions with the guanine nucleobases of the conserved G • A pairs. A hydrophobic loop makes close contact with the L1 and L2 bases, and a glutamate side chain hydrogen bonds with L1. Taken together, these interactions are highly selective for the structure of the k-turn and suggest how conformational selection of the folded k-turn occurs.

Propensities for loop structures of RNA & DNA backbones

RNA oligonucleotides exhibit a large tendency to bend and form a loop conformation which is a major motif contributing to their complex three-dimensional structure. This is in contrast to DNA molecules that predominantly form the double-helix structure. In this paper we investigate by molecular dynamics simulation, as well as, by its combination with the replica-exchange method, the propensity of RNA chains containing the GCUAA pentaloop to form spontaneously a hairpin conformation. The results were then compared with those of analogous hybrid oligonucleotides in which the ribose groups in the loop-region were substituted by deoxyriboses. We find that the RNA oligomers exhibit a marginal excess stability to form loop structures. The equilibrium constant for opening the loop to an extended conformation is twice as large in the hybrid than it is in the RNA chain. Analyses of the hydrogen bonds indicate that the excess stability for forming a hairpin is a result of hydrogen bonds the 2'-hydroxyls in the loop region form with other groups in the loop. Of these hydrogen bonds, the most important is the hydrogen bond donated from the 2'-OH at the first position of the loop to N7 of adenine at the forth position. RNA and DNA backbones are characterized by different backbone dihedral angles and sugar puckering that can potentially facilitate or hamper the hydrogen bonds involving the 2'-OH. Nevertheless, the sugar puckerings of all the pentaloop nucleotides were not significantly different between the two chains displaying the C3'-endo conformation characteristic to the A-form double helix. All of the other backbone dihedrals also did not show any considerable difference in the loop-region except of the δ-dihedral. In this case, the RNA loop exhibited bimodal distributions corresponding to, both, the RNA and DNA backbones, whereas the loop of the hybrid chain behaved mostly as that of a DNA backbone. Thus, it is possible that the behavior of the δ-dihedrals in the loop-region of the RNA adopts conformations that facilitate the intra-nucleotide hydrogen bondings of the 2'-hydroxyls, and consequently renders loop structures in RNA more stable.

Bioinformatics and molecular dynamics simulation study of L1 stalk non-canonical rRNA elements: kink-turns, loops, and tetraloops

The L1 stalk is a prominent mobile element of the large ribosomal subunit. We explore the structure and dynamics of its non-canonical rRNA elements, which include two kink-turns, an internal loop, and a tetraloop. We use bioinformatics to identify the L1 stalk RNA conservation patterns and carry out over 11.5 μs of MD simulations for a set of systems ranging from isolated RNA building blocks up to complexes of L1 stalk rRNA with the L1 protein and tRNA fragment. We show that the L1 stalk tetraloop has an unusual GNNA or UNNG conservation pattern deviating from major GNRA and YNMG RNA tetraloop families. We suggest that this deviation is related to a highly conserved tertiary contact within the L1 stalk. The available X-ray structures contain only UCCG tetraloops which in addition differ in orientation (anti vs syn) of the guanine. Our analysis suggests that the anti orientation might be a mis-refinement, although even the anti interaction would be compatible with the sequence pattern and observed tertiary interaction. Alternatively, the anti conformation may be a real substate whose population could be pH-dependent, since the guanine syn orientation requires protonation of cytosine in the tertiary contact. In absence of structural data, we use molecular modeling to explore the GCCA tetraloop that is dominant in bacteria and suggest that the GCCA tetraloop is structurally similar to the YNMG tetraloop. Kink-turn Kt-77 is unusual due to its 11-nucleotide bulge. The simulations indicate that the long bulge is a stalk-specific eight-nucleotide insertion into consensual kink-turn only subtly modifying its structural dynamics. We discuss a possible evolutionary role of helix H78 and a mechanism of L1 stalk interaction with tRNA. We also assess the simulation methodology. The simulations provide a good description of the studied systems with the latest bsc0χOL3 force field showing improved performance. Still, even bsc0χOL3 is unable to fully stabilize an essential sugar-edge H-bond between the bulge and non-canonical stem of the kink-turn. Inclusion of Mg(2+) ions may deteriorate the simulations. On the other hand, monovalent ions can in simulations readily occupy experimental Mg(2+) binding sites.

The plasticity of a structural motif in RNA: structural polymorphism of a kink turn as a function of its environment

The k-turn is a widespread structural motif that introduces a tight kink into the helical axis of double-stranded RNA. The adenine bases of consecutive G•A pairs are directed toward the minor groove of the opposing helix, hydrogen bonding in a typical A-minor interaction. We show here that the available structures of k-turns divide into two classes, depending on whether N3 or N1 of the adenine at the 2b position accepts a hydrogen bond from the O2' at the -1n position. There is a coordinated structural change involving a number of hydrogen bonds between the two classes. We show here that Kt-7 can adopt either the N3 or N1 structures depending on environment. While it has the N1 structure in the ribosome, on engineering it into the SAM-I riboswitch, it changes to the N3 structure, resulting in a significant alteration in the trajectory of the helical arms.

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.

Single-molecule observation of the induction of k-turn RNA structure on binding L7Ae protein

The k-turn is a commonly occurring structural motif that introduces a tight kink into duplex RNA. In free solution, it can exist in an extended form, or by folding into the kinked structure. Binding of proteins including the L7Ae family can induce the formation of the kinked geometry, raising the question of whether this occurs by passive selection of the kinked structure, or a more active process in which the protein manipulates the RNA structure. We have devised a single-molecule experiment whereby immobilized L7Ae protein binds Cy3-Cy5-labeled RNA from free solution. We find that all bound RNA is in the kinked geometry, with no evidence for transitions to an extended form at the millisecond timescale of the camera. Furthermore, real-time binding experiments provide no evidence for a more extended intermediate even at the earliest times, at a time resolution of 16 ms. The data support a passive conformational selection model by which the protein selects a fraction of RNA that is already in the kinked conformation, thereby drawing the equilibrium into this form.

The structure and folding of kink turns in RNA

The kink turn (k-turn) is a widespread structural motif that introduces a tight kink into the axis of double-stranded RNA, with an included angle ∼60°. A standard k-turn comprises a three-nucleotide bulge followed on the 3' side by a G•A pair, an A•G pair, and usually further non-Watson-Crick pairs. The kinked conformation may be stabilized by three processes. These are the addition of metal ions, the binding of proteins such as the L7Ae family, and by the formation of tertiary interactions. The structure is characterized by specific A-minor interactions with the adenine nucleobases of the G•A pairs, and some very well-conserved hydrogen bonds involving 2'-hydroxyl groups. We can identify two classes of k-turns, that differ in the manner of the hydrogen bonding at the adenine of the bulge-distal G•A pair.

Structure and folding of a rare, natural kink turn in RNA with an A*A pair at the 2b*2n position

The kink turn (k-turn) is a frequently occurring motif, comprising a bulge followed by G•A and A•G pairs that introduces a sharp axial bend in duplex RNA. Natural k-turn sequences exhibit significant departures from the consensus, including the A•G pairs that form critical interactions stabilizing the core of the structure. Kt-23 found in the small ribosomal subunit differs from the consensus in many organisms, particularly in the second A•G pair distal to the bulge (2b•2n). Analysis of many Kt-23 sequences shows that the frequency of occurrence at the 2n position (i.e., on the nonbulged strand, normally G in standard k-turns) is U>C>G>A. Less than 1% of sequences have A at the 2n position, but one such example occurs in Thelohania solenopsae Kt-23. This sequence folds only weakly in the presence of Mg²⁺ ions but is induced to fold normally by the binding of L7Ae protein. Introduction of this sequence into the SAM-I riboswitch resulted in normal binding of SAM ligand, indicating that tertiary RNA contacts have resulted in k-turn folding. X-ray crystallography shows that the T. solenopsae Kt-23 adopts a standard k-turn geometry, making the key, conserved hydrogen bonds in the core and orienting the 1n (of the bulge-proximal A•G pair) and 2b adenine nucleobases in position facing the opposing minor groove. The 2b and 2n adenine nucleobases are not directly hydrogen bonded, but each makes hydrogen bonds to their opposing strands.

YbxF and YlxQ are bacterial homologs of L7Ae and bind K-turns but not K-loops

The archaeal protein L7Ae and eukaryotic homologs such as L30e and 15.5kD comprise the best characterized family of K-turn-binding proteins. K-turns are an RNA motif comprised of a bulge flanked by canonical and noncanonical helices. They are widespread in cellular RNAs, including bacterial gene-regulatory RNAs such as the c-di-GMP-II, lysine, and SAM-I riboswitches, and the T-box. The existence in bacteria of K-turn-binding proteins of the L7Ae family has not been proven, although two hypothetical proteins, YbxF and YlxQ, have been proposed to be L7Ae homologs based on sequence conservation. Using purified, recombinant proteins, we show that Bacillus subtilis YbxF and YlxQ bind K-turns (K(d) ~270 nM and ~2300 nM, respectively). Crystallographic structure determination demonstrates that both YbxF and YlxQ adopt the same overall fold as L7Ae. Unlike the latter, neither bacterial protein recognizes K-loops, a structural motif that lacks the canonical helix of the K-turn. This property is shared between the bacterial and eukaryal family members. Comparison of our structure of YbxF in complex with the K-turn of the SAM-I riboswitch and previously determined structures of archaeal and eukaryal homologs bound to RNA indicates that L7Ae approaches the K-turn at a unique angle, which results in a considerably larger RNA-protein interface dominated by interactions with the noncanonical helix of the K-turn. Thus, the inability of the bacterial and eukaryal L7Ae homologs to bind K-loops probably results from their reliance on interactions with the canonical helix. The biological functions of YbxF and YlxQ remain to be determined.

An energetically beneficial leader-linker interaction abolishes ligand-binding cooperativity in glycine riboswitches

Comprised of two aptamers connected by a short nucleotide linker, the glycine riboswitch was the first example of naturally occurring RNA elements reported to bind small organic molecules cooperatively. Earlier works have shown binding of glycine to the second aptamer allows tertiary interactions to be made between the two aptamers, which facilitates binding of a separate glycine molecule to the first aptamer, leading to glycine-binding cooperativity. Prompted by a distinctive protection pattern in the linker region of a minimal glycine riboswitch construct, we have identified a highly conserved (>90%) leader-linker duplex involving leader nucleotides upstream of the previously reported consensus glycine riboswitch sequences. In >50% of the glycine riboswitches, the leader-linker interaction forms a kink-turn motif. Characterization of three glycine ribsowitches showed that the leader-linker interaction improved the glycine-binding affinities by 4.5- to 86-fold. In-line probing and native gel assays with two aptamers in trans suggested synergistic action between glycine-binding and interaptamer interaction during global folding of the glycine riboswitch. Mutational analysis showed that there appeared to be no ligand-binding cooperativity in the glycine riboswitch when the leader-linker interaction is present, and the previously measured cooperativity is simply an artifact of a truncated construct missing the leader sequence.

RNA tertiary interactions in a riboswitch stabilize the structure of a kink turn

The kink turn is a widespread RNA motif that introduces an acute kink into the axis of duplex RNA, typically comprising a bulge followed by a G⋅A and A⋅G pairs. The kinked conformation is stabilized by metal ions, or the binding of proteins including L7Ae. We now demonstrate a third mechanism for the stabilization of k-turn structure, involving tertiary interactions within a larger RNA structure. The SAM-I riboswitch contains an essential standard k-turn sequence that kinks a helix so that its terminal loop can make a long-range interaction. We find that some sequence variations in the k-turn within the riboswitch do not prevent SAM binding, despite preventing the folding of the k-turn in isolation. Furthermore, two crystal structures show that the sequence-variant k-turns are conventionally folded within the riboswitch. This study shows that the folded structure of the k-turn can be stabilized by tertiary interactions within a larger RNA structure.

A-minor tertiary interactions in RNA kink-turns. Molecular dynamics and quantum chemical analysis

The RNA kink-turn is an important recurrent RNA motif, an internal loop with characteristic consensus sequence forming highly conserved three-dimensional structure. Functional arrangement of RNA kink-turns shows a sharp bend in the phosphodiester backbone. Among other signature interactions, kink-turns form A-minor interaction between their two stems. Most kink-turns possess extended A-minor I (A-I) interaction where adenine of the second A•G base pair of the NC-stem interacts with the first canonical pair of the C-stem (i.e., the receptor pair) via trans-sugar-edge/sugar-edge (tSS) and cis-sugar-edge/sugar-edge (cSS) interactions. The remaining kink-turns have less compact A-minor 0 (A-0) interaction with just one tSS contact. We show that kink-turns with A-I in ribosomal X-ray structures keep G═C receptor base pair during evolution while the inverted pair (C═G) is not realized. In contrast, kink-turns with A-0 in the observed structures alternate G═C and C═G base pairs in sequences. We carried out an extended set (~5 μs) of explicit-solvent molecular dynamics simulations of kink-turns to rationalize this structural/evolutionary pattern. The simulations were done using a net-neutral Na(+) cation atmosphere (with ~0.25 M cation concentration) supplemented by simulations with either excess salt KCl atmosphere or inclusion of Mg(2+). The results do not seem to depend on the treatment of ions. The simulations started with X-ray structures of several kink-turns while we tested the response of the simulated system to base substitutions, modest structural perturbations and constraints. The trends seen in the simulations reveal that the A-I/G═C arrangement is preferred over all three other structures. The A-I/C═G triple appears structurally entirely unstable, consistent with the covariation patterns seen during the evolution. The A-0 arrangements tend to shift toward the A-I pattern in simulations, which suggests that formation of the A-0 interaction is likely supported by the surrounding protein and RNA molecules. A-0 may also be stabilized by additional kink-turn nucleotides not belonging to the kink-turn consensus, as shown for the kink-turn from ribosomal Helix 15. Quantum-chemical calculations on all four A-minor triples suggest that there is a different balance of electrostatic and dispersion stabilization in the A-I/G═C and A-I/C═G triples, which may explain different behavior of these otherwise isosteric triples in the context of kink-turns.

Understanding RNA Flexibility Using Explicit Solvent Simulations: The Ribosomal and Group I Intron Reverse Kink-Turn Motifs

Reverse kink-turn is a recurrent elbow-like RNA building block occurring in the ribosome and in the group I intron. Its sequence signature almost matches that of the conventional kink-turn. However, the reverse and conventional kink-turns have opposite directions of bending. The reverse kink-turn lacks basically any tertiary interaction between its stems. We report unrestrained, explicit solvent molecular dynamics simulations of ribosomal and intron reverse kink-turns (54 simulations with 7.4 μs of data in total) with different variants (ff94, ff99, ff99bsc0, ff99χOL, and ff99bsc0χOL) of the Cornell et al. force field. We test several ion conditions and two water models. The simulations characterize the directional intrinsic flexibility of reverse kink-turns pertinent to their folded functional geometries. The reverse kink-turns are the most flexible RNA motifs studied so far by explicit solvent simulations which are capable at the present simulation time scale to spontaneously and reversibly sample a wide range of geometries from tightly kinked ones through flexible intermediates up to extended, unkinked structures. A possible biochemical role of the flexibility is discussed. Among the tested force fields, the latest χOL variant is essential to obtaining stable trajectories while all force field versions lacking the χ correction are prone to a swift degradation toward senseless ladder-like structures of stems, characterized by high-anti glycosidic torsions. The type of explicit water model affects the simulations considerably more than concentration and the type of ions.

Structural dynamics of the box C/D RNA kink-turn and its complex with proteins: the role of the A-minor 0 interaction, long-residency water bridges, and structural ion-binding sites revealed by molecular simulations

Kink-turns (K-turns) are recurrent elbow-like RNA motifs that participate in protein-assisted RNA folding and contribute to RNA dynamics. We carried out a set of molecular dynamics (MD) simulations using parm99 and parmbsc0 force fields to investigate structural dynamics of the box C/D RNA and its complexes with two proteins: native archaeal L7ae protein and human 15.5 kDa protein, originally bound to very similar structure of U4 snRNA. The box C/D RNA forms K-turn with A-minor 0 tertiary interaction between its canonical (C) and noncanonical (NC) stems. The local K-turn architecture is thus different from the previously studied ribosomal K-turns 38 and 42 having A-minor I interaction. The simulations reveal visible structural dynamics of this tertiary interaction involving altogether six substates which substantially contribute to the elbow-like flexibility of the K-turn. The interaction can even temporarily shift to the A-minor I type pattern; however, this is associated with distortion of the G/A base pair in the NC-stem of the K-turn. The simulations show reduction of the K-turn flexibility upon protein binding. The protein interacts with the apex of the K-turn and with the NC-stem. The protein-RNA interface includes long-residency hydration sites. We have also found long-residency hydration sites and major ion-binding sites associated with the K-turn itself. The overall topology of the K-turn remains stable in all simulations. However, in simulations of free K-turn, we observed instability of the key C16(O2')-A7(N1) H-bond, which is a signature interaction of K-turns and which was visibly more stable in simulations of K-turns possessing A-minor I interaction. It may reflect either some imbalance of the force field or it may be a correct indication of early stages of unfolding since this K-turn requires protein binding for its stabilization. Interestingly, the 16(O2')-7(N1) H- bond is usually not fully lost since it is replaced by a water bridge with a tightly bound water, which is adenine-specific similarly as the original interaction. The 16(O2')-7(N1) H-bond is stabilized by protein binding. The stabilizing effect is more visible with the human 15.5 kDa protein, which is attributed to valine to arginine substitution in the binding site. The behavior of the A-minor interaction is force-field-dependent because the parmbsc0 force field attenuates the A-minor fluctuations compared to parm99 simulations. Behavior of other regions of the box C/D RNA is not sensitive to the force field choice. Simulation with net-neutralizing Na(+) and 0.2 M excess salt conditions appear in all aspects equivalent. The simulations show loss of a hairpin tetraloop, which is not part of the K-turn. This was attributed to force field limitations.

A structural database for k-turn motifs in RNA

The kink-turn (k-turn) is a common structural motif in RNA that introduces a tight kink into the helical axis. k-turns play an important architectural role in RNA structures and serve as binding sites for a number of proteins. We have created a database of known and postulated k-turn sequences and three-dimensional (3D) structures, available via the internet. This site provides (1) a database of sequence and structure, as a resource for the RNA community, and (2) a tool to enable the manipulation and comparison of 3D structures where known.

Plasticity of the RNA kink turn structural motif

The kink turn (K-turn) is an RNA structural motif found in many biologically significant RNAs. While most examples of the K-turn have a similar fold, the crystal structure of the Azoarcus group I intron revealed a novel RNA conformation, a reverse kink turn bent in the direction opposite that of a consensus K-turn. The reverse K-turn is bent toward the major grooves rather than the minor grooves of the flanking helices, yet the sequence differs from the K-turn consensus by only a single nucleotide. Here we demonstrate that the reverse bend direction is not solely defined by internal sequence elements, but is instead affected by structural elements external to the K-turn. It bends toward the major groove under the direction of a tetraloop-tetraloop receptor. The ability of one sequence to form two distinct structures demonstrates the inherent plasticity of the K-turn sequence. Such plasticity suggests that the K-turn is not a primary element in RNA folding, but instead is shaped by other structural elements within the RNA or ribonucleoprotein assembly.

Evaluation of parameters critical for observing nucleic acids inside living Xenopus laevis oocytes by in-cell NMR spectroscopy

In-cell NMR spectroscopy of proteins in different cellular environments is a well-established technique that, however, has not been applied to nucleic acids so far. Here, we show that isotopically labeled DNA and RNA can be observed inside the eukaryotic environment of Xenopus laevis oocytes by in-cell NMR spectroscopy. One limiting factor for the observation of nucleic acids in Xenopus oocytes is their reduced stability. We demonstrate that chemical modification of DNA and RNA can protect them from degradation and can significantly enhance their lifetime. Finally, we show that the imino region of the NMR spectrum is devoid of any oocyte background signals enabling the detection even of isotopically nonlabeled molecules.

Ion-induced folding of a kink turn that departs from the conventional sequence

Kink turns (k-turns) are important structural motifs that create a sharp axial bend in RNA. Most conform to a consensus in which a three-nucleotide bulge is followed by consecutive G*A and A*G base pairs, and when these G*A pairs are modified in vitro this generally leads to a failure to adopt the k-turn conformation. Kt-23 in the 30S ribosomal subunit of Thermus thermophilus is a rare exception in which the bulge-distal A*G pair is replaced by a non-Watson-Crick A*U pair. In the context of the ribosome, Kt-23 adopts a completely conventional k-turn geometry. We show here that this sequence is induced to fold into a k-turn structure in an isolated RNA duplex by Mg(2+) or Na(+) ions. Therefore, the Kt-23 is intrinsically stable despite lacking the key A*G pair; its formation requires neither tertiary interactions nor protein binding. Moreover, the Kt-23 k-turn is stabilized by the same critical hydrogen-bonding interactions within the core of the structure that are found in more conventional sequences such as the near-consensus Kt-7. T. thermophilus Kt-23 has two further non-Watson-Crick base pairs within the non-canonical helix, three and four nucleotides from the bulge, and we find that the nature of these pairs influences the ability of the RNA to adopt k-turn conformation, although the base pair adjacent to the A*U pair is more important than the other.

Discrimination between closely related cellular metabolites by the SAM-I riboswitch

The SAM-I riboswitch is a cis-acting element of genetic control found in bacterial mRNAs that specifically binds S-adenosylmethionine (SAM). We previously determined the 2.9-A X-ray crystal structure of the effector-binding domain of this RNA element, revealing details of RNA-ligand recognition. To improve this structure, variations were made to the RNA sequence to alter lattice contacts, resulting in a 0.5-A improvement in crystallographic resolution and allowing for a more accurate refinement of the crystallographic model. The basis for SAM specificity was addressed by a structural analysis of the RNA complexed to S-adenosylhomocysteine (SAH) and sinefungin and by measuring the affinity of SAM and SAH for a series of mutants using isothermal titration calorimetry. These data illustrate the importance of two universally conserved base pairs in the RNA that form electrostatic interactions with the positively charged sulfonium group of SAM, thereby providing a basis for discrimination between SAM and SAH.

Elbow flexibility of the kt38 RNA kink-turn motif investigated by free-energy molecular dynamics simulations

Kink-turns (K-turns) are common structural motifs that can introduce sharp kinks into double-stranded RNA, and have been proposed to mediate large-scale motions in the ribosome. K-turns consist of a bulge loop region flanked by trans sugar-Hoogsteen G:A pairs, and the sharp kink conformation is stabilized by A-minor interactions (adenine contacting a G:C basepair in the minor groove). Umbrella-sampling molecular dynamics simulations were used to disrupt an A-minor interaction in the ribosomal kt38 turn and to calculate the associated free-energy change. Coupling of umbrella sampling with replica exchanges between neighboring umbrella-sampling intervals could further improve the convergence of the free-energy calculations. The simulations revealed a coupled A-minor disruption and global opening of the K-turn motif, and allowed us to characterize several intermediate A-minor conformations. The calculated free-energy profile indicated a meta-stable, semi-open structure of slightly higher free energy ( approximately 1 kcal mol(-1)), separated by a small free-energy barrier ( approximately 1.5 kcal mol(-1)) from the closed (highly kinked) form. Both K-turn states are stabilized by distinct variants of the A-minor interaction. Further opening of the K-turn structure required significantly larger free-energy changes. The semi-open form had a reduced kink angle compatible with experimental data on K-turn solution structures, and opening was coupled to a continuous global unwinding of the K-turn motif. The range of free-energy changes associated with kt38 opening and unwinding are compatible with the idea that K-turns may facilitate biologically relevant motions during large-scale ribosome dynamics.

Dynamics of the base of ribosomal A-site finger revealed by molecular dynamics simulations and Cryo-EM

Helix 38 (H38) of the large ribosomal subunit, with a length of 110 A, reaches the small subunit through intersubunit bridge B1a. Previous cryo-EM studies revealed that the tip of H38 moves by more than 10 A from the non-ratcheted to the ratcheted state of the ribosome while mutational studies implicated a key role of flexible H38 in attenuation of translocation and in dynamical signaling between ribosomal functional centers. We investigate a region including the elbow-shaped kink-turn (Kt-38) in the Haloarcula marismortui archaeal ribosome, and equivalently positioned elbows in three eubacterial species, located at the H38 base. We performed explicit solvent molecular dynamics simulations on the H38 elbows in all four species. They are formed by at first sight unrelated sequences resulting in diverse base interactions but built with the same overall topology, as shown by X-ray crystallography. The elbows display similar fluctuations and intrinsic flexibilities in simulations indicating that the eubacterial H38 elbows are structural and dynamical analogs of archaeal Kt-38. We suggest that this structural element plays a pivotal role in the large motions of H38 and may act as fulcrum for the abovementioned tip motion. The directional flexibility inferred from simulations correlates well with the cryo-EM results.

Comparative gel electrophoresis analysis of helical junctions in RNA

Comparative gel electrophoresis provides information on the relative angles subtended between helical arms at a branchpoint in RNA. It is based upon the comparison of electrophoretic mobility in polyacrylamide gels of species containing two long arms, with the remaining one(s) being significantly shorter. Although the method currently lacks a really well-established basis of physical theory, it is very powerful, yet simple to apply. It has had a number of significant successes in RNA, DNA and DNA-protein complexes, and in all cases to date the results have stood the test of time and eventual comparison with crystallographic analysis.