Preparation of selective and segmentally labeled single-stranded DNA for NMR by self-primed PCR and asymmetrical endonuclease double digestion

We demonstrate a new, efficient and easy-to-use method for enzymatic synthesis of (stereo-)specific and segmental ¹³C/¹⁵N/²H isotope-labeled single-stranded DNA in amounts sufficient for NMR, based on the highly efficient self-primed PCR. To achieve this, new approaches are introduced and combined. (i) Asymmetric endonuclease double digestion of tandem-repeated PCR product. (ii) T4 DNA ligase mediated ligation of two ssDNA segments. (iii) In vitro dNTP synthesis, consisting of in vitro rNTP synthesis followed by enzymatic stereo-selective reduction of the C2′ of the rNTP, and a one-pot add-up synthesis of dTTP from dUTP. The method is demonstrated on two ssDNAs: (i) a 36-nt three-way junction, selectively ¹³C₉/¹⁵N₃/²H_{(1′,2″,3′,4′,5′,5″)}-dC labeled and (ii) a 39-nt triple-repeat three-way junction, selectively ¹³C₉/¹⁵N₃/²H_{(1′,2″,3′,4′,5′,5″)}-dC and ¹³C₉/¹⁵N₂/²H_{(1′,2″,3′,4′,5′,5″)}-dT labeled in segment C20-C39. Their NMR spectra show the spectral simplification, while the stereo-selective ²H-labeling in the deoxyribose of the dC-residues, straightforwardly provided assignment of their C1′–H2′ and C2′–H2′ resonances. The labeling protocols can be extended to larger ssDNA molecules and to more than two segments.     

Structural basis for uracil DNA glycosylase interaction with uracil: NMR study

Two dimensional (2D) NMR and molecular dynamics simulations have been used to determine the three dimensional (3D) structure of a hairpin DNA, d-CTA-GAGGATCC-TUTT-GGATCCT (22mer; abbreviated as U2-hairpin), which has uracil at the second position from the 5′ end of the tetraloop. The ¹H resonances of this hairpin have been assigned almost completely. NMR restrained molecular dynamics and energy minimization procedures have been used to describe the 3D structure of U2-hairpin. This study establishes that the stem of the hairpin adopts a right-handed B-DNA conformation, while the T₁₂ and T₁₅ nucleotides stack upon 3′ and 5′ ends of the stem, respectively. Further, T₁₄ stacks upon both T₁₂ and T₁₅. Though U₁₃ partially stacks upon T₁₄, no stacking interaction is observed between U₁₃ and T₁₂. All the individual nucleotide bases belonging to the stem and T₁₂ and T₁₅ of the loop adopt ‘anti’ conformation with respect to their sugar moiety, while the U₁₃ and T₁₄ of the loop are in ‘syn’ conformation. The turning phosphate in the loop is located between T₁₃ and T₁₄. This study and a concurrent NMR structural study on yet another hairpin DNA d-CTAGAGGAATAA-TTTU-GGATCCT (22mer; abbreviated as U4-hairpin), with uracil at the fourth position from the 5′ end of the tetraloop throw light upon various interactions which have been reported between Escherichia coli uracil DNA glycosylase (UDG) and uracil containing DNA. The of T₁₂ and α, β, γ, and ζ of U₁₃ and γ of T₁₄, which partially influence the local conformation of U₁₃ in U2-hairpin are all locked in ‘trans’ conformation. Such stretched out backbone conformation in the vicinity of U₁₃ could be the reason as to why the U2-hairpin is found to be the poor substrate for its interaction with UDG compared to the other substrates in which the uracil is at first, third and fourth positions of the tetraloop from its 5′ end, as reported earlier by Vinay and Varshney. This study shows that UDG actively promotes the flipping of uracil from a stacked conformation and rules out the possibility of UDG recognizing the flipped out uracil bases.     

Unique structural and stabilizing roles for the individual pseudouridine residues in the 1920 region of Escherichia coli 23S rRNA

