Crystal structure of an RNA duplex r(gugucgcac)(2) with uridine bulges

Increasing the length of poly-pyrimidine bulges broadens RNA conformational ensembles with minimal impact on stacking energetics

Helical elements separated by bulges frequently undergo transitions between unstacked and coaxially stacked conformations during the folding and function of noncoding RNAs. Here, we examine the dynamic properties of poly-pyrimidine bulges of varying length (n = 1-4, 7) across a range of Mg²⁺ concentrations using HIV-1 TAR RNA as a model system and solution NMR spectroscopy. In the absence of Mg²⁺, helices linked by bulges with n ≥ 3 residues adopt predominantly unstacked conformations (stacked population <15%), whereas one-bulge and two-bulge motifs adopt predominantly stacked conformations (stacked population >74%). In the presence of 3 mM Mg²⁺, the helices predominantly coaxially stack (stacked population >84%), regardless of bulge length, and the midpoint for the Mg²⁺-dependent stacking transition is within threefold regardless of bulge length. In the absence of Mg²⁺, the difference between free energy of interhelical coaxial stacking across the bulge variants is estimated to be ∼2.9 kcal/mol, based on an NMR chemical shift mapping with stacking being more energetically disfavored for the longer bulges. This difference decreases to ∼0.4 kcal/mol in the presence of Mg²⁺ NMR RDCs and resonance intensity data show increased dynamics in the stacked state with increasing bulge length in the presence of Mg²⁺ We propose that Mg²⁺ helps to neutralize the growing electrostatic repulsion in the stacked state with increasing bulge length thereby increasing the number of coaxial conformations that are sampled. Energetically compensated interhelical stacking dynamics may help to maximize the conformational adaptability of RNA and allow a wide range of conformations to be optimally stabilized by proteins and ligands.

An asymmetric mesoscopic model for single bulges in RNA

Simple one-dimensional DNA or RNA mesoscopic models are of interest for their computational efficiency while retaining the key elements of the molecular interactions. However, they only deal with perfectly formed DNA or RNA double helices and consider the intra-strand interactions to be the same on both strands. This makes it difficult to describe highly asymmetric structures such as bulges and loops and, for instance, prevents the application of mesoscopic models to determine RNA secondary structures. Here we derived the conditions for the Peyrard-Bishop mesoscopic model to overcome these limitations and applied it to the calculation of single bulges, the smallest and simplest of these asymmetric structures. We found that these theoretical conditions can indeed be applied to any situation where stacking asymmetry needs to be considered. The full set of parameters for group I RNA bulges was determined from experimental melting temperatures using an optimization procedure, and we also calculated average opening profiles for several RNA sequences. We found that guanosine bulges show the strongest perturbation on their neighboring base pairs, considerably reducing the on-site interactions of their neighboring base pairs.

Distinctive structural motifs of RNA G-quadruplexes composed of AGG, CGG and UGG trinucleotide repeats

Trinucleotide repeats are microsatellite sequences that are polymorphic in length. Their expansion in specific genes underlies a number of neurodegenerative disorders. Using ultraviolet-visible, circular dichroism, nuclear magnetic resonance (NMR) spectroscopies and electrospray ionization mass spectrometry, the structural preferences of RNA molecules composed of two and four repeats of AGG, CGG and UGG in the presence of K⁺, Na⁺ and NH₄⁺ were analysed. (AGG)₂A, (AGG)₄A, p(UGG)₂U and p(UGG)₄U strongly prefer folding into G-quadruplexes, whereas CGG-containing sequences can adopt different types of structure depending on the cation and on the number of repeats. In particular, the two-repeat CGG sequence folds into a G-quadruplex in potassium buffer. We also found that each G-quadruplex fold is different: A:(G:G:G:G)A hexads were found for (AGG)₂A, whereas mixed G:C:G:C tetrads and U-tetrads were observed in the NMR spectra of G(CGG)₂C and p(UGG)₂U, respectively. Finally, our NMR study highlights the influence of the strand sequence on the structure formed, and the influence of the intracellular environment on the folding. Importantly, we highlight that although potassium ions are prevalent in cells, the structures observed in the HeLa cell extract are not always the same as those prevailing in biophysical studies in the presence of K⁺ ions.

Non-nearest-neighbor dependence of the stability for RNA bulge loops based on the complete set of group I single-nucleotide bulge loops

