Base Pairs In Helix Research Articles

Crystal structure of a poly(rA) staggered zipper at acidic pH: evidence that adenine N1 protonation mediates parallel double helix formation

We have solved at 1.07 Å resolution the X-ray crystal structure of a polyriboadenylic acid (poly(rA)) parallel and continuous double helix. Fifty-nine years ago, double helices of poly(rA) were first proposed to form at acidic pH. Here, we show that 7-mer oligo(rA), i.e. rA7, hybridizes and overlaps in all registers at pH 3.5 to form stacked double helices that span the crystal. Under these conditions, rA7 forms well-ordered crystals, whereas rA6 forms fragile crystalline-like structures, and rA5, rA8 and rA11 fail to crystallize. Our findings support studies from ∼50 years ago: one showed using spectroscopic methods that duplex formation at pH 4.5 largely starts with rA7 and begins to plateau with rA8; another proposed a so-called ‘staggered zipper’ model in which oligo(rA) strands overlap in multiple registers to extend the helical duplex. While never shown, protonation of adenines at position N1 has been hypothesized to be critical for helix formation. Bond angles in our structure suggest that N1 is protonated on the adenines of every other rAMP−rAMP helix base pair. Our data offer new insights into poly(rA) duplex formation that may be useful in developing a pH sensor.

Nuclear Magnetic Resonance-Assisted Prediction of Secondary Structure for RNA: Incorporation of Direction-Dependent Chemical Shift Constraints.

Knowledge of RNA structure is necessary to determine structure–function relationships and to facilitate design of potential therapeutics. RNA secondary structure prediction can be improved by applying constraints from nuclear magnetic resonance (NMR) experiments to a dynamic programming algorithm. Imino proton walks from NOESY spectra reveal double-stranded regions. Chemical shifts of protons in GH1, UH3, and UH5 of GU pairs, UH3, UH5, and AH2 of AU pairs, and GH1 of GC pairs were analyzed to identify constraints for the 5′ to 3′ directionality of base pairs in helices. The 5′ to 3′ directionality constraints were incorporated into an NMR-assisted prediction of secondary structure (NAPSS-CS) program. When it was tested on 18 structures, including nine pseudoknots, the sensitivity and positive predictive value were improved relative to those of three unrestrained programs. The prediction accuracy for the pseudoknots improved the most. The program also facilitates assignment of chemical shifts to individual nucleotides, a necessary step for determining three-dimensional structure.

Sequence and secondary structure of the mitochondrial 16S ribosomal RNA gene of Ixodes scapularis

The complete DNA sequences and secondary structure of the mitochondrial (mt) 16S ribosomal (r) RNA gene were determined for six Ixodes scapularis adults. There were 44 variable nucleotide positions in the 1252 bp sequence alignment. Most (95%) nucleotide alterations did not affect the integrity of the secondary structure of the gene because they either occurred at unpaired positions or represented compensatory changes that maintained the base pairing in helices. A large proportion (75%) of the intraspecific variation in DNA sequence occurred within Domains I, II and VI of the 16S gene. Therefore, several regions within this gene may be highly informative for studies of the population genetics and phylogeography of I. scapularis, a major vector of pathogens of humans and domestic animals in North America.

Identification, molecular characterization, and evolution of group I introns at the expansion segment D11 of 28S rDNA in Rhizoctonia species

The nuclear ribosomal DNA of Rhizoctonia species is polymorphic in terms of the nucleotide composition and length. Insertions of 349–410 nucleotides in length with characteristics of group I introns were detected at a single insertion point at the expansion segment D11 of 28S rDNA in 12 out of 64 isolates. Eleven corresponded to Rhizoctonia solani (teleomorph: Thanatephorous) and one (AG-Q) to Rhizoctonia spp. (teleomorph: Ceratobasidium). Sequence data showed that all but AG-Q contained conserved DNA catalytic core regions (P, Q, R, and S) essential for selfsplicing. The predicted secondary structure revealed that base-paired helices corresponded to subgroup IC1. Isolates from same anastomosis group and even subgroups within R. solani were variable with regard to possession of introns. Phylogenetic analyses indicated that introns were vertically transmitted. Unfortunately, sequence data from the conserved region from all 64 isolates were not useful for delimiting species. Analyses with IC1 introns at same insertion point, of both Ascomycota and Basidiomycota indicated the possibility of horizontal transfer at this site. The present study uncovered new questions on evolutionary pattern of change of these introns within Rhizoctonia species.

