The snRNP 15.5K protein folds its cognate K-turn RNA: a combined theoretical and biochemical study

SlRBP1 promotes translational efficiency via SleIF4A2 to maintain chloroplast function in tomato

Many glycine-rich RNA-binding proteins (GR-RBPs) have critical functions in RNA processing and metabolism. Here, we describe a role for the tomato (Solanum lycopersicum) GR-RBP SlRBP1 in regulating mRNA translation. We found that SlRBP1 knockdown mutants (slrbp1) displayed reduced accumulation of total chlorophyll and impaired chloroplast ultrastructure. These phenotypes were accompanied by deregulation of the levels of numerous key transcripts associated with chloroplast functions in slrbp1. Furthermore, native RNA immunoprecipitation-sequencing (nRIP-seq) recovered 61 SlRBP1-associated RNAs, most of which are involved in photosynthesis. SlRBP1 binding to selected target RNAs was validated by nRIP-qPCR. Intriguingly, the accumulation of proteins encoded by SlRBP1-bound transcripts, but not the mRNAs themselves, was reduced in slrbp1 mutants. Polysome profiling followed by RT-qPCR assays indicated that the polysome occupancy of target RNAs was lower in slrbp1 plants than in wild-type. Furthermore, SlRBP1 interacted with the eukaryotic translation initiation factor SleIF4A2. Silencing of SlRBP1 significantly reduced SleIF4A2 binding to SlRBP1-target RNAs. Taking these observations together, we propose that SlRBP1 binds to and channels RNAs onto the SleIF4A2 translation initiation complex and promotes the translation of its target RNAs to regulate chloroplast functions.

RNA kink-turns are highly anisotropic with respect to lateral displacement of the flanking stems

Kink-turns are highly bent internal loop motifs commonly found in the ribosome and other RNA complexes. They frequently act as binding sites for proteins and mediate tertiary interactions in larger RNA structures. Kink-turns have been a topic of intense research, but their elastic properties in the folded state are still poorly understood. Here we use extensive all-atom molecular dynamics simulations to parameterize a model of kink-turn in which the two flanking helical stems are represented by effective rigid bodies. Time series of the full set of six interhelical coordinates enable us to extract minimum energy shapes and harmonic stiffness constants for kink-turns from different RNA functional classes. The analysis suggests that kink-turns exhibit isotropic bending stiffness but are highly anisotropic with respect to lateral displacement of the stems. The most flexible lateral displacement mode is perpendicular to the plane of the static bend. These results may help understand the structural adaptation and mechanical signal transmission by kink-turns in complex natural and artificial RNA structures.

Classification and function of RNA-protein interactions

Almost all RNAs need to interact with proteins to fully exert their functions, and proteins also bind to RNAs to act as regulators. It has now become clear that RNA-protein interactions play important roles in many biological processes among organisms. Despite the great progress that has been made in the field, there is still no precise classification system for RNA-protein interactions, which makes it challenging to further decipher the functions and mechanisms of these interactions. In this review, we propose four different categories of RNA-protein interactions according to their basic characteristics: RNA motif-dependent RNA-protein interactions, RNA structure-dependent RNA-protein interactions, RNA modification-dependent RNA-protein interactions, and RNA guide-based RNA-protein interactions. Moreover, the integration of different types of RNA-protein interactions and the regulatory factors implicated in these interactions are discussed. Furthermore, we emphasize the functional diversity of these four types of interactions in biological processes and disease development and assess emerging trends in this exciting research field. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA Interactions with Proteins and Other Molecules > Protein-RNA Recognition RNA Processing > RNA Editing and Modification.

The yeast C/D box snoRNA U14 adopts a "weak" K-turn like conformation recognized by the Snu13 core protein in solution

