The hepatitis delta virus is a single-stranded circular RNA virus, which is characterized by high self-complementarity. About 70% of the genome sequences can form base-pairs with internal nucleotides. There are many studies on the evolution of the hepatitis delta virus. However, the secondary structure has not been taken into account in these studies. In this study, we developed a method to examine the effect of base pairing as a constraint on the nucleotide substitutions during the evolution of the hepatitis delta virus. The method revealed that the base pairing can reduce the evolutionary rate in the non-coding region of the virus. In addition, it is suggested that the non-coding nucleotides without base pairing may be under some constraint, and that the intensity of the constraint is weaker than that by the base pairing but stronger than that on the synonymous site.

A-to-I editing catalyzed by adenosine deaminase acting on RNAs impacts numerous physiological and biochemical processes that are essential for cellular functions and is a big contributor to the infectivity of certain RNA viruses. The outcome of this deamination leads to changes in the eukaryotic transcriptome functionally resembling A-G transitions since inosine preferentially pairs with cytosine. Moreover, hyper-editing or multiple A to G transitions in clusters were detected in measles virus. Inosine modifications either directly on viral RNA or on cellular RNA can have antiviral or pro-viral repercussions. While many of the significant roles of inosine in cellular RNAs are well understood, the effects of hyper-editing of A to I on viral polymerase activity during RNA replication remain elusive. Moreover, biological strategies such as molecular cloning and RNA-seq for transcriptomic interrogation rely on RT-polymerase chain reaction with little to no emphasis placed on the first step, reverse transcription, which may reshape the sequencing results when hypermodification is present. In this study, we systematically explore the influence of inosine modification, varying the number and position of inosines, on decoding outcomes using three different reverse transcriptases (RTs) followed by standard Sanger sequencing. We find that inosine alone or in clusters can differentially affect the RT activity. To gain structural insights into the accommodation of inosine in the polymerase site of HIV-1 reverse transcriptase (HIV-1-RT) and how this structural context affects the base pairing rules for inosine, we performed molecular dynamics simulations of the HIV-1-RT. The simulations highlight the importance of the protein-nucleotide interaction as a critical factor in deciphering the base pairing behavior of inosine clusters. This effort sets the groundwork for decrypting the physiological significance of inosine and linking the fidelity of reverse transcriptase and the possible diverse transcription outcomes of cellular RNAs and/or viral RNAs where hyper-edited inosines are present in the transcripts.

Natural RNAs utilize extensive chemical modifications to diversify their structures and functions. 2-Thiouridine geranylation is a special hydrophobic tRNA modification that has been discovered very recently in several bacteria, such as Escherichia coli, Enterobacter aerogenes, Pseudomonas aeruginosa and Salmonella Typhimurium The geranylated residues are located in the first anticodon position of tRNAs specific for lysine, glutamine and glutamic acid. This big hydrophobic terpene functional group affects the codon recognition patterns and reduces frameshifting errors during translation. We aimed to systematically study the structure, function and biosynthesis mechanism of this geranylation pathway, as well as answer the question of why nature uses such a hydrophobic modification in hydrophilic RNA systems. Recently, we have synthesized the deoxy-analog of S-geranyluridine and showed the geranylated T-G pair is much stronger than the geranylated T-A pair and other mismatched pairs in the B-form DNA duplex context, which is consistent with the observation that the geranylated tRNA(Glu) UUC recognizes GAG more efficiently than GAA. In this manuscript we report the synthesis and base pairing specificity studies of geranylated RNA oligos. We also report extensive molecular simulation studies to explore the structural features of the geranyl group in the context of A-form RNA and its effect on codon-anticodon interaction during ribosome binding.

Hoogsteen (HG) base pairing is characterized by a 180° rotation of the purine base with respect to the Watson-Crick-Franklin (WCF) motif. Recently, it has been found that both conformations coexist in a dynamical equilibrium and that several biological functions require HG pairs. This relevance has motivated experimental and computational investigations of the base-pairing transition. However, a systematic simulation of sequence variations has remained out of reach. Here, we employ advanced path-based methods to perform unprecedented free-energy calculations. Our methodology enables us to study the different mechanisms of purine rotation, either remaining inside or after flipping outside of the double helix. We study seven different sequences, which are neighbor variations of a well-studied A⋅T pair in A6-DNA. We observe the known effect of A⋅T steps favoring HG stability, and find evidence of triple-hydrogen-bonded neighbors hindering the inside transition. More importantly, we identify a dominant factor: the direction of the A rotation, with the 6-ring pointing either towards the longer or shorter segment of the chain, respectively relating to a lower or higher barrier. This highlights the role of DNA's relative flexibility as a modulator of the WCF/HG dynamic equilibrium. Additionally, we provide a robust methodology for future HG proclivity studies.

