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 (∼1 kcal mol⁻¹), separated by a small free-energy barrier (∼1.5 kcal mol⁻¹) 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.     

Ribosomal RNA kink-turn motif--a flexible molecular hinge

Ribosomal RNA K-turn motifs are asymmetric internal loops characterized by a sharp bend in the phosphodiester backbone resulting in "V" shaped structures, recurrently observed in ribosomes and showing a high degree of sequence conservation. We have carried out extended explicit solvent molecular dynamics simulations of selected K-turns, in order to investigate their intrinsic structural and dynamical properties. The simulations reveal an unprecedented dynamical flexibility of the K-turns around their X-ray geometries. The K-turns sample, on the nanosecond timescale, different conformational substates. The overall behavior of the simulations suggests that the sampled geometries are essentially isoenergetic and separated by minimal energy barriers. The nanosecond dynamics of isolated K-turns can be qualitatively considered as motion of two rigid helix stems controlled by a very flexible internal loop which then leads to substantial hinge-like motions between the two stems. This internal dynamics of K-turns is strikingly different for example from the bacterial 5S rRNA Loop E motif or BWYV frameshifting pseudoknot which appear to be rigid in the same type of simulations. Bistability and flexibility of K-turns was also suggested by several recent biochemical studies. Although the results of MD simulations should be considered as a qualitative picture of the K-turn dynamics due to force field and sampling limitations, the main advantage of the MD technique is its ability to investigate the region close to K-turn ribosomal-like geometries. This part of the conformational space is not well characterized by the solution experiments due to large-scale conformational changes seen in the experiments. We suggest that K-turns are well suited to act as flexible structural elements of ribosomal RNA. They can for example be involved in mediation of large-scale motions or they can allow a smooth assembling of the other parts of the ribosome.     

Ribosomal RNA Kink-turn Motif—A Flexible Molecular Hinge

Abstract

Ribosomal RNA K-turn motifs are asymmetric internal loops characterized by a sharp bend in the phosphodiester backbone resulting in “V” shaped structures, recurrently observed in ribosomes and showing a high degree of sequence conservation. We have carried out extended explicit solvent molecular dynamics simulations of selected K-turns, in order to investigate their intrinsic structural and dynamical properties. The simulations reveal an unprecedented dynamical flexibility of the K-turns around their X-ray geometries. The K-turns sample, on the nanosecond timescale, different conformational substates. The overall behavior of the simulations suggests that the sampled geometries are essentially isoenergetic and separated by minimal energy barriers. The nanosecond dynamics of isolated K-turns can be qualitatively considered as motion of two rigid helix stems controlled by a very flexible internal loop which then leads to substantial hinge-like motions between the two stems. This internal dynamics of K-turns is strikingly different for example from the bacterial 5S rRNA Loop E motif or BWYV frameshifting pseudoknot which appear to be rigid in the same type of simulations. Bistability and flexibility of K-turns was also suggested by several recent biochemical studies. Although the results of MD simulations should be considered as a qualitative picture of the K-turn dynamics due to force field and sampling limitations, the main advantage of the MD technique is its ability to investigate the region close to K-turn riboso- mal-like geometries. This part of the conformational space is not well characterized by the solution experiments due to large-scale conformational changes seen in the experiments. We suggest that K-turns are well suited to act as flexible structural elements of ribosomal RNA. They can for example be involved in mediation of large-scale motions or they can allow a smooth assembling of the other parts of the ribosome.     

The kink-turn: a new RNA secondary structure motif

Analysis of the Haloarcula marismortui large ribosomal subunit has revealed a common RNA structure that we call the kink-turn, or K-turn. The six K-turns in H.marismortui 23S rRNA superimpose with an r.m.s.d. of 1.7 A. There are two K-turns in the structure of Thermus thermophilus 16S rRNA, and the structures of U4 snRNA and L30e mRNA fragments form K-turns. The structure has a kink in the phosphodiester backbone that causes a sharp turn in the RNA helix. Its asymmetric internal loop is flanked by C-G base pairs on one side and sheared G-A base pairs on the other, with an A-minor interaction between these two helical stems. A derived consensus secondary structure for the K-turn includes 10 consensus nucleotides out of 15, and predicts its presence in the 5'-UTR of L10 mRNA, helix 78 in Escherichia coli 23S rRNA and human RNase MRP. Five K-turns in 23S rRNA interact with nine proteins. While the observed K-turns interact with proteins of unrelated structures in different ways, they interact with L7Ae and two homologous proteins in the same way.     

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 Å, reaches the small subunit through intersubunit bridge B1a. Previous cryo-EM studies revealed that the tip of H38 moves by more than 10 Å 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.     

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–44rRNA 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.     

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.     

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.     

The K-loop, a general feature of the Pyrococcus C/D guide RNAs, is an RNA structural motif related to the K-turn