Helix untwisting and bubble formation in circular DNA

The base pair fluctuations and helix untwisting are examined for a circular molecule. A realistic mesoscopic model including twisting degrees of freedom and bending of the molecular axis is proposed. The computational method, based on path integral techniques, simulates a distribution of topoisomers with various twist numbers and finds the energetically most favorable molecular conformation as a function of temperature. The method can predict helical repeat, openings loci, and bubble sizes for specific sequences in a broad temperature range. Some results are presented for a short DNA circle recently identified in mammalian cells.

Design of hairpin ribozyme variants with improved activity for poorly processed substrates

Application of ribozymes for knockdown of RNA targets requires the identification of suitable target sites according to the consensus sequence. For the hairpin ribozyme, this was originally defined as Y⁻² N⁻¹ *G+¹ U+² Y+³ B+⁴, with Y = U or C, and B = U, C or G, and C being the preferred nucleobase at positions -2 and +4. In the context of development of ribozymes for destruction of an oncogenic mRNA, we have designed ribozyme variants that efficiently process RNA substrates at U⁻² G⁻¹ *G+¹ U+² A+³ A+⁴ sites. Substrates with G⁻¹ *G+¹ U+² A+³ sites were previously shown to be processed by the wild-type hairpin ribozyme. However, our study demonstrates that, in the specific sequence context of the substrate studied herein, compensatory base changes in the ribozyme improve activity for cleavage (eight-fold) and ligation (100-fold). In particular, we show that A+³ and A+⁴ are well tolerated if compensatory mutations are made at positions 6 and 7 of the ribozyme strand. Adenine at position +4 is neutralized by G⁶ →U, owing to restoration of a Watson-Crick base pair in helix 1. In this ribozyme-substrate complex, adenine at position +3 is also tolerated, with a slightly decreased cleavage rate. Additional substitution of A⁷ with uracil doubled the cleavage rate and restored ligation, which was lost in variants with A⁷, C⁷ and G⁷. The ability to cleave, in conjunction with the inability to ligate RNA, makes these ribozyme variants particularly suitable candidates for RNA destruction.

Variable helix elongation as a tool to modulate RNA alignment and motional couplings

The application of residual dipolar couplings (RDCs) in studies of RNA structure and dynamics can be complicated by the presence of couplings between collective helix motions and overall alignment and by the inability to modulate overall alignment of the molecule by changing the ordering medium. Here, we show for a 27-nt TAR RNA construct that variable levels of helix elongation can be used to alter both overall alignment and couplings to collective helix motions in a semi-predictable manner. In the absence of elongation, a four base-pair helix II capped by a UUCG apical loop exhibits a higher degree of order compared to a six base-pair helix I ( ϑ I / ϑ II = 0.56 ± 0.1 ). The principal S zz direction is nearly parallel to the axis of helix II but deviates by ∼40° relative to the axis of helix I. Elongating helix I by three base-pairs equalizes the alignment of the two helices and pushes the RNA into the motional coupling limit such that the two helices have comparable degrees of order ( ϑ I / ϑ II = 0.92 ± 0.04 ) and orientations relative to S zz (∼17°). Increasing the length of elongation further to 22 base-pairs pushes the RNA into the motional decoupling limit in which helix I dominates alignment ( ϑ II / ϑ I = 0.45 ± 0.05 ), with S zz orientated nearly parallel to its helix axis. Many of these trends can be rationalized using PALES simulations that employ a previously proposed three-state dynamic ensemble of TAR. Our results provide new insights into motional couplings, offer guidelines for assessing their extent, and suggest that variable degrees of helix elongation can allow access to independent sets of RDCs for characterizing RNA structural dynamics.

Novel structural determinants in human SECIS elements modulate the translational recoding of UGA as selenocysteine