Non-coding RNAs associate with proteins to form ribonucleoproteins (RNPs), such as ribosome, box C/D snoRNPs, H/ACA snoRNPs, ribonuclease P, telomerase and spliceosome to ensure cell viability. The assembly of these RNA-protein complexes relies on the ability of the RNA to adopt the correct bound conformation. K-turn motifs represent ubiquitous binding platform for proteins found in several cellular environment. This structural motif has an internal three-nucleotide bulge flanked on its 3' side by a G•A/A•G tandem pairs followed by one or two non-Watson-Crick pairs, and on its 5' side by a classical RNA helix. This peculiar arrangement induces a strong curvature of the phosphodiester backbone, which makes it conducive to multiple tertiary interactions. SNU13/Snu13p (Human/Yeast) binds specifically the U14 C/D box snoRNA K-turn sequence motif. This event is the prerequisite to promote the assembly of the RNP, which contains NOP58/Nop58 and NOP56/Nop56 core proteins and the 2'-O-methyl-transferase, Fibrillarin/Nop1p. The U14 small nucleolar RNA is a conserved non-coding RNA found in yeast and vertebrates required for the pre-rRNA maturation and ribose methylation. Here, we report the solution structure of the free U14 snoRNA K-turn motif (kt-U14) as determined by Nuclear Magnetic Resonance. We demonstrate that a major fraction of free kt-U14 adopts a pre-folded conformation similar to protein bound K-turn, even in the absence of divalent ions. In contrast to the kt-U4 or tyrS RNA, kt-U14 displays a sharp bent in the phosphodiester backbone. The U•U and G•A tandem base pairs are formed through weak hydrogen bonds. Finally, we show that the structure of kt-U14 is stabilized upon Snu13p binding. The structure of the free U14 RNA is the first reference example for the canonical motifs of the C/D box snoRNA family.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Cross-species conservation of complementary amino acid-ribonucleobase interactions and their potential for ribosome-free encoding

The role of amino acid-RNA nucleobase interactions in the evolution of RNA translation and protein-mRNA autoregulation remains an open area of research. We describe the inference of pairwise amino acid-RNA nucleobase interaction preferences using structural data from known RNA-protein complexes. We observed significant matching between an amino acid’s nucleobase affinity and corresponding codon content in both the standard genetic code and mitochondrial variants. Furthermore, we showed that knowledge of nucleobase preferences allows statistically significant prediction of protein primary sequence from mRNA using purely physiochemical information. Interestingly, ribosomal primary sequences were more accurately predicted than non-ribosomal sequences, suggesting a potential role for direct amino acid-nucleobase interactions in the genesis of amino acid-based ribosomal components. Finally, we observed matching between amino acid-nucleobase affinities and corresponding mRNA sequences in 35 evolutionarily diverse proteomes. We believe these results have important implications for the study of the evolutionary origins of the genetic code and protein-mRNA cross-regulation.

RNA structure determination by solid-state NMR spectroscopy

Knowledge of the RNA three-dimensional structure, either in isolation or as part of RNP complexes, is fundamental to understand the mechanism of numerous cellular processes. Because of its flexibility, RNA represents a challenge for crystallization, while the large size of cellular complexes brings solution-state NMR to its limits. Here, we demonstrate an alternative approach on the basis of solid-state NMR spectroscopy. We develop a suite of experiments and RNA labeling schemes and demonstrate for the first time that ssNMR can yield a RNA structure at high-resolution. This methodology allows structural analysis of segmentally labelled RNA stretches in high-molecular weight cellular machines—independent of their ability to crystallize— and opens the way to mechanistic studies of currently difficult-to-access RNA-protein assemblies.

The determination of RNA structures within high-molecular weight protein-RNA complexes in non-crystalline state is technically challenging. Here, the authors describe a solid-state NMR protocol for the determination of RNA structures at high resolution.

Molecular dynamic simulations of protein/RNA complexes: CRISPR/Csy4 endoribonuclease

BACKGROUND

Many prokaryotic genomes comprise Clustered Regularly Interspaced Short Palindromic Repeats (CRISPRs) offering defense against foreign nucleic acids. These immune systems are conditioned by the production of small CRISPR-derived RNAs matured from long RNA precursors. This often requires a Csy4 endoribonuclease cleaving the RNA 3'-end.

METHODS

We report extended explicit solvent molecular dynamic (MD) simulations of Csy4/RNA complex in precursor and product states, based on X-ray structures of product and inactivated precursor (55 simulations; ~3.7μs in total).

RESULTS

The simulations identify double-protonated His29 and deprotonated terminal phosphate as the likely dominant protonation states consistent with the product structure. We revealed potential substates consistent with Ser148 and His29 acting as the general base and acid, respectively. The Ser148 could be straightforwardly deprotonated through solvent and could without further structural rearrangements deprotonate the nucleophile, contrasting similar studies investigating the general base role of nucleobases in ribozymes. We could not locate geometries consistent with His29 acting as general base. However, we caution that the X-ray structures do not always capture the catalytically active geometries and then the reactive structures may be unreachable by the simulation technique.

CONCLUSIONS