There are >170 naturally occurring RNA chemical modifications, with both known and unknown biological functions. Analytical methods for detecting chemical modifications and for analyzing their effects are relatively limited and have had difficulty keeping pace with the demand for RNA chemical biology and biochemistry research. Some modifications can affect the ability of RNA to hybridize with its complementary sequence or change the selectivity of base pairing. Here, we investigate the use of affinity-based DNA nanoswitches to resolve energetic differences in hybridization. We found that a single m3C modification can sufficiently destabilize hybridization to abolish a detection signal, while an s4U modification can selectively hybridize with G over A. These results establish proof of concept for using DNA nanoswitches to detect certain RNA modifications and analyzing their effects in base pairing stability and specificity.

The reverse transcriptase telomerase adds telomeric repeats to chromosome ends to counteract telomere shortening and thereby assures genomic stability in dividing human cells. Key parameters in telomere homeostasis are the frequency with which telomerase engages the chromosome end and the number of telomeric repeats it adds during each association event. To study telomere elongation in vivo, we have established a live-cell imaging assay to track individual telomerase ribonucleoproteins in CRISPR-edited HeLa cells. Using this assay and the drug imetelstat, which is a competitive inhibitor of telomeric DNA binding, we demonstrate that stable association of telomerase with the single-stranded overhang of the chromosome end requires telomerase-DNA base pairing. Furthermore, we show that telomerase processivity contributes to telomere elongation in vivo. Together, these findings provide new insight into the dynamics of telomerase recruitment and the importance of processivity in maintaining telomere length in human cancer cells.

The ribonuclease Dicer produces microRNAs (miRNAs) and small interfering RNAs that are handed over to Ago proteins to control gene expression by targeting complementary sequences within transcripts. Interestingly, a growing number of reports have demonstrated that the activity of Dicer may extend beyond the biogenesis of small regulatory RNAs. Among them, a report from our latest studies revealed that human Dicer facilitates base pairing of complementary sequences present in two nucleic acids, thus acting as a nucleic acid annealer. Accordingly, in this manuscript, we address how RNA structure influences the annealing activity of human Dicer. We show that Dicer supports hybridization between a small RNA and a complementary sequence of a longer RNA in vitro, even when both complementary sequences are trapped within secondary structures. Moreover, we show that under applied conditions, human Ago2, a core component of RNA-induced silencing complex, displays very limited annealing activity. Based on the available data from new-generation sequencing experiments regarding the RNA pool bound to Dicer in vivo, we show that multiple Dicer-binding sites within mRNAs also contain miRNA targets. Subsequently, we demonstrate in vitro that Dicer but not Ago2 can anneal miRNA to its target present within mRNA. We hypothesize that not all miRNA duplexes are handed over to Ago proteins. Instead, miRNA-Dicer complexes could target specific sequences within transcripts and either compete or cooperate for binding sites with miRNA-Ago complexes. Thus, not only Ago but also Dicer might be directly involved in the posttranscriptional control of gene expression.

Bone growth is linked to expansion of nearby organs, as is the case for the cranial base and the brain. Here, we focused on development of the mouse clivus, a sloping surface of the basioccipital bone, to define mechanisms underlying morphological changes in bone in response to brain enlargement. Histological analysis indicated that both endocranial and ectocranial cortical bone layers in the basioccipital carry the osteoclast surface dorsally and the osteoblast surface ventrally. Finite element analysis of mechanical stress on the clivus revealed that compressive and tensile stresses appeared mainly on respective dorsal and ventral surfaces of the basioccipital bone. Osteoclastic bone resorption occurred primarily in the compression area, whereas areas of bone formation largely coincided with the tension area. These data collectively suggest that compressive and tensile stresses govern respective localization of osteoclasts and osteoblasts. Developmental analysis of the basioccipital bone revealed the clivus to be angled in early postnatal wild-type mice, whereas its slope was less prominent in Tnfsf11-/- mice, which lack osteoclasts. We propose that osteoclast-osteoblast "trans-pairing" across cortical bone is primarily induced by mechanical stress from growing organs and regulates shape and size of bones that encase the brain.