The synthesis of a 5′-O-BzH–2′-O-ACE-protected pseudouridine phosphoramidite is reported [BzH, benzhydryloxy-bis(trimethylsilyloxy)silyl; ACE, bis(2-acetoxyethoxy)methyl]. The availability of the phosphoramidite allows for reliable and efficient syntheses of hairpin RNAs containing single or multiple pseudouridine modifications in the stem or loop regions. Five 19-nt hairpin RNAs representing the 1920-loop region (G₁₉₀₆–C₁₉₂₄) of Escherichia coli 23S rRNA were synthesized with pseudouridine residues located at positions 1911, 1915 and 1917. Thermodynamic parameters, circular dichroism spectra and NMR data are presented for all five RNAs. Overall, three different structural contexts for the pseudouridine residues were examined and compared with the unmodified RNA. Our main findings are that pseudouridine modifications exhibit a range of effects on RNA stability and structure, depending on their locations. More specifically, pseudouridines in the single-stranded loop regions of the model RNAs are slightly destabilizing, whereas a pseudouridine at the stem–loop junction is stabilizing. Furthermore, the observed effects on stability are approximately additive when multiple pseudouridine residues are present. The possible relationship of these results to RNA function is discussed.     

The determination of tRNALeu recognition nucleotides for Escherichia coli L/F transferase

Escherichia coli leucyl/phenylalanyl-tRNA protein transferase catalyzes the tRNA-dependent post-translational addition of amino acids onto the N-terminus of a protein polypeptide substrate. Based on biochemical and structural studies, the current tRNA recognition model by L/F transferase involves the identity of the 3′ aminoacyl adenosine and the sequence-independent docking of the D-stem of an aminoacyl-tRNA to the positively charged cluster on L/F transferase. However, this model does not explain the isoacceptor preference observed 40 yr ago. Using in vitro-transcribed tRNA and quantitative MALDI-ToF MS enzyme activity assays, we have confirmed that, indeed, there is a strong preference for the most abundant leucyl-tRNA, tRNA^Leu (anticodon 5′-CAG-3′) isoacceptor for L/F transferase activity. We further investigate the molecular mechanism for this preference using hybrid tRNA constructs. We identified two independent sequence elements in the acceptor stem of tRNA^Leu (CAG)—a G₃:C₇₀ base pair and a set of 4 nt (C₇₂, A₄:U₆₉, C₆₈)—that are important for the optimal binding and catalysis by L/F transferase. This maps a more specific, sequence-dependent tRNA recognition model of L/F transferase than previously proposed.     

Structural domains within the 3' untranslated region of Turnip crinkle virus

The genomes of positive-strand RNA viruses undergo conformational shifts that complicate efforts to equate structures with function. We have initiated a detailed analysis of secondary and tertiary elements within the 3′ end of Turnip crinkle virus (TCV) that are required for viral accumulation in vivo. MPGAfold, a massively parallel genetic algorithm, suggested the presence of five hairpins (H4a, H4b, and previously identified hairpins H4, H5, and Pr) and one H-type pseudoknot (Ψ₃) within the 3′-terminal 194 nucleotides (nt). In vivo compensatory mutagenesis analyses confirmed the existence of H4a, H4b, Ψ₃ and a second pseudoknot (Ψ₂) previously identified in a TCV satellite RNA. In-line structure probing of the 194-nt fragment supported the coexistence of H4, H4a, H4b, Ψ₃ and a pseudoknot that connects H5 and the 3′ end (Ψ₁). Stepwise replacements of TCV elements with the comparable elements from Cardamine chlorotic fleck virus indicated that the complete 142-nt 3′ end, and subsets containing Ψ₃, H4a, and H4b or Ψ₃, H4a, H4b, H5, and Ψ₂, form functional domains for virus accumulation in vivo. A new 3-D molecular modeling protocol (RNA2D3D) predicted that H4a, H4b, H5, Ψ₃, and Ψ₂ are capable of simultaneous existence and bears some resemblance to a tRNA. The related Japanese iris necrotic ring virus does not have comparable domains. These results provide a framework for determining how interconnected elements participate in processes that require 3′ untranslated region sequences such as translation and replication.     

Reproducible features of small RNAs in C. elegans reveal NU RNAs and provide insights into 22G RNAs and 26G RNAs