We identified potential catalytic arrangement of the Csy4/RNA complex but we also report limitations of the simulation technique. Even for the dominant protonation state we could not achieve full agreement between the simulations and the structural data.

GENERAL SIGNIFICANCE

Potential catalytic arrangement of the Csy4/RNA complex is found. Further, we provide unique insights into limitations of simulations of protein/RNA complexes, namely, the influence of the starting experimental structures and force field limitations. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

Biochemical defects in minor spliceosome function in the developmental disorder MOPD I

Biallelic mutations of the human RNU4ATAC gene, which codes for the minor spliceosomal U4atac snRNA, cause the developmental disorder, MOPD I/TALS. To date, nine separate mutations in RNU4ATAC have been identified in MOPD I patients. Evidence suggests that all of these mutations lead to abrogation of U4atac snRNA function and impaired minor intron splicing. However, the molecular basis of these effects is unknown. Here, we use a variety of in vitro and in vivo assays to address this question. We find that only one mutation, 124G>A, leads to significantly reduced expression of U4atac snRNA, whereas four mutations, 30G>A, 50G>A, 50G>C and 51G>A, show impaired binding of essential protein components of the U4atac/U6atac di-snRNP in vitro and in vivo. Analysis of MOPD I patient fibroblasts and iPS cells homozygous for the most common mutation, 51G>A, shows reduced levels of the U4atac/U6atac.U5 tri-snRNP complex as determined by glycerol gradient sedimentation and immunoprecipitation. In this report, we establish a mechanistic basis for MOPD I disease and show that the inefficient splicing of genes containing U12-dependent introns in patient cells is due to defects in minor tri-snRNP formation, and the MOPD I-associated RNU4ATAC mutations can affect multiple facets of minor snRNA function.

Interactions of noncanonical motifs with hnRNP A2 promote activity-dependent RNA transport in neurons

Ca²⁺-dependent RNA–protein interactions enable activity-inducible RNA transport in dendrites.

A key determinant of neuronal functionality and plasticity is the targeted delivery of select ribonucleic acids (RNAs) to synaptodendritic sites of protein synthesis. In this paper, we ask how dendritic RNA transport can be regulated in a manner that is informed by the cell’s activity status. We describe a molecular mechanism in which inducible interactions of noncanonical RNA motif structures with targeting factor heterogeneous nuclear ribonucleoprotein (hnRNP) A2 form the basis for activity-dependent dendritic RNA targeting. High-affinity interactions between hnRNP A2 and conditional GA-type RNA targeting motifs are critically dependent on elevated Ca²⁺ levels in a narrow concentration range. Dendritic transport of messenger RNAs that carry such GA motifs is inducible by influx of Ca²⁺ through voltage-dependent calcium channels upon β-adrenergic receptor activation. The combined data establish a functional correspondence between Ca²⁺-dependent RNA–protein interactions and activity-inducible RNA transport in dendrites. They also indicate a role of genomic retroposition in the phylogenetic development of RNA targeting competence.

Kink turn sRNA folding upon L7Ae binding using molecular dynamics simulations

The kink-turn sRNA motif in archaea, whose combination with protein L7Ae initializes the assembly of small ribonucleoprotein particles (sRNPs), plays a key role in ribosome maturation and the translation process. Although many studies have been reported on this motif, the mechanism of sRNA folding coupled with protein binding is still poorly understood. Here, room and high temperature molecular dynamics (MD) simulations were performed on the complex of 25-nt kink-turn sRNA and L7Ae. The average RMSD values between the bound and corresponding apo structures and Kolmogorov-Smirnov P test analysis indicate that sRNA may follow an induced fit mechanism upon binding with L7Ae, both locally and globally. These conclusions are further supported by high-temperature unfolding kinetic analysis. Principal component analysis (PCA) found both closing and opening motions of the kink-turn sRNA. This might play a key role in the sRNP assembly and methylation catalysis. These combined computational methods can be used to study the specific recognition of other sRNAs and proteins.

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.

Structural basis for the dual U4 and U4atac snRNA-binding specificity of spliceosomal protein hPrp31