The chemical synthesis of ribonucleic acids (RNA) with novel chemical modifications is largely driven by the motivation to identify eligible functional probes for the various applications in life sciences. To this end, we have a strong focus on the development of novel fluorinated RNA derivatives that are powerful in NMR spectroscopic analysis of RNA folding and RNA ligand interactions. Here, we report on the synthesis of 2'-SCF3 pyrimidine nucleoside containing oligoribonucleotides and the comprehensive investigation of their structure and base pairing properties. While this modification has a modest impact on thermodynamic stability when it resides in single-stranded regions, it was found to be destabilizing to a surprisingly high extent when located in double helical regions. Our NMR spectroscopic investigations on short single-stranded RNA revealed a strong preference for C2'-endo conformation of the 2'-SCF3 ribose unit. Together with a recent computational study (L. Li, J. W. Szostak, J. Am. Chem. Soc. 2014, 136, 2858-2865) that estimated the extent of destabilization caused by a single C2'-endo nucleotide within a native RNA duplex to amount to 6 kcal mol(-1) because of disruption of the planar base pair structure, these findings support the notion that the intrinsic preference for C2'-endo conformation of 2'-SCF3 nucleosides is most likely responsible for the pronounced destabilization of double helices. Importantly, we were able to crystallize 2'-SCF3 modified RNAs and solved their X-ray structures at atomic resolution. Interestingly, the 2'-SCF3 containing nucleosides that were engaged in distinct mismatch arrangements, but also in a standard Watson-Crick base pair, adopted the same C3'-endo ribose conformations as observed in the structure of the unmodified RNA. Likely, strong crystal packing interactions account for this observation. In all structures, the fluorine atoms made surprisingly close contacts to the oxygen atoms of the corresponding pyrimidine nucleobase (O2), and the 2'-SCF3 moieties participated in defined water-bridged hydrogen-bonding networks in the minor groove. All these features allow a rationalization of the structural determinants of the 2'-SCF3 nucleoside modification and correlate them to base pairing properties.

Pairs of nucleotides within functional nucleic acid secondary structures often display evidence of coevolution that is consistent with the maintenance of base-pairing. Here, we introduce a sequence evolution model, MESSI (Modeling the Evolution of Secondary Structure Interactions), that infers coevolution associated with base-paired sites in DNA or RNA sequence alignments. MESSI can estimate coevolution while accounting for an unknown secondary structure. MESSI can also use graphics processing unit parallelism to increase computational speed. We used MESSI to infer coevolution associated with GC, AU (AT in DNA), GU (GT in DNA) pairs in noncoding RNA alignments, and in single-stranded RNA and DNA virus alignments. Estimates of GU pair coevolution were found to be higher at base-paired sites in single-stranded RNA viruses and noncoding RNAs than estimates of GT pair coevolution in single-stranded DNA viruses. A potential biophysical explanation is that GT pairs do not stabilize DNA secondary structures to the same extent that GU pairs do in RNA. Additionally, MESSI estimates the degrees of coevolution at individual base-paired sites in an alignment. These estimates were computed for a SHAPE-MaP-determined HIV-1 NL4-3 RNA secondary structure. We found that estimates of coevolution were more strongly correlated with experimentally determined SHAPE-MaP pairing scores than three nonevolutionary measures of base-pairing covariation. To assist researchers in prioritizing substructures with potential functionality, MESSI automatically ranks substructures by degrees of coevolution at base-paired sites within them. Such a ranking was created for an HIV-1 subtype B alignment, revealing an excess of top-ranking substructures that have been previously identified as having structure-related functional importance, among several uncharacterized top-ranking substructures.

HSUR2 is a viral non-coding RNA (ncRNA) that functions as a microRNA (miRNA) adaptor. HSUR2 inhibits apoptosis in infected cells by recruiting host miRNAs miR-142-3p and miR-16 to mRNAs encoding apoptotic factors. HSUR2's target recognition mechanism is not understood. It is also unknown why HSUR2 utilizes miR-16 to downregulate only a subset of transcripts. We developed a general method for individual-nucleotide resolution RNA-RNA interaction identification by crosslinking and capture (iRICC) to identify sequences mediating interactions between HSUR2 and target mRNAs in vivo. Mutational analyses confirmed identified HSUR2-mRNA interactions and validated iRICC as a method that confidently determines sequences mediating RNA-RNA interactions in vivo. We show that HSUR2 does not display a 'seed' region to base-pair with most target mRNAs, but instead uses different regions to interact with different transcripts. We further demonstrate that this versatile mode of interaction via variable base-pairing provides HSUR2 with a mechanism for differential miRNA recruitment.