The selenocysteine insertion sequence (SECIS) element directs the translational recoding of UGA as selenocysteine. In eukaryotes, the SECIS is located downstream of the UGA codon in the 3′-UTR of the selenoprotein mRNA. Despite poor sequence conservation, all SECIS elements form a similar stem-loop structure containing a putative kink-turn motif. We functionally characterized the 26 SECIS elements encoded in the human genome. Surprisingly, the SECIS elements displayed a wide range of UGA recoding activities, spanning several 1000-fold in vivo and several 100-fold in vitro. The difference in activity between a representative strong and weak SECIS element was not explained by differential binding affinity of SECIS binding Protein 2, a limiting factor for selenocysteine incorporation. Using chimeric SECIS molecules, we identified the internal loop and helix 2, which flank the kink-turn motif, as critical determinants of UGA recoding activity. The simultaneous presence of a GC base pair in helix 2 and a U in the 5′-side of the internal loop was a statistically significant predictor of weak recoding activity. Thus, the SECIS contains intrinsic information that modulates selenocysteine incorporation efficiency.

Control of translation initiation involves a factor‐induced rearrangement of helix 44 of 16S ribosomal RNA

Initiation of translation involves recognition of the start codon by the initiator tRNA in the 30S subunit. To investigate the role of ribosomal RNA (rRNA) in this process, we isolated a number of 16S rRNA mutations that increase translation from the non-canonical start codon AUC. These mutations cluster to distinct regions that overlap remarkably well with previously identified class III protection sites and implicate both IF1 and IF3 in start codon selection. Two mutations map to the 790 loop and presumably act by inhibiting IF3 binding. Another cluster of mutations surrounds the conserved A1413(o)G1487 base pair of helix 44 in a region known to be distorted by IF1 and IF3. Site-directed mutagenesis in this region confirmed that this factor-induced rearrangement of helix 44 helps regulate initiation fidelity. A third cluster of mutations maps to the neck of the 30S subunit, suggesting that the dynamics of the head domain influences translation initiation. In addition to identifying mutations that decrease fidelity, we found that many P-site mutations increase the stringency of start codon selection. These data provide evidence that the interaction between the initiator tRNA and the 30S P site is tuned to balance efficiency and accuracy during initiation.

Genetic identification of yeast 18S rRNA residues required for efficient recruitment of initiator tRNAMet and AUG selection

High-resolution structures of bacterial 70S ribosomes have provided atomic details about mRNA and tRNA binding to the decoding center during elongation, but such information is lacking for preinitiation complexes (PICs). We identified residues in yeast 18S rRNA critical in vivo for recruiting methionyl tRNA(i)(Met) to 40S subunits during initiation by isolating mutations that derepress GCN4 mRNA translation. Several such Gcd(-) mutations alter the A928:U1389 base pair in helix 28 (h28) and allow PICs to scan through the start codons of upstream ORFs that normally repress GCN4 translation. The A928U substitution also impairs TC binding to PICs in a reconstituted system in vitro. Mutation of the bulge G926 in h28 and certain other residues corresponding to direct contacts with the P-site codon or tRNA in bacterial 70S complexes confer Gcd(-) phenotypes that (like A928 substitutions) are suppressed by overexpressing tRNA(i)(Met). Hence, the nonconserved 928:1389 base pair in h28, plus conserved 18S rRNA residues corresponding to P-site contacts in bacterial ribosomes, are critical for efficient Met-tRNA(i)(Met) binding and AUG selection in eukaryotes.

Identification of the RNA Binding Regions of SRP68/72 and SRP72 by Systematic Mutagenesis of Human SRP RNA

Within the large domain of the human signal recognition particle (SRP), 18 mutant SRP RNAs were constructed to disrupt Watson-Crick and G-U basepairs in helices 5, 6, and 8. Using a double-filter assay, the competitive binding of the mutant RNAs to purified human SRP68/72 or to a 7.4 kDa RNA-binding fragment of SRP72 (72frg) was measured. Binding of SRP68/72 was impaired by several mutations in the large domain with the most pronounced effects caused by changes in helix 5 (residues 222-231) and helix 8 (residues 176-191 and 202-214). Binding of the 72frg was diminished prominently by altering helix 5, in particular residues 120-128, and was unaffected by deleting helices 6 and 8. Deleting helix 8 diminished binding of SRP68/72 to a greater extent than deleting helix 6. The data suggest that nucleotide residues throughout most of the large SRP domain are directly and/or indirectly engaged in the binding of SRP68. In contrast, SRP72 binds only to a portion of the 5ef region.

Nudged Elastic Band Calculation of Minimal Energy Paths for the Conformational Change of a GG Non-canonical Pair