Small RNAs regulate gene expression and most genes in the worm Caenorhabditis elegans are subject to their regulation. Here, we analyze small RNA data sets and use reproducible features of RNAs present in multiple data sets to discover a new class of small RNAs and to reveal insights into two known classes of small RNAs—22G RNAs and 26G RNAs. We found that reproducibly detected 22-nt RNAs, although are predominantly RNAs with a G at the 5′ end, also include RNAs with A, C, or U at the 5′ end. These RNAs are synthesized downstream from characteristic sequence motifs on mRNA and have U-tailed derivatives. Analysis of 26G RNAs revealed that they are processed from a blunt end of double-stranded RNAs and that production of one 26G RNA generates a hotspot immediately downstream for production of another. To our surprise, analysis of RNAs shorter than 18 nt revealed a new class of RNAs, which we call NU RNAs (pronounced “new RNAs”) because they have a NU bias at the 5′ end, where N is any nucleotide. NU RNAs are antisense to genes and originate downstream from U bases on mRNA. Although many genes have complementary NU RNAs, their genome-wide distribution is distinct from that of previously known classes of small RNAs. Our results suggest that current approaches underestimate reproducibly detected RNAs that are shorter than 18 nt, and theoretical considerations suggest that such shorter RNAs could be used for sequence-specific gene regulation in organisms like C. elegans that have small genomes.     

Locating the uracil-5-yl radical formed upon photoirradiation of 5-bromouracil-substituted DNA

In a previous study, we found that 2-deoxyribonolactone is effectively generated in the specific 5-bromouracil (^BrU)-substituted sequence 5′-(G/C)[A]_{n = 1,2}^BrU^BrU-3′ and proposed that a formed uracil-5-yl radical mainly abstracts the C1′ hydrogen from the 5′-side of ^BrU^BrU under 302-nm irradiation condition. In the present work, we performed photoirradiation of ^BrU-substituted DNA in the presence of a hydrogen donor, tetrahydrofuran, to quench the uracil-5-yl radical to uracil and then subjected the sample to uracil DNA glycosylase digestion. Slab gel sequence analysis indicated that uracil residues were formed at the hot-spot sequence of 5′-(G/C)[A]_{n = 1,2}^BrU^BrU-3′ in 302-nm irradiation of ^BrU-substituted DNA. Furthermore, we found that the uracil residue was also formed at the reverse sequence 5′-^BrU^BrU[A]_{n = 1,2}(G/C)-3′, which suggests that both 5′-(G/C)[A]_{n = 1,2}^BrU^BrU-3′ and 5′-^BrU^BrU[A]_{n = 1,2}(G/C)-3′ are hot-spot sequences for the formation of the uracil-5-yl radical.     

Inhibition of translation by IFIT family members is determined by their ability to interact selectively with the 5′-terminal regions of cap0-, cap1- and 5′ppp- mRNAs

Ribosomal recruitment of cellular mRNAs depends on binding of eIF4F to the mRNA’s 5′-terminal ‘cap’. The minimal ‘cap0’ consists of N7-methylguanosine linked to the first nucleotide via a 5′-5′ triphosphate (ppp) bridge. Cap0 is further modified by 2′-O-methylation of the next two riboses, yielding ‘cap1’ (m7GpppNmN) and ‘cap2’ (m7GpppNmNm). However, some viral RNAs lack 2′-O-methylation, whereas others contain only ppp- at their 5′-end. Interferon-induced proteins with tetratricopeptide repeats (IFITs) are highly expressed effectors of innate immunity that inhibit viral replication by incompletely understood mechanisms. Here, we investigated the ability of IFIT family members to interact with cap1-, cap0- and 5′ppp- mRNAs and inhibit their translation. IFIT1 and IFIT1B showed very high affinity to cap-proximal regions of cap0-mRNAs (K_1/2,app ∼9 to 23 nM). The 2′-O-methylation abrogated IFIT1/mRNA interaction, whereas IFIT1B retained the ability to bind cap1-mRNA, albeit with reduced affinity (K_1/2,app ∼450 nM). The 5′-terminal regions of 5′ppp-mRNAs were recognized by IFIT5 (K_1/2,app ∼400 nM). The activity of individual IFITs in inhibiting initiation on a specific mRNA was determined by their ability to interact with its 5′-terminal region: IFIT1 and IFIT1B efficiently outcompeted eIF4F and abrogated initiation on cap0-mRNAs, whereas inhibition on cap1- and 5′ppp- mRNAs by IFIT1B and IFIT5 was weaker and required higher protein concentrations.     

Regulation of small RNA stability: methylation and beyond