Human proteins 15.5K and hPrp31 are components of the major spliceosomal U4 snRNP and of the minor spliceosomal U4atac snRNP. The two proteins bind to related 5'-stem loops (5'SLs) of the U4 and U4atac snRNAs in a strictly sequential fashion. The primary binding 15.5K protein binds at K-turns that exhibit identical sequences in the two snRNAs. However, RNA sequences contacted by the secondary binding hPrp31 differ in U4 and U4atac snRNAs, and the mechanism by which hPrp31 achieves its dual specificity is presently unknown. We show by crystal structure analysis that the capping pentaloops of the U4 and U4atac 5'SLs adopt different structures in the ternary hPrp31-15.5K-snRNA complexes. In U4atac snRNA, a noncanonical base pair forms across the pentaloop, based on which the RNA establishes more intimate interactions with hPrp31 compared with U4 snRNA. Stacking of hPrp31-His270 on the noncanonical base pair at the base of the U4atac pentaloop recapitulates intramolecular stabilizing principles known from the UUCG and GNRA families of RNA tetraloops. Rational mutagenesis corroborated the importance of the noncanonical base pair and the U4atac-specific hPrp31-RNA interactions for complex stability. The more extensive hPrp31-U4atac snRNA interactions are in line with a higher stability of the U4atac compared with the U4-based ternary complex seen in gel-shift assays, which may explain how U4atac snRNA can compete with the more abundant U4 snRNA for the same protein partners in vivo.

Spatial code recognition in neuronal RNA targeting: role of RNA-hnRNP A2 interactions

Recognition of non-canonical purine•purine RNA motifs by hnRNP A2 mediates targeted delivery of neuronal RNAs to dendrites.

In neurons, regulation of gene expression occurs in part through translational control at the synapse. A fundamental requirement for such local control is the targeted delivery of select neuronal mRNAs and regulatory RNAs to distal dendritic sites. The nature of spatial RNA destination codes, and the mechanism by which they are interpreted for dendritic delivery, remain poorly understood. We find here that in a key dendritic RNA transport pathway (exemplified by BC1 RNA, a dendritic regulatory RNA, and protein kinase M ζ [PKMζ] mRNA, a dendritic mRNA), noncanonical purine•purine nucleotide interactions are functional determinants of RNA targeting motifs. These motifs are specifically recognized by heterogeneous nuclear ribonucleoprotein A2 (hnRNP A2), a trans-acting factor required for dendritic delivery. Binding to hnRNP A2 and ensuing dendritic delivery are effectively competed by RNAs with CGG triplet repeat expansions. CGG repeats, when expanded in the 5′ untranslated region of fragile X mental retardation 1 (FMR1) mRNA, cause fragile X–associated tremor/ataxia syndrome. The data suggest that cellular dysregulation observed in the presence of CGG repeat RNA may result from molecular competition in neuronal RNA transport pathways.

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.

Structure of the K-turn U4 RNA: a combined NMR and SANS study

K-turn motifs are universal RNA structural elements providing a binding platform for proteins in several cellular contexts. Their characteristic is a sharp kink in the phosphate backbone that puts the two helical stems of the protein-bound RNA at an angle of 60°. However, to date no high-resolution structure of a naked K-turn motif is available. Here, we present the first structural investigation at atomic resolution of an unbound K-turn RNA (the spliceosomal U4-Kt RNA) by a combination of NMR and small-angle neutron scattering data. With this study, we wish to address the question whether the K-turn structural motif assumes the sharply kinked conformation in the absence of protein binders and divalent cations. Previous studies have addressed this question by fluorescence resonance energy transfer, biochemical assays and molecular dynamics simulations, suggesting that the K-turn RNAs exist in equilibrium between a kinked conformation, which is competent for protein binding, and a more extended conformation, with the population distribution depending on the concentration of divalent cations. Our data shows that the U4-Kt RNA predominantly assumes the more extended conformation in the absence of proteins and divalent cations. The internal loop region is well structured but adopts a different conformation from the one observed in complex with proteins. Our data suggests that the K-turn consensus sequence does not per se code for the kinked conformation; instead the sharp backbone kink requires to be stabilized by protein binders.

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.

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.

Dynamic guide-target interactions contribute to sequential 2'-O-methylation by a unique archaeal dual guide box C/D sRNP