Single strands of EcoRI-cleaved pCR1 deoxyribonucleic acid were examined by electron microscopy for intrastrand base pairing by using partial denaturing conditions. The locations of three stem and loop structures were mapped relative to the inverted repeat of Tn903. Potential roles and origins of these loops are discussed.

Uveal melanoma (UM), the most common primary intraocular tumor in adults, has been extensively characterized by omics technologies during the last 5 yr. Despite the discovery of gene signatures, the molecular actors driving cancer aggressiveness are not fully understood, and UM is still associated with very poor overall survival (OS) at the metastatic stage. By defining the miR-16 interactome, we revealed that miR-16 mainly interacts via non-canonical base-pairing to a subset of RNAs, promoting their expression levels. Consequently, the canonical miR-16 activity, involved in the RNA decay of oncogenes, such as cyclin D3, is impaired. This non-canonical base-pairing can explain both the derepression of miR-16 targets and the promotion of oncogene expression observed in patients with poor OS in two cohorts. miR-16 activity, assessment using our RNA signature, discriminates the patient's OS as effectively as current methods. To the best of our knowledge, this is the first time that a predictive signature has been composed of genes belonging to the same mechanism (miR-16) in UM. Altogether, our results strongly suggest that UM is a miR-16 disease.

Accurate pre-mRNA splicing is crucial for gene expression. The 5' splice site (5' ss)--the highly diverse element at the 5' end of introns--is initially recognized via base-pairing to the 5' end of the U1 small nuclear RNA (snRNA). However, many natural 5' ss have a poor match to the consensus sequence, and are predicted to be weak. Using genetic suppression experiments in human cells, we demonstrate that some atypical 5' ss are actually efficiently recognized by U1, in an alternative base-pairing register that is shifted by one nucleotide. These atypical 5' ss are phylogenetically widespread, and many of them are conserved. Moreover, shifted base-pairing provides an explanation for the effect of a 5' ss mutation associated with pontocerebellar hypoplasia. The unexpected flexibility in 5' ss-U1 base-pairing challenges an established paradigm and has broad implications for splice-site prediction algorithms and gene-annotation efforts in genome projects.

The synthesis of a caged RNA phosphoramidite building block containing the oxidatively damaged base 5-hydroxycytidine (5-HOrC) has been accomplished. To determine the effect of this highly mutagenic lesion on complementary base recognition and coding properties, this building block was incorporated into a 12-mer oligoribonucleotide for T(m) and CD measurements and a 31-mer template strand for primer extension experiments with HIV-, AMV- and MMLV-reverse transcriptase (RT). In UV-melting experiments, we find an unusual biphasic transition with two distinct T(m)'s when 5-HOrC is paired against a DNA or RNA complement with the base guanine in opposing position. The higher T(m) closely matches that of a C-G base pair while the lower is close to that of a C-A mismatch. In single nucleotide extension reactions, we find substantial misincorporation of dAMP and to a lesser extent dTMP, with dAMP almost equaling that of the parent dGMP in the case of HIV-RT. A working hypothesis for the biphasic melting transition does not invoke tautomeric variability of 5-HOrC but rather local structural perturbations of the base pair at low temperature induced by interactions of the 5-HO group with the phosphate backbone. The properties of this RNA damage is discussed in the context of its putative biological function.

During normal translation, the binding of a release factor to one of the three stop codons (UGA, UAA or UAG) results in the termination of protein synthesis. However, modification of the initial uridine to a pseudouridine (Ψ) allows efficient recognition and read-through of these stop codons by a transfer RNA (tRNA), although it requires the formation of two normally forbidden purine-purine base pairs. Here we determined the crystal structure at 3.1 Å resolution of the 30S ribosomal subunit in complex with the anticodon stem loop of tRNA(Ser) bound to the ΨAG stop codon in the A site. The ΨA base pair at the first position is accompanied by the formation of purine-purine base pairs at the second and third positions of the codon, which show an unusual Watson-Crick/Hoogsteen geometry. The structure shows a previously unsuspected ability of the ribosomal decoding centre to accommodate non-canonical base pairs.

Understanding the dynamic interactions of ligands to DNA is important in DNA-based nanotechnologies. By structurally tracking the dissociation of Hoechst 33258-bound DNA (d(CGCAAATTTGCG)2) complex (H-DNA) with T-jump 2D-IR spectroscopy, the ligand is found to strongly disturb the stability of the three C:G base pairs adjacent to A:T the binding site, with the broken base pairs being more than triple at 100 ns. The strong stabilization effect of the ligand on DNA duplex makes this observation quite striking, which dramatically increases the melting temperature and dissociation time. MD simulations demonstrate an important role of hydration water and counter cations in maintaining the separation of terminal base pairs. The hydrogen bonds between the ligand and thymine carbonyls are crucial in stabilizing H-DNA, whose breaking signal appearing prior to the complete dissociation. Thermodynamic analysis informs us that H-DNA association is a concerted process, where H cooperates with DNA single strands in forming H-DNA.