The C/D guide RNAs predicted from the genomic sequences of three species of Pyrococcus delineate a family of small non-coding archaeal RNAs involved in the methylation of rRNA and tRNA. The C/D guides assemble into ribonucleoprotein (RNP) that contains the methyltransferase. The protein L7Ae, a key structural component of the RNP, binds to a Kink-turn (K-turn) formed by the C/D motif. The K-turn is a structure that consists of two RNA stems separated by a short asymmetric loop with a characteristic sharp bend (kink) between the two stems. The majority of the pyrococcal C/D guides contain a short 3 nt-spacer between the C′/D′ motifs. We show here that conserved terminal stem–loops formed by the C′/D′ motif of the Pyrococcus C/D RNAs are also L7Ae-binding sites. These stem–loops are related to the K-turn by sequence and structure, but they consist of a single stem closed by a terminal loop. We have named this structure the K-loop. We show that conserved non-canonical base pairs in the stem of the K-loop are necessary for L7Ae binding. For the C/D guides with a 3 nt-spacer we show that the sequence and length is also important. The K-loop could improve the stability of the C/D guide RNAs in Pyrococcal species, which are extreme hyperthermophiles.     

The importance of G.A hydrogen bonding in the metal ion- and protein-induced folding of a kink turn RNA

The kink turn (K-turn) is a common motif in RNA structure, found in many RNA species important in translation, RNA modification and splicing, and the control of gene expression. In general the K-turn comprises a three nucleotide bulge followed by trans sugar-Hoogsteen G·A pairs. The RNA adopts a tightly kinked conformation, and is a common target for binding proteins, exemplified by the L7Ae family. We have measured the rates of association and dissociation for the binding of L7Ae to the Kt-7 kink turn, from which we calculate an affinity of K_D = 10 pM. This high affinity is consistent with the role of this binding as the first stage in the assembly of key functional nucleoproteins such as box C/D snoRNP. Kink-turn RNA undergoes a two-state transition between the kinked conformation, and a more extended structure, and folding into the kinked form is induced by divalent metal ions, or by binding of proteins of the L7Ae class. The K-turn provides an excellent, simple model for RNA folding, which can be dissected at the atomic level. We have analyzed the contributions of the hydrogen bonds that form the G·A pairs to the ion- and protein-induced folding of the K-turn. We find that all four hydrogen bonds are important to the stability of the kinked form of the RNA, and we can now define all the important hydrogen bonding interactions that stabilize the K-turn. The high affinity of L7Ae binding is coupled to the induced folding of the K-turn, allowing some sub-optimal variants to adopt the kinked geometry. However, in all such cases the affinity is lowered, and the results underline the importance of both G·A pairs to the stability of the K-turn.     

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.     

Structure of protein L7Ae bound to a K-turn derived from an archaeal box H/ACA sRNA at 1.8 A resolution

The archaeal RNA binding protein L7Ae and its eukaryotic homolog 15.5 kDa/Snu13 recognize K-turns. This structural motif is canonically comprised of two stems (one with tandem A.G base pairs, the other with Watson-Crick pairs) linked by an asymmetric internal loop. L7Ae recognizes conventional K-turns in ribosomal and box C/D RNAs but also binds specifically to some box H/ACA RNAs at terminal stem loops. These have the A.G paired stem, but lack the Watson-Crick stem. The structure of Methanococcus jannaschii L7Ae bound to a symmetric duplex RNA without Watson-Crick stems demonstrates how a binding site for this component of diverse ribonucleoprotein complexes can be constructed with only the A.G stem and the loop. The RNA adopts a functional conformation with the aid of a base triple and tight binding of divalent cations. Comparison with the 15.5 kDa/Snu13-RNA complex structure suggests why the eukaryotic homolog does not recognize terminal stem loop L7Ae binding sites.     

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

The role of 2'-hydroxyl groups in stabilizing the tightly kinked geometry of the kink-turn (K-turn) has been investigated. Individual 2'-OH groups have been removed by chemical synthesis, and the kinking of the RNA has been studied by gel electrophoresis and fluorescence resonance energy transfer. The results have been analyzed by reference to a database of 11 different crystallographic structures of K-turns. The potential hydrogen bonds fall into several classes. The most important are those in the core of the turn and ribose-phosphate interactions around the bulge. Of these the single most important hydrogen bond is one donated from the 2'-OH of the 5' nucleotide of the bulge to the N1 of the adenine of the kink-proximal A*G pair. This is present in all known K-turn structures, and removal of the 2'-OH completely prevents metal ion-induced folding. Hydrogen bonds formed in the minor grooves of the helical stems are less important, and removal of the participating 2'-OH groups leads to reduced impairment of folding. These interactions are generally more polymorphic, and hydrogen bonds probably form where possible, as permitted by the global structure.     

A structural,phylogenetic, and functional study of 15.5-kD/Snu13 protein binding on U3 small nucleolar RNA