The nudged elastic band (NEB) technique has been implemented in AMBER to calculate low-energy paths for conformational changes. A novel simulated annealing protocol that does not require an initial hypothesis for the path is used to sample low-energy paths. This was used to study the conformational change of an RNA cis Watson-Crick/Hoogsteen GG non-canonical pair, with one G syn around the glycosidic bond and the other anti. A previous solution structure, determined by NMR-constrained modeling, demonstrated that the GG pairs change from (syn)G-(anti)G to (anti)G-(syn)G in the context of duplex r(GCAGGCGUGC) on the millisecond timescale. The set of low-energy paths found by NEB show that each G flips independently around the glycosidic bond, with the anti G flipping to syn first. Guanine bases flip without opening adjacent base-pairs by protruding into the major groove, accommodated by a transient change by the ribose to C2'-exo sugar pucker. Hydrogen bonds between bases and the backbone, which lower the energetic barrier to flipping, are observed along the path. The results show the plasticity of RNA base-pairs in helices, which is important for biological processes, including mismatch repair, protein recognition, and translation. The modeling of the GG conformational change also demonstrates that NEB can be used to discover non-trivial paths for macromolecules and therefore NEB can be used as an exploratory method for predicting putative conformational change paths.

Assessing the odd secondary structural properties of nuclear small subunit ribosomal RNA sequences (18S) of the twisted‐wing parasites (Insecta: Strepsiptera)

We report the entire sequence (2864 nts) and secondary structure of the nuclear small subunit ribosomal RNA (SSU rRNA) gene (18S) from the twisted-wing parasite Caenocholax fenyesi texensis Kathirithamby & Johnston (Strepsiptera: Myrmecolacidae). The majority of the base pairings in this structural model map on to the SSU rRNA secondary and tertiary helices that were previously predicted with comparative analysis. These regions of the core rRNA were unambiguously aligned across all Arthropoda. In contrast, many of the variable regions, as previously characterized in other insect taxa, had very large insertions in C. f. texensis. The helical base pairs in these regions were predicted with a comparative analysis of a multiple sequence alignment (that contains C. f. texensis and 174 published arthropod 18S rRNA sequences, including eleven strepsipterans) and thermodynamic-based algorithms. Analysis of our structural alignment revealed four unusual insertions in the core rRNA structure that are unique to animal 18S rRNA and in general agreement with previously proposed insertion sites for strepsipterans. One curious result is the presence of a large insertion within a hairpin loop of a highly conserved pseudoknot helix in variable region 4. Despite the extraordinary variability in sequence length and composition, this insertion contains the conserved sequences 5'-AUUGGCUUAAA-3' and 5'-GAC-3' that immediately flank a putative helix at the 5'- and 3'-ends, respectively. The longer sequence has the potential to form a nine base pair helix with a sequence in the variable region 2, consistent with a recent study proposing this tertiary interaction. Our analysis of a larger set of arthropod 18S rRNA sequences has revealed possible errors in some of the previously published strepsipteran 18S rRNA sequences. Thus we find no support for the previously recovered heterogeneity in the 18S molecules of strepsipterans. Our findings lend insight to the evolution of RNA structure and function and the impact large insertions pose on genome size. We also provide a novel alignment template that will improve the phylogenetic placement of the Strepsiptera among other insect taxa.

Influence of Helix Length on Cleavage Efficiency of Hammerhead Ribozymes

The cleavage rates of RNA substrates by trans-acting, hammerhead ribozymes are controlled by interactions between helices I and II. The interactions are affected by the relative lengths of these two double helices and by unpaired nucleotides protruding beyond helix I, either in the substrate or the ribozyme strand. Maximum cleavage rates are observed for ribozyme–substrate complexes with three or more base pairs in helix II and six or less base pairs in helix I. However, for these helix combinations, rates fall sharply with unpaired nucleotides at the end of helix I. Cleavage rates by ribozymes with one or two base pairs in helix II increase as helix I is lengthened, and are unaffected by unpaired nucleotides on the end. Since miniribozymes, with one base pair in helix II, efficiently cleave long RNA transcripts under physiological conditions, they represent the optimal design for the simple hammerheads for application in vivo.

MRNA Helicase Activity of the Ribosome