Fifty-nine RNA duplexes containing single-nucleotide bulge loops were optically melted in 1 M NaCl, and the thermodynamic parameters DeltaH degrees, DeltaS degrees, DeltaG 37 degrees, and TM for each sequence were determined. Sequences from this study were combined with sequences from previous studies [Longfellow, C. E., et al. (1990) Biochemistry 29, 278-285; Znosko, B. M., et al. (2002) Biochemistry 41, 10406-10417], thus examining all possible group I single-nucleotide bulge loop and nearest-neighbor sequence combinations. The free energy increments at 37 degrees C for the introduction of a group I single-nucleotide bulge loop range between 1.3 and 5.2 kcal/mol. The combined data were used to develop a model for predicting the free energy of a RNA duplex containing a single-nucleotide bulge. For bulge loops with adjacent Watson-Crick base pairs, neither the identity of the bulge nor the nearest-neighbor base pairs had an effect on the influence of the bulge loop on duplex stability. The proposed model for prediction of the stability of a duplex containing a bulged nucleotide was primarily affected by non-nearest-neighbor interactions. The destabilization of the duplex by the bulge was related to the stability of the stems adjacent to the bulge. Specifically, there was a direct correlation between the destabilization of the duplex and the stability of the less stable duplex stem. The stability of a duplex containing a bulged nucleotide adjacent to a wobble base pair also was primarily affected by non-nearest-neighbor interactions. Again, there was a direct correlation between the destabilization of the duplex and the stability of the less stable duplex stem. However, when one or both of the bulge nearest neighbors was a wobble base pair, the free energy increment for insertion of a bulge loop is dependent upon the position and orientation of the wobble base pair relative the bulged nucleotide. Bulge sequences of the type ((5'UBX)(3'GY)), ((5'GBG)(3'UU)) and ((5'UBU)(3'GG)) are less destabilizing by 0.6 kcal/mol, and bulge sequences of the type ((5'GBX)(3'UY)) and ((5'XBU)(3'YG)) are more destabilizing by 0.4 kcal/mol than bulge loops adjacent to Watson-Crick base pairs.

Base-tetrad swapping results in dimerization of RNA quadruplexes: implications for formation of the i-motif RNA octaplex

Nucleic acids adopt different multistranded helical architectures to perform various biological functions. Here, we report a crystal structure of an RNA quadruplex containing "base-tetrad swapping" and bulged nucleotide at 2.1-Angstroms resolution. The base-tetrad swapping results in a dimer of quadruplexes with an intercalated octaplex fragment at the 5' end junction. The intercalated base tetrads provide the basic repeat unit for constructing a model of intercalated RNA octaplex. The model we obtained shows fundamentally different characteristics from duplex, triplex, and quadruplex. We also observed two different orientations of bulged uridine residues that are related to the interaction with surroundings. This structural evidence reflects the conformational flexibility of bulged nucleotides in RNA quadruplexes and implies the potential roles of bulged nucleotides as recognition and interaction sites in RNA-protein and RNA-RNA interactions.

The prediction of the wild-type telomerase RNA pseudoknot structure and the pivotal role of the bulge in its formation

In this study, the three-dimensional structure of the wild-type human telomerase RNA pseudoknot was predicted via molecular modeling. The wild-type pseudoknot structure is then compared to the recent NMR solution structure of the telomerase pseudoknot, which does not contain the U177 bulge. The removal of the bulge from the pseudoknot structure results in higher stability and significant reduction of activity of telomerase. We show that the effect of the bulge on the structure results in a significant transformation of the pseudoknot junction region where the starting base pairs are disrupted and unique triple base pairs are formed. We found that the formation of the junction region is greatly influenced by interactions of the U177 bulge with loop residues and rotation of residue A174. Moreover, this is the first study to our knowledge where a structure as complex as the pseudoknot has been solved by purely theoretical methods.

RNA binding and thiolytic stability of a quinoline-containing helix-threading peptide

Helix-threading peptides (HTPs) bind selectively to sites predisposed to intercalation in folded RNA molecules placing peptide functional groups into the dissimilar grooves of the duplex. Here we report the design and synthesis of new HTPs with quinoline as the intercalation domain. A quinoline-containing HTP is shown to bind selectively to duplex RNA binding sites. Furthermore, the affinity cleavage pattern generated using an EDTA.Fe modified derivative is consistent with minor groove localization of its N-terminus. This compound binds base-pair steps flanked by single nucleotide bulges on the 3' side on both strands, whereas bulges on the 5' side of the intercalation site do not support binding. Furthermore, unlike acridine HTPs, the quinoline compound is resistant to thiolytic degradation that leads to loss of RNA-binding activity. The RNA-binding selectivity and stability observed for quinoline-containing HTPs make them excellent candidates for further development as regulators of intracellular RNA function.