As central components of RNA silencing, small RNAs play diverse and important roles in many biological processes in eukaryotes. Aberrant reduction or elevation in the levels of small RNAs is associated with many developmental and physiological defects. The in vivo levels of small RNAs are precisely regulated through modulating the rates of their biogenesis and turnover. 2′-O-methylation on the 3′ terminal ribose is a major mechanism that increases the stability of small RNAs. The small RNA methyltransferase HUA ENHANCER1 (HEN1) and its homologs methylate microRNAs and small interfering RNAs (siRNAs) in plants, Piwi-interacting RNAs (piRNAs) in animals, and siRNAs in Drosophila. 3′ nucleotide addition, especially uridylation, and 3′-5′ exonucleolytic degradation are major mechanisms that turnover small RNAs. Other mechanisms impacting small RNA stability include complementary RNAs, cis-elements in small RNA sequences and RNA-binding proteins. Investigations are ongoing to further understand how small RNA stability impacts their accumulation in vivo in order to improve the utilization of RNA silencing in biotechnology and therapeutic applications.     

siRNAs compete with miRNAs for methylation by HEN1 in Arabidopsis

Plant microRNAs (miRNAs) and small interfering RNAs (siRNAs) bear a 2′-O-methyl group on the 3′-terminal nucleotide. This methyl group is post-synthetically added by the methyltransferase protein HEN1 and protects small RNAs from enzymatic activities that target the 3′-OH. A mutagenesis screen for suppressors of the partial loss-of-function hen1-2 allele in Arabidopsis identified second-site mutations that restore miRNA methylation. These mutations affect two subunits of the DNA-dependent RNA polymerase IV (Pol IV), which is essential for the biogenesis of 24 nt endogenous siRNAs. A mutation in RNA-dependent RNA polymerase 2, another essential gene for the biogenesis of endogenous 24-nt siRNAs, also rescued the defects in miRNA methylation of hen1-2, revealing a previously unsuspected, negative influence of siRNAs on HEN1-mediated miRNA methylation. In addition, our findings imply the existence of a negative modifier of HEN1 activity in the Columbia genetic background.     

Position-dependent effects of locked nucleic acid (LNA) on DNA sequencing and PCR primers

Genomes are becoming heavily annotated with important features. Analysis of these features often employs oligonucleotides that hybridize at defined locations. When the defined location lies in a poor sequence context, traditional design strategies may fail. Locked Nucleic Acid (LNA) can enhance oligonucleotide affinity and specificity. Though LNA has been used in many applications, formal design rules are still being defined. To further this effort we have investigated the effect of LNA on the performance of sequencing and PCR primers in AT-rich regions, where short primers yield poor sequencing reads or PCR yields. LNA was used in three positional patterns: near the 5′ end (LNA-5′), near the 3′ end (LNA-3′) and distributed throughout (LNA-Even). Quantitative measures of sequencing read length (Phred Q30 count) and real-time PCR signal (cycle threshold, C_T) were characterized using two-way ANOVA. LNA-5′ increased the average Phred Q30 score by 60% and it was never observed to decrease performance. LNA-5′ generated cycle thresholds in quantitative PCR that were comparable to high-yielding conventional primers. In contrast, LNA-3′ and LNA-Even did not improve read lengths or C_T. ANOVA demonstrated the statistical significance of these results and identified significant interaction between the positional design rule and primer sequence.     

Shared RNA-binding sites for interacting members of the Drosophila ELAV family of neuronal proteins

The product of the Drosophila embryonic lethal abnormal visual system is a conserved protein (ELAV) necessary for normal neuronal differentiation and maintenance. It possesses three RNA-binding domains and is involved in the regulation of RNA metabolism. The long elav 3′-untranslated region (3′-UTR) is necessary for autoregulation. We used RNA-binding assays and in vitro selection to identify the ELAV best binding site in the elav 3′-UTR. This site resembles ELAV-binding sites identified previously in heterologous targets, both for its nucleotide sequence and its significant affinity for ELAV (K_d 40 nM). This finding supports our model that elav autoregulation depends upon direct interaction between ELAV and elav RNA. We narrowed down the best binding site to a 20 nt long sequence A(U₅)A(U₃)G(U₂)A(U₆) in an alternative 3′ exon. We propose and test a model in which the regulated use of this alternative 3′ exon is involved in normal elav regulation. Found in NEurons (FNE), another neuronal RNA-binding protein paralogous to ELAV, also binds this site. These observations provide a molecular basis for the in vivo interactions reported previously between elav and fne.     