Genetic information is encoded in the DNA double helix, which, in its physiological milieu, is characterized by the iconical Watson-Crick nucleo-base pairing. Recent NMR relaxation experiments revealed the transient presence of an alternative, Hoogsteen (HG) base pairing pattern in naked DNA duplexes, and estimated its relative stability and lifetime. In contrast with DNA, such structures were not observed in RNA duplexes. Understanding HG base pairing is important because the underlying "breathing" motion between the two conformations can significantly modulate protein binding. However, a detailed mechanistic insight into the transition pathways and kinetics is still missing. We performed enhanced sampling simulation (with combined metadynamics and adaptive force-bias method) and Markov state modeling to obtain accurate free energy, kinetics, and the intermediates in the transition pathway between Watson-Crick and HG base pairs for both naked B-DNA and A-RNA duplexes. The Markov state model constructed from our unbiased MD simulation data revealed previously unknown complex extrahelical intermediates in the seemingly simple process of base flipping in B-DNA. Extending our calculation to A-RNA, for which HG base pairing is not observed experimentally, resulted in relatively unstable, single-hydrogen-bonded, distorted Hoogsteen-like bases. Unlike B-DNA, the transition pathway primarily involved base paired and intrahelical intermediates with transition timescales much longer than that of B-DNA. The seemingly obvious flip-over reaction coordinate (i.e., the glycosidic torsion angle) is unable to resolve the intermediates. Instead, a multidimensional picture involving backbone dihedral angles and distance between hydrogen bond donor and acceptor atoms is required to gain insight into the molecular mechanism.

Programmability of DNA sequences enables the formation of synthetic DNA nanostructures and their macromolecular assemblies such as DNA hydrogels. The base pair-level interaction of DNA is a foundational and powerful mechanism to build DNA structures at the nanoscale; however, its temperature sensitivity and weak interaction force remain a barrier for the facile and scalable assembly of DNA structures toward higher-order structures. We conducted this study to provide an alternative, non-base-pairing approach to connect nanoscale DNA units to yield micrometer-sized gels based on the sequential phase transition of amphiphilic unit structures. Strong electrostatic interactions between DNA nanostructures and polyelectrolyte spermines led to the formation of giant phase-separated aggregates of monomer units. Gelation could be initiated by the addition of NaCl, which weakened the electrostatic DNA-spermine interaction while attractive interactions between cholesterols created stable networks by crosslinking DNA monomers. In contrast to the conventional DNA gelation techniques, our system used solid aggregates as a precursor for DNA microgels. Therefore, in situ gelation could be achieved by depositing aggregates on the desired substrate and subsequently initiating a phase transition. Our approach can expand the utility and functionality of DNA hydrogels by using more complex nucleic acid assemblies as unit structures and combining the technique with top-down microfabrication methods.

8-Nitro-2'-deoxyguanosine (8-nitrodG) is a relatively unstable, mutagenic lesion of DNA that is increasingly believed to be associated with tissue inflammation. Due to the lability of the glycosidic bond, 8-nitrodG cannot be incorporated into oligodeoxynucleotides (ODNs) by chemical DNA synthesis and thus very little is known about its physicochemical properties and base-pairing preferences. Here we describe the synthesis of 8-nitro-2'-O-methylguanosine, a ribonucleoside analogue of this lesion, which is sufficiently stable to be incorporated into ODNs. Physicochemical studies demonstrated that 8-nitro-2'-O-methylguanosine adopts a syn conformation about the glycosidic bond; thermal melting studies and molecular modelling suggest a relatively stable syn-8-nitroG·anti-G base pair. Interestingly, when this lesion analogue was placed in a primer-template system, extension of the primer by either avian myeloblastosis virus reverse transcriptase (AMV-RT) or human DNA polymerase β (pol β), was significantly impaired, but where incorporation opposite 8-nitroguanine did occur, pol β showed a 2:1 preference to insert dA over dC, while AMV-RT incorporated predominantly dC. The fact that no 8-nitroG·G base pairing is seen in the primer extension products suggests that the polymerases may discriminate against this pairing system on the basis of its poor geometric match to a Watson-Crick pair.