Hinge-like motions in RNA kink-turns: the role of the second a-minor motif and nominally unpaired bases

Kink-turn (K-turn) motifs are asymmetric internal loops found at conserved positions in diverse RNAs, with sharp bends in phosphodiester backbones producing V-shaped structures. Explicit-solvent molecular dynamics simulations were carried out for three K-turns from 23S rRNA, i.e., Kt-38 located at the base of the A-site finger, Kt-42 located at the base of the L7/L12 stalk, and Kt-58 located in domain III, and for the K-turn of human U4 snRNA. The simulations reveal hinge-like K-turn motions on the nanosecond timescale. The first conserved A-minor interaction between the K-turn stems is entirely stable in all simulations. The angle between the helical arms of Kt-38 and Kt-42 is regulated by local variations of the second A-minor (type I) interaction between the stems. Its variability ranges from closed geometries to open ones stabilized by insertion of long-residency waters between adenine and cytosine. The simulated A-minor geometries fully agree with x-ray data. Kt-58 and Kt-U4 exhibit similar elbow-like motions caused by conformational change of the adenosine from the nominally unpaired region. Despite the observed substantial dynamics of K-turns, key tertiary interactions are stable and no sign of unfolding is seen. We suggest that some K-turns are flexible elements mediating large-scale ribosomal motions during the protein synthesis cycle.     

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.     

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.     

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.     

Biochemical characterization of the kink-turn RNA motif

RNA, which acts as a medium for transmitting genetic information, plays a variety of roles in a cell. As with proteins, elucidation of the three- dimensional (3D) structures of RNAs is important for understanding their various roles. Determination of the atomic structures of crystallized ribosome has enabled the identification of previously unknown RNA structural motifs. The kink-turn (K-turn or GA) motif, which causes a sharp bend in an RNA double helix, was identified as one of these structural motifs. To biochemically characterize the K-turn, the motif was inserted into a hinge region of P4-P6 RNA, which is the most extensively studied self-folding RNA, and its properties were investigated. The stability and metal ion requirement of the constructs containing three different K-turn motifs were analyzed using native PAGE and dimethyl sulfate (DMS) modification. The formation of the sharp bending structure depends on the presence of divalent cation like Mg²⁺ or Ca²⁺, although its required concentration is different for each motif.     

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.     

The kink-turn motif in RNA is dimorphic, and metal ion-dependent

The kink-turn (K-turn) is a new motif in RNA structure that was identified by examination of the crystal structures of the ribosome. We examined the structural and dynamic properties of this element in free solution. The K-turn RNA exists in a dynamic equilibrium between a tightly kinked conformation and a more open structure similar to a simple bulge bend. The highly kinked form is stabilized by the noncooperative binding of metal ions, but a significant population of the less-kinked form is present even in the presence of relatively high concentrations of divalent metal ions. The conformation of the tightly kinked population is in excellent agreement with that of the K-turn structures observed in the ribosome by crystallography. The end-to-end FRET efficiency of this species agrees closely with that of the ribosomal K-turn, and the direction of the bend measured by comparative gel electrophoresis also corresponds very well. These results show that the tightly kinked conformation of the K-turn requires stabilization by other factors, possibly by protein binding, for example. The K-turn is therefore unlikely to be of itself a primary organizing feature in RNA.     

The Contribution of a Covalently Bound Cofactor to the Folding and Thermodynamic Stability of an Integral Membrane Protein

The factors controlling the stability, folding, and dynamics of integral membrane proteins are not fully understood. The high stability of the membrane protein bacteriorhodopsin (bR), an archetypal member of the rhodopsin photoreceptor family, has been ascribed to its covalently bound retinal cofactor. We investigate here the role of this cofactor in the thermodynamic stability and folding kinetics of bR. Multiple spectroscopic probes were used to determine the kinetics and energetics of protein folding in mixed lipid/detergent micelles in the presence and absence of retinal. The presence of retinal increases extrapolated values for the overall unfolding free energy from 6.3 ± 0.4 kcal mol^− 1 to 23.4 ± 1.5 kcal mol^− 1 at zero denaturant, suggesting that the cofactor contributes 17.1 kcal mol^− 1 towards the overall stability of bR. In addition, the cooperativity of equilibrium unfolding curves is markedly reduced in the absence of retinal with overall m-values decreasing from 31.0 ± 2.0 kcal mol^− 1 to 10.9 ± 1.0 kcal mol^− 1, indicating that the folded state of the apoprotein is less compact than the equivalent for the holoprotein. This change in the denaturant response means that the difference in the unfolding free energy at a denaturant concentration midway between the two unfolding curves is only ca 3-6 kcal mol^− 1. Kinetic data show that the decrease in stability upon removal of retinal is associated with an increase in the apparent intrinsic rate constant of unfolding, k_u^H2O, from ~1 × 10^− 16 s^− 1 to ~1 × 10^− 4 s^− 1 at 25 °C. This correlates with a decrease in the unfolding activation energy by 16.3 kcal mol^− 1 in the apoprotein, extrapolated to zero SDS. These results suggest that changes in bR stability induced by retinal binding are mediated solely by changes in the activation barrier for unfolding. The results are consistent with a model in which bR is kinetically stabilized via a very slow rate of unfolding arising from protein-retinal interactions that increase the rigidity and compactness of the polypeptide chain.     

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.     

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.     

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.     

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.     

SEA domain autoproteolysis accelerated by conformational strain: energetic aspects