RNase III cleavage demonstrates a long range RNA: RNA duplex element flanking the hepatitis C virus internal ribosome entry site

Here, we show that Escherichia coli Ribonuclease III cleaves specifically the RNA genome of hepatitis C virus (HCV) within the first 570 nt with similar efficiency within two sequences which are ~400 bases apart in the linear HCV map. Demonstrations include determination of the specificity of the cleavage sites at positions C²⁷ and U³³ in the first (5′) motif and G⁴³⁹ in the second (3′) motif, complete competition inhibition of 5′ and 3′ HCV RNA cleavages by added double-stranded RNA in a 1:6 to 1:8 weight ratio, respectively, 50% reverse competition inhibition of the RNase III T7 R1.1 mRNA substrate cleavage by HCV RNA at 1:1 molar ratio, and determination of the 5′ phosphate and 3′ hydroxyl end groups of the newly generated termini after cleavage. By comparing the activity and specificity of the commercial RNase III enzyme, used in this study, with the natural E.coli RNase III enzyme, on the natural bacteriophage T7 R1.1 mRNA substrate, we demonstrated that the HCV cuts fall into the category of specific, secondary RNase III cleavages. This reaction identifies regions of unusual RNA structure, and we further showed that blocking or deletion of one of the two RNase III-sensitive sequence motifs impeded cleavage at the other, providing direct evidence that both sequence motifs, besides being far apart in the linear RNA sequence, occur in a single RNA structural motif, which encloses the HCV internal ribosome entry site in a large RNA loop.     

Sequence specific integration by the family 1 casposase from Candidatus Nitrosopumilus koreensis AR1

Casposase, a homolog of Cas1 integrase, is encoded by a superfamily of mobile genetic elements known as casposons. While family 2 casposase has been well documented in both function and structure, little is known about the other three casposase families. Here, we studied the family 1 casposase lacking the helix-turn-helix (HTH) domain from Candidatus Nitrosopumilus koreensis AR1 (Ca. N. koreensis). The determinants for integration by Ca. N. koreensis casposase were extensively investigated, and it was found that a 13-bp target site duplication (TSD) sequence, a minimal 3-bp leader and three different nucleotides of the TSD sequences are indispensable for target specific integration. Significantly, the casposase can site-specifically integrate a broad range of terminal inverted repeat (TIR)-derived oligonucleotides ranging from 7-nt to ∼4000-bp, and various oligonucleotides lacking the 5′-TTCTA-3′ motif at the 3′ end of TIR sequence can be integrated efficiently. Furthermore, similar to some Cas1 homologs, the casposase utilizes a 5′-ATAA-3′ motif in the TSD as a molecular ruler to dictate nucleophilic attack at 9-bp downstream of the end of the ruler during the spacer-side integration. By characterizing the family 1 Ca. N. koreensis casposase, we have extended our understanding on mechanistic similarities and evolutionary connections between casposons and the adaptation elements of CRISPR-Cas immunity.     

In vitro selection of a 5'-purine ribonucleotide transferase ribozyme

Here we report in vitro selection of a novel ribozyme that catalyzes the 5′-nucleotidyl transfer reaction forming the 2′–5′ phosphodiester bond. This ribozyme was retrieved as a sole sequence in the pool enriched for the 5′-triphosphate-dependent activities in incorporating ATP-γS. The originally selected ribozyme consisting of 109-nucleotide (nt) was miniaturized to 45-nt M4 ribozyme via a series of mutation studies, and based on this mini-ribozyme a trans-acting system was constructed. One of the most challenging tasks in our study was to determine the chemistry occurring at the 5′-ppp site. We utilized various analytical methods including MALDI-TOF analysis of the product generated by the trans-acting system and elucidated the chemistry to be 3′→5′ mononucleotide extension forming the 2′–5′ phosphodiester bond. Interestingly, M4 ribozyme promiscuously accepts a variety of purine nucleotides bearing 5′-mono-, di- and triphosphates as substrates. This remarkable ability of M4 ribozyme would lead us to the development of a new tool for the 5′-modification of RNAs with unique chemical groups.