The 15.5-kD protein and its yeast homolog Snu13p bind U4 snRNA, U3 snoRNA, and the C/D box snoRNAs. In U4 snRNA, they associate with a helix-bulge-helix (K-turn) structure. U3 snoRNA contains two conserved pairs of boxes, C'/D and B/C, which were both expected to bind the 15.5-kD/Snu13 protein. Only binding to the B/C motif was experimentally demonstrated. Here, by chemical probing of in vitro reconstituted RNA/protein complexes, we demonstrate the independent binding of the 15.5-kD/Snu13 protein to each of the two motifs. Due to a highly reduced stem I (1 bp), the K-turn structure is not formed in the naked B/C motif. However, gel-shift experiments revealed a higher affinity of Snu13p for the B/C motif, compared to the C'/D motif. A phylogenetic analysis of U3 snoRNA, coupled with an analysis of Snu13p affinity for variant yeast C'/D and B/C motifs, and a study of the functionality of a truncated yeast U3 snoRNA carrying base substitutions in the C'/D and B/C motifs, revealed that conservation of the identities of residues 2 and 3 in the B/C K-turn is more important for Snu13p binding and U3 snoRNA function, than conservation of the identities of corresponding residues in the C'/D K-turn. This suggests that binding of Snu13p to K-turns with a very short helix I imposes sequence constraints in the bulge. Altogether, the data demonstrate the strong importance of the binding of the 15.5-kD/Snu13 protein to the C'/D and B/C motifs for both U3 snoRNP assembly and activity.     

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.     

Functional implications for a prototypical K-turn binding protein from structural and dynamical studies of 15.5K

The kink-turn (K-turn) motif is recognized and bound by a family of proteins that act as nucleation factors for ribonucleoparticle assembly. The binding of various proteins to a conserved RNA structural motif known as the K-turn has been shown to be an important component of regulation in the ribosome, in the spliceosome, and in RNA modification. 15.5K is a prototypical example of a K-turn binding protein, which has been shown to bind the 5'-U4 stem-loop of the spliceosome and the box C/D motif. We describe the solution NMR structure of free 15.5K, as well as studies of conformational flexibility from 15N NMR relaxation and H/D exchange experiments. The protein appears well-structured aside from conformational fluctuation in alpha3. Flexibility in fast time scale motions and the observation of limited intermediate and slow motions further characterize the free protein and may suggest local contributions to recognition and binding.     

A molecular dynamics study of the C-terminal fragment of the L7/L12 ribosomal protein. Secondary structure motion in a 150 picosecond trajectory

A 150 picosecond molecular dynamics computer simulation of the C-terminal fragment of the L7/L12 ribosomal protein from Escherichia coli is reported. The molecular dynamics results are compared with the available high-resolution X-ray data in terms of atomic positions, distances and positional fluctuations. Good agreement is found between the molecular dynamics results and the X-ray data. The form and parameters of the interaction potential energy function and the procedures for deriving it are discussed. Some current misunderstandings concerning the ways of evaluating the efficiency of molecular dynamics algorithms and of application of bond-length constraints in protein simulations are cleared up. The 150 picosecond trajectory has been scanned in a search for correlated motions within and between secondary structure elements. The beta-strands have diffusional stretching modes, and uncorrelated transversal displacements. The dynamic analysis of alpha-helices shows a variety of features. The atomic fluctuations differ between the helix ends; this effect reflects long time-scale motions. Two alpha-helices, alpha A and alpha C, show diffusive longitudinal stretching modes. The third helix, alpha B, has a correlated asymmetric longitudinal stretching; the N-terminal part dominates this behaviour. Furthermore, alpha B presents a librational motion with respect to the other parts of the molecule with a frequency of approximately 5 cm-1. This motion is coupled to helix stretching. Interestingly, the regions of highly conserved residues contain the most mobile parts of the molecule.     

Inflammation in the avian spleen: timing is everything

Background

Human cells depend critically on the signal recognition particle (SRP) for the sorting and delivery of their proteins. The SRP is a ribonucleoprotein complex which binds to signal sequences of secretory polypeptides as they emerge from the ribosome. Among the six proteins of the eukaryotic SRP, the largest protein, SRP72, is essential for protein targeting and possesses a poorly characterized RNA binding domain.

Results

We delineated the minimal region of SRP72 capable of forming a stable complex with an SRP RNA fragment. The region encompassed residues 545 to 585 of the full-length human SRP72 and contained a lysine-rich cluster (KKKKKKKKGK) at postions 552 to 561 as well as a conserved Pfam motif with the sequence PDPXRWLPXXER at positions 572 to 583. We demonstrated by site-directed mutagenesis that both regions participated in the formation of a complex with the RNA. In agreement with biochemical data and results from chymotryptic digestion experiments, molecular modeling of SRP72 implied that the invariant W577 was located inside the predicted structure of an RNA binding domain. The 11-nucleotide 5e motif contained within the SRP RNA fragment was shown by comparative electrophoresis on native polyacrylamide gels to conform to an RNA kink-turn. The model of the complex suggested that the conserved A240 of the K-turn, previously identified as being essential for the binding to SRP72, could protrude into a groove of the SRP72 RNA binding domain, similar but not identical to how other K-turn recognizing proteins interact with RNA.

Conclusions

The results from the presented experiments provided insights into the molecular details of a functionally important and structurally interesting RNA-protein interaction. A model for how a ligand binding pocket of SRP72 can accommodate a new RNA K-turn in the 5e region of the eukaryotic SRP RNA is proposed.