Assembly and guide-target interaction of an archaeal box C/D-guide sRNP was investigated under various conditions by analyzing the lead (II)-induced cleavage of the guide RNA. Guide and target RNAs derived from Haloferax volcanii pre-tRNA(Trp) were used with recombinant Methanocaldococcus jannaschii core proteins in the reactions. Core protein L7Ae binds differentially to C/D and C'/D' motifs of the guide RNA, and interchanging the two motifs relative to the termini of the guide RNA did not affect L7Ae binding or sRNA function. L7Ae binding to the guide RNA exposes its D'-guide sequence first followed by the D guide. These exposures are reduced when aNop5p and aFib proteins are added. The exposed guide sequences did not pair with the target sequences in the presence of L7Ae alone. The D-guide sequence could pair with the target in the presence of L7Ae and aNop5p, suggesting a role of aNop5p in target recruitment and rearrangement of sRNA structure. aFib binding further stabilizes this pairing. After box C/D-guided modification, target-guide pairing at the D-guide sequence is disrupted, suggesting that each round of methylation may require some conformational change or reassembly of the RNP. Asymmetric RNPs containing only one L7Ae at either of the two box motifs can be assembled, but a functional RNP requires L7Ae at the box C/D motif. This arrangement resembles the asymmetric eukaryal snoRNP. Observations of initial D-guide-target pairing and the functional requirement for L7Ae at the box C/D motif are consistent with our previous report of the sequential 2'-O-methylations of the target RNA.

Conformations of flanking bases in HIV-1 RNA DIS kissing complexes studied by molecular dynamics

Explicit solvent molecular dynamics simulations (in total almost 800 ns including locally enhanced sampling runs) were applied with different ion conditions and with two force fields (AMBER and CHARMM) to characterize typical geometries adopted by the flanking bases in the RNA kissing-loop complexes. We focus on flanking base positions in multiple x-ray and NMR structures of HIV-1 DIS kissing complexes and kissing complex from the large ribosomal subunit of Haloarcula marismortui. An initial x-ray open conformation of bulged-out bases in HIV-1 DIS complexes, affected by crystal packing, tends to convert to a closed conformation formed by consecutive stretch of four stacked purine bases. This is in agreement with those recent crystals where the packing is essentially avoided. We also observed variants of the closed conformation with three stacked bases, while nonnegligible populations of stacked geometries with bulged-in bases were detected, too. The simulation results reconcile differences in positions of the flanking bases observed in x-ray and NMR studies. Our results suggest that bulged-out geometries are somewhat more preferred, which is in accord with recent experiments showing that they may mediate tertiary contacts in biomolecular assemblies or allow binding of aminoglycoside antibiotics.

Elastic properties of ribosomal RNA building blocks: molecular dynamics of the GTPase-associated center rRNA

Explicit solvent molecular dynamics (MD) was used to describe the intrinsic flexibility of the helix 42–44 portion of the 23S rRNA (abbreviated as Kt-42+rGAC; kink-turn 42 and GTPase-associated center rRNA). The bottom part of this molecule consists of alternating rigid and flexible segments. The first flexible segment (Hinge1) is the highly anharmonic kink of Kt-42. The second one (Hinge2) is localized at the junction between helix 42 and helices 43/44. The rigid segments are the two arms of helix 42 flanking the kink. The whole molecule ends up with compact helices 43/44 (Head) which appear to be modestly compressed towards the subunit in the Haloarcula marismortui X-ray structure. Overall, the helix 42–44 rRNA is constructed as a sophisticated intrinsically flexible anisotropic molecular limb. The leading flexibility modes include bending at the hinges and twisting. The Head shows visible internal conformational plasticity, stemming from an intricate set of base pairing patterns including dynamical triads and tetrads. In summary, we demonstrate how rRNA building blocks with contrasting intrinsic flexibilities can form larger architectures with highly specific patterns of preferred low-energy motions and geometries.

Analysis of sequence and structural features that identify the B/C motif of U3 small nucleolar RNA as the recognition site for the Snu13p-Rrp9p protein pair

The eukaryal Snu13p/15.5K protein binds K-turn motifs in U4 snRNA and snoRNAs. Two Snu13p/15.5K molecules bind the nucleolar U3 snoRNA required for the early steps of preribosomal processing. Binding of one molecule on the C'/D motif allows association of proteins Nop1p, Nop56p, and Nop58p, whereas binding of the second molecule on the B/C motif allows Rrp9p recruitment. To understand how the Snu13p-Rrp9p pair recognizes the B/C motif, we first improved the identification of RNA determinants required for Snu13p binding by experiments using the systematic evolution of ligands by exponential enrichment. This demonstrated the importance of a U.U pair stacked on the sheared pairs and revealed a direct link between Snu13p affinity and the stability of helices I and II. Sequence and structure requirements for efficient association of Rrp9p on the B/C motif were studied in yeast cells by expression of variant U3 snoRNAs and immunoselection assays. A G-C pair in stem II, a G residue at position 1 in the bulge, and a short stem I were found to be required. The data identify the in vivo function of most of the conserved residues of the U3 snoRNA B/C motif. They bring important information to understand how different K-turn motifs can recruit different sets of proteins after Snu13p association.