Crystal structure of a bulged RNA tetraplex at 1.1 a resolution: implications for a novel binding site in RNA tetraplex

Bulges are an important structural motif in RNA and can be used as recognition and interaction sites in RNA-protein interaction and RNA-RNA interaction. Here we report the first crystal structure of a bulged RNA tetraplex at 1.1 A resolution. The hexamer r(U)(BrdG)r(UGGU) forms a parallel tetraplex with the uridine sandwiched by guanines bulging out. The bulged uridine adopts the syn glycosidic conformation and its O2 and N3 atoms face outwards, serving as an effective recognition and interaction site. The bulge formation both widens the groove width and changes the groove hydrogen-bonding pattern on its 5' side. However, the bulge does not make any bends or kinks in the tetraplex structure. The present study demonstrates the dramatic difference between uridine and guanine in forming tetraplex structure. In addition, both G(syn) tetrad and G(anti) tetrad have been observed. They display the same base-pairing pattern and similar C1'-C1' distance but different hydrogen-bonding patterns in the groove.

Solution structures of stem-loop RNAs that bind to the two N-terminal RNA-binding domains of nucleolin

Nucleolin, a multi-domain protein involved in ribosome biogenesis, has been shown to bind the consensus sequence (U/G)CCCG(A/G) in the context of a hairpin loop structure (nucleolin recognition element; NRE). Previous studies have shown that the first two RNA-binding domains in nucleolin (RBD12) are responsible for the interaction with the in vitro selected NRE (sNRE). We have previously reported the structures of nucleolin RBD12, sNRE and nucleolin RBD12-sNRE complex. A comparison of free and bound sNRE shows that the NRE loop becomes structured upon binding. From this observation, we hypothesized that the disordered hairpin loop of sNRE facilitates conformational rearrangements when the protein binds. Here, we show that nucleolin RBD12 is also sufficient for sequence- specific binding of two NRE sequences found in pre-rRNA, b1NRE and b2NRE. Structural investigations of the free NREs using NMR spectroscopy show that the b1NRE loop is conformationally heterogeneous, while the b2NRE loop is structured. The b2NRE forms a hairpin capped by a YNMG-like tetraloop. Comparison of the chemical shifts of sNRE and b2NRE in complex with nucleolin RBD12 suggests that the NRE consensus nucleotides adopt a similar conformation. These results show that a disordered NRE consensus sequence is not a prerequisite for nucleolin RBD12 binding.

Global mapping of nucleic acid conformational space: dinucleoside monophosphate conformations and transition pathways among conformational classes

A global conformational space of 6253 dinucleoside monophosphate (DMP) units consisting of RNA and DNA (free and protein/drug-bound) was 'mapped' using high resolution crystal structures cataloged in the Nucleic Acid Database (NDB). The torsion angles of each DMP were clustered in a reduced three-dimensional space using a classical multi-dimensional scaling method. The mapping of the conformational space reveals nine primary clusters which distinguish among the common A-, B- and Z-forms and their various substates, plus five secondary clusters for kinked or bent structures. Conformational relationships and possible transitional pathways among the substates are also examined using the conformational states of DNA and RNA bound with proteins or drugs as potential pathway intermediates.

The Mg2+ binding sites of the 5S rRNA loop E motif as investigated by molecular dynamics simulations

Molecular dynamics simulations have been used to investigate the binding of Mg(2+) ions to the deep groove of the eubacterial 5S rRNA loop E. The simulations suggest that long-lived and specific water-mediated interactions established between the hydrated ions and the RNA atoms lining up the binding sites contribute to the stabilization of this motif. The Mg(2+) binding specificity is modulated by two factors: (i) a required electrostatic complementarity and (ii) a structural correspondence between the hydrated ion and its binding pocket that can be estimated by its degree of dehydration and the resulting number and lifetime of the intervening water-mediated contacts. Two distinct binding modes for pentahydrated Mg(2+) ions that result in a significant freezing of the tumbling motions of the ions are described, and mechanistic details related to the stabilization of nucleic acids by divalent ions are provided.