A subclass of proteins with the SEA (sea urchin sperm protein, enterokinase, and agrin) domain fold exists as heterodimers generated by autoproteolytic cleavage within a characteristic G^− 1S^+ 1VVV sequence. Autoproteolysis occurs by a nucleophilic attack of the serine hydroxyl on the vicinal glycine carbonyl followed by an N → O acyl shift and hydrolysis of the resulting ester. The reaction has been suggested to be accelerated by the straining of the scissile peptide bond upon protein folding. In an accompanying article, we report the mechanism; in this article, we provide further key evidence and account for the energetics of coupled protein folding and autoproteolysis. Cleavage of the GPR116 domain and that of the MUC1 SEA domain occur with half-life (t_½) values of 12 and 18 min, respectively, with lowering of the free energy of the activation barrier by ∼ 10 kcal mol^− 1 compared with uncatalyzed hydrolysis. The free energies of unfolding of the GPR116 and MUC1 SEA domains were measured to ∼ 11 and ∼ 15 kcal mol^− 1, respectively, but ∼ 7 kcal mol^− 1 of conformational energy is partitioned as strain over the scissile peptide bond in the precursor to catalyze autoproteolysis by substrate destabilization. A straining energy of ∼ 7 kcal mol^− 1 was measured by using both a pre-equilibrium model to analyze stability and cleavage kinetics data obtained with the GPR116 SEA domain destabilized by core mutations or urea addition, as well as the difference in thermodynamic stabilities of the MUC1 SEA precursor mutant S1098A (with a G^− 1A^+ 1VVV motif) and the wild-type protein. The results imply that cleavage by N → O acyl shift alone would proceed with a t_½ of ∼ 2.3 years, which is too slow to be biochemically effective. A subsequent review of structural data on other self-cleaving proteins suggests that conformational strain of the scissile peptide bond may be a common mechanism of autoproteolysis.     

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.     

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.     

Unfolding thermodynamics of the Delta-domain in the prohead I subunit of phage HK97: determination by factor analysis of Raman spectra

An early step in the morphogenesis of the double-stranded DNA (dsDNA) bacteriophage HK97 is the assembly of a precursor shell (prohead I) from 420 copies of a 384-residue subunit (gp5). Although formation of prohead I requires direct participation of gp5 residues 2-103 (Δ-domain), this domain is eliminated by viral protease prior to subsequent shell maturation and DNA packaging. The prohead I Δ-domain is thought to resemble a phage scaffolding protein, by virtue of its highly α-helical secondary structure and a tertiary fold that projects inward from the interior surface of the shell. Here, we employ factor analysis of temperature-dependent Raman spectra to characterize the thermostability of the Δ-domain secondary structure and to quantify the thermodynamic parameters of Δ-domain unfolding. The results are compared for the Δ-domain within the prohead I architecture (in situ) and for a recombinantly expressed 111-residue peptide (in vitro). We find that the α-helicity (∼ 70%), median melting temperature (T_m = 58 °C), enthalpy (ΔH_m = 50 ± 5 kcal mol^− 1), entropy (ΔS_m = 150 ± 10 cal mol^− 1 K^− 1), and average cooperative melting unit (〈n_c〉 ∼ 3.5) of the in situ Δ-domain are altered in vitro, indicating specific interdomain interactions within prohead I. Thus, the in vitro Δ-domain, despite an enhanced helical secondary structure (∼ 90% α-helix), exhibits diminished thermostability (T_m = 40 °C; ΔH_m = 27 ± 2 kcal mol^− 1; ΔS_m = 86 ± 6 cal mol^− 1 K^− 1) and noncooperative unfolding (〈n_c〉 ∼ 1) vis-à-vis the in situ Δ-domain. Temperature-dependent Raman markers of subunit side chains, particularly those of Phe and Trp residues, also confirm different local interactions for the in situ and in vitro Δ-domains. The present results clarify the key role of the gp5 Δ-domain in prohead I architecture by providing direct evidence of domain structure stabilization and interdomain interactions within the assembled shell.     

The energetic contribution of induced electrostatic asymmetry to DNA bending by a site-specific protein

DNA bending can be promoted by reducing the net negative electrostatic potential around phosphates on one face of the DNA, such that electrostatic repulsion among phosphates on the opposite face drives bending toward the less negative surface. To provide the first assessment of energetic contribution to DNA bending when electrostatic asymmetry is induced by a site-specific DNA binding protein, we manipulated the electrostatics in the EcoRV endonuclease-DNA complex by mutation of cationic side chains that contact DNA phosphates and/or by replacement of a selected phosphate in each strand with uncharged methylphosphonate. Reducing the net negative charge at two symmetrically located phosphates on the concave DNA face contributes − 2.3 kcal mol^− 1 to − 0.9 kcal mol^− 1 (depending on position) to complex formation. In contrast, reducing negative charge on the opposing convex face produces a penalty of + 1.3 kcal mol^− 1. Förster resonance energy transfer experiments show that the extent of axial DNA bending (about 50°) is little affected in modified complexes, implying that modification affects the energetic cost but not the extent of DNA bending. Kinetic studies show that the favorable effects of induced electrostatic asymmetry on equilibrium binding derive primarily from a reduced rate of complex dissociation, suggesting stabilization of the specific complex between protein and markedly bent DNA. A smaller increase in the association rate may suggest that the DNA in the initial encounter complex is mildly bent. The data imply that protein-induced electrostatic asymmetry makes a significant contribution to DNA bending but is not itself sufficient to drive full bending in the specific EcoRV-DNA complex.