Structure, dynamics, and elasticity of free 16s rRNA helix 44 studied by molecular dynamics simulations

Molecular dynamics (MD) simulations were employed to investigate the structure, dynamics, and local base-pair step deformability of the free 16S ribosomal helix 44 from Thermus thermophilus and of a canonical A-RNA double helix. While helix 44 is bent in the crystal structure of the small ribosomal subunit, the simulated helix 44 is intrinsically straight. It shows, however, substantial instantaneous bends that are isotropic. The spontaneous motions seen in simulations achieve large degrees of bending seen in the X-ray structure and would be entirely sufficient to allow the dynamics of the upper part of helix 44 evidenced by cryo-electron microscopic studies. Analysis of local base-pair step deformability reveals a patch of flexible steps in the upper part of helix 44 and in the area proximal to the bulge bases, suggesting that the upper part of helix 44 has enhanced flexibility. The simulations identify two conformational substates of the second bulge area (bottom part of the helix) with distinct base pairing. In agreement with nuclear magnetic resonance (NMR) and X-ray studies, a flipped out conformational substate of conserved 1492A is seen in the first bulge area. Molecular dynamics (MD) simulations reveal a number of reversible alpha-gamma backbone flips that correspond to transitions between two known A-RNA backbone families. The flipped substates do not cumulate along the trajectory and lead to a modest transient reduction of helical twist with no significant influence on the overall geometry of the duplexes. Despite their considerable flexibility, the simulated structures are very stable with no indication of substantial force field inaccuracies.

RNA kink-turns as molecular elbows: hydration, cation binding, and large-scale dynamics

The presence of Kink-turns (Kt) at key functional sites in the ribosome (e.g., A-site finger and L7/L12 stalk) suggests that some Kink-turns can confer flexibility on RNA protuberances that regulate the traversal of tRNAs during translocation. Explicit solvent molecular dynamics demonstrates that Kink-turns can act as flexible molecular elbows. Kink-turns are associated with a unique network of long-residency static and dynamical hydration sites that is intimately involved in modulating their conformational dynamics. An implicit solvent conformational search confirms the flexibility of Kink-turns around their X-ray geometries and identifies a second low-energy region with open structures that could correspond to Kink-turn geometries seen in solution experiments. An extended simulation of Kt-42 with the factor binding site (helices 43 and 44) shows that the local Kt-42 elbow-like motion fully propagates beyond the Kink-turn, and that there is no other comparably flexible site in this rRNA region. Kink-turns could mediate large-scale adjustments of distant RNA segments.

Consecutive GA pairs stabilize medium-size RNA internal loops

Internal loops in RNA are important for folding and function. Consecutive noncanonical pairs can form in internal loops having at least two nucleotides on each side. Thermodynamic and structural insights into such internal loops should improve approximations for their stabilities and predictions of secondary and three-dimensional structures. Most natural internal loops are purine rich. A series of oligoribonucleotides that form purine-rich internal loops of 5-10 nucleotides, including kink-turn loops, were studied by UV melting, exchangeable proton and phosphorus NMR. Three consecutive GA pairs with the motif 5' Y GGA/3' R AAG or GGA R 3'/AAG Y 5' (i.e., 5' GGA 3'/3' AAG 5' closed on at least one side with a CG, UA, or UG pair with Y representing C or U and R representing A or G) stabilize internal loops having 6-10 nucleotides. Certain motifs with two consecutive GA pairs are also stabilizing. In internal loops with three or more nucleotides on each side, the motif 5' U G/3' G A has stability similar to 5' C G/3' G A. A revised model for predicting stabilities of internal loops with 6-10 nucleotides is derived by multiple linear regression. Loops with 2 x 3 nucleotides are predicted well by a previous thermodynamic model.

Molecular dynamics simulations of sarcin-ricin rRNA motif