Most mRNAs contain secondary structure, yet their codons must be in single-stranded form to be translated. Until now, no helicase activity has been identified which could account for the ability of ribosomes to translate through downstream mRNA secondary structure. Using an oligonucleotide displacement assay, together with a stepwise in vitro translation system made up of purified components, we show that ribosomes are able to disrupt downstream helices, including a perfect 27 base pair helix of predicted T m = 70°. Using helices of different lengths and registers, the helicase active site can be localized to the middle of the downstream tunnel, between the head and shoulder of the 30S subunit. Mutation of residues in proteins S3 and S4 that line the entry to the tunnel impairs helicase activity. We conclude that the ribosome itself is an mRNA helicase and that proteins S3 and S4 may play a role in its processivity.

An emerging consensus for telomerase RNA structure.

Telomerase maintains telomere length by adding G-rich telomeric repeats to the ends of linear eukaryotic chromosomes. The core telomerase enzyme consists of two components: an essential RNA component and a catalytic protein component. The catalytic component, telomerase reverse transcriptase (TERT), contains sequence motifs homologous to those in the catalytic domain of reverse transcriptase enzymes. Although the TERT component of telomerase is fairly conserved among eukaryotes, the RNA component varies dramatically in sequence composition and in size. The structure of the telomerase RNA has been well established in ciliates and vertebrates (1–4). These structural models have led to experiments that provided fundamental insight into telomerase RNA function. However, it has been difficult to determine the global secondary structure of the very large yeast RNA. In the past month, there have been four new reports proposing a secondary structure for the yeast telomerase RNA (5–8). In this issue of PNAS, Blackburn and colleagues (7) propose a core structure for the yeast telomerase RNA that resembles the structure present in ciliate and vertebrate telomerase RNAs. Telomerase RNA is very divergent across eukaryotic species. The size ranges from ≈150 nt in ciliates to >1,300 nt in fungi. Although there is some sequence similarity within each group of closely related eukaryotes, there is no sequence similarity between groups. Thus, the telomerase RNA sequences from vertebrates cannot be aligned with the RNA sequences from ciliates. The rapid divergence in RNA size and sequence has made it difficult to identify common secondary structures for telomerase RNAs. Phylogenetic comparative analysis has proven to be the most powerful approach for inferring higher-order RNA structures. In this method, base paired helices in RNA secondary structure are determined by nucleotide covariation between species (9). For example, an A:U base pair in one species may change to a …

Determination of thermodynamic parameters for HIV DIS type loop-loop kissing complexes.

The HIV-1 type dimerization initiation signal (DIS) loop was used as a starting point for the analysis of the stability of Watson-Crick (WC) base pairs in a tertiary structure context. We used ultraviolet melting to determine thermodynamic parameters for loop-loop tertiary interactions and compared them with regular secondary structure RNA helices of the same sequences. In 1 M Na+ the loop-loop interaction of a HIV-1 DIS type pairing is 4 kcal/mol more stable than its sequence in an equivalent regular and isolated RNA helix. This difference is constant and sequence independent, suggesting that the rules governing the stability of WC base pairs in the secondary structure context are also valid for WC base pairs in the tertiary structure context. Moreover, the effect of ion concentration on the stability of loop-loop tertiary interactions differs considerably from that of regular RNA helices. The stabilization by Na+ and Mg2+ is significantly greater if the base pairing occurs within the context of a loop-loop interaction. The dependence of the structural stability on salt concentration was defined via the slope of a T(m)/log [ion] plot. The short base-paired helices are stabilized by 8 degrees C/log [Mg2+] or 11 degrees C/log [Na+], whereas base-paired helices forming tertiary loop-loop interactions are stabilized by 16 degrees C/log [Mg2+] and 26 degrees C/log [Na+]. The different dependence on ionic strength that is observed might reflect the contribution of specific divalent ion binding to the preformation of the hairpin loops poised for the tertiary kissing loop-loop contacts.

Kinetic analysis of ribozyme-substrate complex formation in yeast.