Explicit solvent molecular dynamics (MD) simulations were carried out for sarcin-ricin domain (SRD) motifs from 23S (Escherichia coli) and 28S (rat) rRNAs. The SRD motif consists of GAGA tetraloop, G-bulged cross-strand A-stack, flexible region and duplex part. Detailed analysis of the overall dynamics, base pairing, hydration, cation binding and other SRD features is presented. The SRD is surprisingly static in multiple 25 ns long simulations and lacks any non-local motions, with root mean square deviation (r.m.s.d.) values between averaged MD and high-resolution X-ray structures of 1-1.4 A. Modest dynamics is observed in the tetraloop, namely, rotation of adenine in its apex and subtle reversible shift of the tetraloop with respect to the adjacent base pair. The deformed flexible region in low-resolution rat X-ray structure is repaired by simulations. The simulations reveal few backbone flips, which do not affect positions of bases and do not indicate a force field imbalance. Non-Watson-Crick base pairs are rigid and mediated by long-residency water molecules while there are several modest cation-binding sites around SRD. In summary, SRD is an unusually stiff rRNA building block. Its intrinsic structural and dynamical signatures seen in simulations are strikingly distinct from other rRNA motifs such as Loop E and Kink-turns.

Small non-coding RNAs in Archaea

Biochemical and informatics analyses conducted over the past few years have revealed the presence of a plethora of small non-coding RNAs in various species of Archaea. A large proportion of these RNAs contain a common structural motif called the RNA kink turn (K-turn). The best-characterized are the C/D box and the H/ACA box guide small (s)RNAs. Both contain the K-turn fold and require the binding of the L7Ae protein to stabilize the structure of this crucial motif. These sRNAs assemble with L7Ae and several other proteins into complex and dynamic ribonucleoprotein machines that mediate guide-directed ribose methylation or pseudouridylation to specific locations in ribosomal or transfer RNA. Analyses of new archaeal sRNA libraries have identified additional classes of novel sRNAs; many of these contain the RNA K-turn motif and suggest that the RNAs might function as ribonucleoprotein complexes. Some have characteristics of small interfering RNAs or of micro RNAs that have been implicated in the post-transcriptional control of gene expression, whereas others appear to be involved in protein translocation or in ribosomal RNA processing and ribosome assembly. A complete understanding of the structure of the K-turn motif and its contribution to various RNA-RNA and RNA-protein interactions will be absolutely essential to fully elucidate the biological organization, activity and function of these novel archaeal ribonucleoprotein machines.

Conformational dynamics of RNA-peptide binding: a molecular dynamics simulation study

Molecular dynamics simulations of the binding of the heterochiral tripeptide KkN to the transactivation responsive (TAR) RNA of HIV-1 is presented, using an all-atom force field with explicit water. To obtain starting structures for the TAR-KkN complex, semirigid docking calculations were performed that employ an NMR structure of free TAR RNA. The molecular dynamics simulations show that the starting structures in which KkN binds to the major groove of TAR (as it is the case for the Tat-TAR complex of HIV-1) are unstable. On the other hand, the minor-groove starting structures are found to lead to several binding modes, which are stabilized by a complex interplay of stacking, hydrogen bonding, and electrostatic interactions. Although the ligand does not occupy the binding position of Tat protein, it is shown to hinder the interhelical motion of free TAR RNA. The latter is presumably necessary to achieve the conformational change of TAR RNA to bind Tat protein. Considering the time evolution of the trajectories, the binding process is found to be ligand-induced and cooperative. That is, the conformational rearrangement only occurs in the presence of the ligand and the concerted motion of the ligand and a large part of the RNA binding site is necessary to achieve the final low-energy binding state.

Detecting protein-induced folding of the U4 snRNA kink-turn by single-molecule multiparameter FRET measurements

The kink-turn (k-turn), a new RNA structural motif found in the spliceosome and the ribosome, serves as a specific protein recognition element and as a structural building block. While the structure of the spliceosomal U4 snRNA k-turn/15.5K complex is known from a crystal structure, it is unclear whether the k-turn also exists in this folded conformation in the free U4 snRNA. Thus, we investigated the U4 snRNA k-turn by single-molecule FRET measurements in the absence and presence of the 15.5K protein and its dependence on the Na(+) and Mg(2+) ion concentration. We show that the unfolded U4 snRNA k-turn introduces a kink of 85 degrees +/- 15 degrees in an RNA double helix. While Na(+) and Mg(2+) ions induce this more open conformation of the k-turn, binding of the 15.5K protein was found to induce the tightly kinked conformation in the RNA that increases the kink to 52 degrees +/- 15 degrees . By comparison of the measured FRET distances with a computer-modeled structure, we show that this strong kink is due to the k-turn motif adopting its folded conformation. Thus, in the free U4 snRNA, the k-turn exists only in an unfolded conformation, and its folding is induced by binding of the 15.5K protein.