Many RNA-mediated reactions in transcription, translation, RNA processing, and transport require assembly of RNA complexes, yet assembly pathways remain poorly understood. Assembly mechanisms can be difficult to assess in a biological context because many components interact in complex pathways and individual steps are difficult to isolate experimentally. Our previous studies of self-cleaving hairpin ribozymes showed that kinetic and equilibrium parameters measured in yeast agree well with parameters measured in vitro under ionic conditions that mimic the intracellular environment. We now report studies of intermolecular reactions with ribozyme and target sequences expressed in yeast as separate chimeric U3 snoRNAs. In this system, intracellular cleavage rates reflect the kinetics of ribozyme-substrate complex formation through annealing of base-paired helices. Second-order rate constants increased with increasing helix length for in vitro reactions with 2 mM MgCl(2) and 150 mM NaCl and in vivo but not in reactions with 10 mM MgCl(2). Thus, efficient RNA complex formation required a larger extent of complementarity in vivo than in vitro under conditions with high concentrations of divalent cations. The most efficient intracellular cleavage reactions exhibited second-order rate constants that were 15- to 30-fold below rate constants for cleavage of oligonucleotides in vitro. Careful analysis of structural features that influence cleavage efficiency points to substrate binding as the rate-determining step in the intracellular cleavage pathway. Second-order rate constants for intermolecular cleavage agree well with diffusion coefficients reported for U3 snoRNPs in vivo suggesting that complex formation between chimeric ribozyme and substrate snoRNPs in yeast nuclei is diffusion limited.

Study of the Functional Interaction of the 900 Tetraloop of 16 S Ribosomal RNA with Helix 24 within the Bacterial Ribosome

The 900 tetraloop that caps helix 27 of 16 S ribosomal RNA (rRNA) is amongst the most conserved regions of rRNA. This tetraloop forms a GNRA motif that docks into the minor groove of three base-pairs at the bottom of helix 24 of 16 S rRNA in the 30 S subunit. Both the tetraloop and its receptor in helix 24 contact the 23 S rRNA, forming the intersubunit bridge B2c. Here, we investigated the interaction between the 900 tetraloop and its receptor by genetic complementation. We used a specialized ribosome system in combination with an in vivo instant evolution approach to select mutations in helix 24 compensating for a mutation in the 900 tetraloop (A900G) that severely decreases ribosomal activity, impairing subunit association and translational fidelity. We selected two mutants where the G769-C810 base-pair of helix 24 was substituted with either U-A or C·A. When these mutations in helix 24 were investigated in the context of a wild-type 900 tetraloop, the C·A but not the U-A mutation severely impaired ribosome activity, interfering with subunit association and decreasing translational fidelity. In the presence of the A900G mutation, both mutations in helix 24 increased the ribosome activity to the same extent. Subunit association and translational fidelity were increased to the same level. Computer modeling was used to analyze the effect of the mutations in helix 24 on the interaction between the tetraloop and its receptor. This study demonstrates the functional importance of the interaction between the 900 tetraloop and helix 24.

Two strategically placed base pairs in helix 8 of mammalian signal recognition particle RNA are crucial for the SPR19-dependent binding of protein SRP54

Signal recognition particle (SRP) guides secretory proteins to biological membranes in all organisms. Assembly of the large domain of mammalian SRP requires binding of SRP19 prior to the binding of protein SRP54 to SRP RNA. The crystal structure of the ternary complex reveals the parallel arrangement of RNA helices 6 and 8, a bridging of the helices via a hydrogen bonded A149-A201 pair and protein SRP19, and two A minor motifs between the asymmetric loop of helix 8 (A213 and A214) and helix 6. We investigated which residues in helix 8 are responsible for the SRP19-dependent binding of SRP54 by taking advantage of the finding that binding of human SRP54 to Methanococcus jannaschii SRP RNA is independent of SRP19. Chimeric human/M. jannaschii SRP RNA molecules were synthesized containing predominantly human SRP RNA but possessing M. jannaschii SRP RNA-derived substitutions. Activities of the chimeric RNAs were measured with respect to protein SRP19 and the methionine-rich RNA-binding domain of protein SRP54 (SRP54M). Changing A213 and A214 to a uridine has no effect on the SRP19-dependent binding of SRP54M. Instead, the two base pairs C189-G210 and C190-G209, positioned between the conserved binding site of SRP54 and the asymmetric loop, are critical for conveying SRP19 dependency. Furthermore, the nucleotide composition of five base pairs surrounding the asymmetric loop affects binding of SRP54M significantly. These results demonstrate that subtle, and not easily perceived, structural differences are of crucial importance in the assembly of mammalian SRP.