Molecular basis for RNA kink-turn recognition by the h15.5K small RNP protein

The interaction between box C/D small nucleolar (sno)RNAs and the 15.5K protein nucleates snoRNP assembly. Many eukaryotic snoRNAs contain two potential binding sites for this protein, only one of which appears to be utilized in vivo. The binding site conforms to the consensus for a kink-turn motif. We have investigated the molecular basis for selection of one potential site over the other using in vitro mobility shift assays and nucleotide analog interference mapping of Xenopus U25 snoRNA and of a circularly permuted form. We find that preferential binding of human 15.5K is not dependent on the proximity of RNA ends, but instead appears to require a structural context beyond the kink-turn itself. Direct analysis of the energetic contributions to binding made by 18 functional groups within the kink-turn identified both backbone atoms and base functionalities as key for interaction. An intramolecular RNA-RNA contact via a 2'-hydroxyl may supercede a putative Type I A-minor interaction in stabilizing the RNA-protein complex.

Induced fit of RNA on binding the L7Ae protein to the kink-turn motif

The kink-turn is a widespread motif in RNA consisting of a three-nucleotide bulge flanked on one side by consecutive A3G mismatches. Important examples are found in the ribosome, U4 RNA, and in snoRNAs involved in RNA modification. The motif is a common protein binding site, and the RNA has been found to adopt a tightly kinked conformation in crystal structures. However, in free solution there is a dynamic exchange between kinked and extended conformations, with the equilibrium driven toward the kinked form by the addition of metal ions. Here we used fluorescence resonance energy transfer (FRET) to show that the L7Ae protein of Archaeoglobus fulgidus binds to RNA containing a kink-turn with nanomolar affinity, and induces folding into the tightly kinked conformation even in the absence of metal ions. Thus this RNA may act as a relatively flexible hinge during RNA folding, until fixed into its ultimate kinked structure by the binding of L7 or related protein.

Loss of G-A base pairs is insufficient for achieving a large opening of U4 snRNA K-turn motif

Upon binding to the 15.5K protein, two tandem-sheared G–A base pairs are formed in the internal loop of the kink-turn motif of U4 snRNA (Kt-U4). We have reported that the folding of Kt-U4 is assisted by protein binding. Unstable interactions that contribute to a large opening of the free RNA (‘k–e motion’) were identified using locally enhanced sampling molecular dynamics simulations, results that agree with experiments. A detailed analysis of the simulations reveals that the k–e motion in Kt-U4 is triggered both by loss of G–A base pairs in the internal loop and backbone flexibility in the stems. Essential dynamics show that the loss of G–A base pairs is correlated along the first mode but anti-correlated along the third mode with the k–e motion. Moreover, when enhanced sampling was confined to the internal loop, the RNA adopted an alternative conformation characterized by a sharper kink, opening of G–A base pairs and modified stacking interactions. Thus, loss of G–A base pairs is insufficient for achieving a large opening of the free RNA. These findings, supported by previously published RNA structure probing experiments, suggest that G–A base pair formation occurs upon protein binding, thereby stabilizing a selective orientation of the stems.

Recurrent structural RNA motifs, Isostericity Matrices and sequence alignments

The occurrences of two recurrent motifs in ribosomal RNA sequences, the Kink-turn and the C-loop, are examined in crystal structures and systematically compared with sequence alignments of rRNAs from the three kingdoms of life in order to identify the range of the structural and sequence variations. Isostericity Matrices are used to analyze structurally the sequence variations of the characteristic non-Watson-Crick base pairs for each motif. We show that Isostericity Matrices for non-Watson-Crick base pairs provide important tools for deriving the sequence signatures of recurrent motifs, for scoring and refining sequence alignments, and for determining whether motifs are conserved throughout evolution. The systematic use of Isostericity Matrices identifies the positions of the insertion or deletion of one or more nucleotides relative to the structurally characterized examples of motifs and, most importantly, specifies whether these changes result in new motifs. Thus, comparative analysis coupled with Isostericity Matrices allows one to produce and refine structural sequence alignments. The analysis, based on both sequence and structure, permits therefore the evaluation of the conservation of motifs across phylogeny and the derivation of rules of equivalence between structural motifs. The conservations observed in Isostericity Matrices form a predictive basis for identifying motifs in sequences.