Simulations of RNA base pairs in a nanodroplet reveal solvation-dependent stability.
Ontology highlight
ABSTRACT:
SUBMITTER: Sykes MT
PROVIDER: S-EPMC1920539 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:In an effort to develop unnatural DNA base pairs we examined six pyridine-based nucleotides, d3MPy, d4MPy, d5MPy, d34DMPy, d35DMPy and d45DMPy. Each bears a pyridyl nucleobase scaffold but they are differentiated by methyl substitution, and were designed to vary both inter- and intra-strand packing within duplex DNA. The effects of the unnatural base pairs on duplex stability demonstrate that the pyridine scaffold may be optimized for stable and selective pairing, and identify one self pair, the pair formed between two d34DMPy nucleotides, which is virtually as stable as a dA:dT base pair in the same sequence context. In addition, we found that the incorporation of either the d34DMPy self pair or a single d34DMPy paired opposite a natural dA significantly increases oligonucleotide hybridization fidelity at other positions within the duplex. Hypersensitization of the duplex to mispairing appears to result from global and interdependent solvation effects mediated by the unnatural nucleotide(s) and the mispair. The results have important implications for our efforts to develop unnatural base pairs and suggest that the unnatural nucleotides might be developed as novel biotechnological tools, diagnostics, or therapeutics for applications where hybridization stringency is important.
| S-EPMC2724283 | biostudies-literature
Project description:Nuclear magnetic resonance spectroscopy and proton exchange have been used to characterize two RNA-DNA hybrids from the tR2 intrinsic transcription terminator site of phage lambda. The hybrids have the same base sequence [5'-GGCGCAGGCC(T/U)(T/U)CC-3'/5'-GGAAGGCC(T/U)GCGCC-3'] but differ from each other by an interchange of DNA and RNA strands. The opening of single base pairs in the two hybrids is characterized by measuring the rates of exchange of imino protons with solvent protons as a function of the concentration of a proton acceptor (ammonia base) at 10 degrees C. The free energy change in the opening reaction provides a measure of the stability of the base pair, while the rates of opening and closing define the base pair dynamics. The results demonstrate that, within the same base sequence context, dA-rU base pairs are less stable than dT-rA base pairs. The differences in stability are enhanced when two dA-rU base pairs are located next to each other in the hybrid structure. For the G-C base pairs, the rates of opening and closing and the stability are affected by the base sequence context and by the nature of the sugar moiety attached to the guanine. The dominant feature of the base sequence is the proximity of the dA-rU base pair, which destabilizes the G-C base pair when the guanine is located on the DNA strand. Two G-C base pairs (namely, those in the fourth and 10th positions) exhibit large differences in their opening and closing rates between the two hybrids, while maintaining the same stability. These results provide the first demonstration that, for RNA-DNA hybrid structures with the same base sequence, the opening dynamics and the stability of individual base pairs are strongly influenced by the chemical nature of each strand.
| S-EPMC2758270 | biostudies-literature
Project description:Most of the hairpin, internal and junction loops that appear single-stranded in standard RNA secondary structures form recurrent 3D motifs, where non-Watson-Crick base pairs play a central role. Non-Watson-Crick base pairs also play crucial roles in tertiary contacts in structured RNA molecules. We previously classified RNA base pairs geometrically so as to group together those base pairs that are structurally similar (isosteric) and therefore able to substitute for each other by mutation without disrupting the 3D structure. Here, we introduce a quantitative measure of base pair isostericity, the IsoDiscrepancy Index (IDI), to more accurately determine which base pair substitutions can potentially occur in conserved motifs. We extract and classify base pairs from a reduced-redundancy set of RNA 3D structures from the Protein Data Bank (PDB) and calculate centroids (exemplars) for each base combination and geometric base pair type (family). We use the exemplars and IDI values to update our online Basepair Catalog and the Isostericity Matrices (IM) for each base pair family. From the database of base pairs observed in 3D structures we derive base pair occurrence frequencies for each of the 12 geometric base pair families. In order to improve the statistics from the 3D structures, we also derive base pair occurrence frequencies from rRNA sequence alignments.
| S-EPMC2673412 | biostudies-literature
Project description:Two identical 5'GACG3' tetra-loop motifs with different stem sequences (called H2 and H3) are found in the 5' end region of Moloney Murine Leukemia Virus (MMLV) genomic RNA. They play important roles in RNA dimerization and encapsidation through two identical tetra-loops (5'GACG3') forming a loop-to-loop kissing complex, the smallest RNA kissing complex ever found in nature. We examined the effects of a loop-closing base pair as well as a stem sequence on the conformational stability of the kissing complex. UV melting analysis and gel electrophoresis were performed on eight RNA sequences mimicking the H2 and H3 hairpin tetra-loops with variation in loop-closing base pairs. Our results show that changing the loop-closing base pair from the wildtype (5'A·U3' for H3, 5'U·A3' for H2) to 5'G·C3'/5'C·G3' has significant effect on the stability of the kissing complexes: the substitution to 5'C·G3' significantly decreases both thermal and mechanical stability, while switching to the 5'G·C3' significantly increases the mechanical stability only. The kissing complexes with the wildtype loop-closing base pairs (5'A·U3' for H3 and 5'U·A3' for H2) show different stability when attached to a different stem sequence (H2 stem vs. H3 stem). This suggests that not only the loop-closing base pair itself, but also the stem sequence, affects the conformational stability of the RNA kissing complex.
| S-EPMC5268073 | biostudies-literature
Project description:Non-Watson-Crick base pairs mediate specific interactions responsible for RNA-RNA self-assembly and RNA-protein recognition. An unambiguous and descriptive nomenclature with well-defined and nonoverlapping parameters is needed to communicate concisely structural information about RNA base pairs. The definitions should reflect underlying molecular structures and interactions and, thus, facilitate automated annotation, classification, and comparison of new RNA structures. We propose a classification based on the observation that the planar edge-to-edge, hydrogen-bonding interactions between RNA bases involve one of three distinct edges: the Watson-Crick edge, the Hoogsteen edge, and the Sugar edge (which includes the 2'-OH and which has also been referred to as the Shallow-groove edge). Bases can interact in either of two orientations with respect to the glycosidic bonds, cis or trans relative to the hydrogen bonds. This gives rise to 12 basic geometric types with at least two H bonds connecting the bases. For each geometric type, the relative orientations of the strands can be easily deduced. High-resolution examples of 11 of the 12 geometries are presently available. Bifurcated pairs, in which a single exocyclic carbonyl or amino group of one base directly contacts the edge of a second base, and water-inserted pairs, in which single functional groups on each base interact directly, are intermediate between two of the standard geometries. The nomenclature facilitates the recognition of isosteric relationships among base pairs within each geometry, and thus facilitates the recognition of recurrent three-dimensional motifs from comparison of homologous sequences. Graphical conventions are proposed for displaying non-Watson-Crick interactions on a secondary structure diagram. The utility of the classification in homology modeling of RNA tertiary motifs is illustrated.
| S-EPMC1370104 | biostudies-other
Project description:Many types of RNAs exist in the human transcriptome, yet only the bacterial ribosome has been exploited as a small molecule drug target. Aside from rRNA, other cellular RNAs such as noncoding RNAs have primarily secondary structure and limited tertiary structure. Within these secondary structures of noncanonically paired and unpaired regions, more than 50% are base paired, with most efforts to target these structures focused on looped regions. A void exists in the availability of small molecules capable of targeting RNA base pairs. Using chemoinformatics, an RNA-focused library enriched for nitrogen-containing heterocycles was developed and tested for binding RNA base pairs, leading to the identification of six selective and previously unknown binders. While all binders were derivatives of benzimidazoles, those with expanded aromatic polycycles bound selectively to AU pairs, while those with flexible urea side chains bound selectively to GC pairs. Two of the three selective GC pair binders can distinguish between two different orientations, 5'GG/3'CC and 5'GC/3'CG pairs. Furthermore, all six molecules showed >50-fold selectivity for RNA over DNA. These studies provide foundational knowledge to better exploit RNA as targets for small molecule chemical probes or lead therapeutics by using modules that target RNA base pairs.
| S-EPMC6325646 | biostudies-literature
Project description:Identification and characterization of base-multiplets, which are essentially mediated by base-pairing interactions, can provide insights into the diversity in the structure and dynamics of complex functional RNAs, and thus facilitate hypothesis driven biological research. The necessary nomenclature scheme, an extension of the geometric classification scheme for base-pairs by Leontis and Westhof, is however available only for base-triplets. In the absence of information on topology, this scheme is not applicable to quartets and higher order multiplets. Here we propose a topology-based classification scheme which, in conjunction with a graph-based algorithm, can be used for the automated identification and characterization of higher order base-multiplets in RNA structures. Here, the RNA structure is represented as a graph, where nodes represent nucleotides and edges represent base-pairing connectivity. Sets of connected components (of n nodes) within these graphs constitute subgraphs representing multiplets of "n" nucleotides. The different topological variants of the RNA multiplets thus correspond to different nonisomorphic forms of these subgraphs. To annotate RNA base-multiplets unambiguously, we propose a set of topology-based nomenclature rules for quartets, which are extendable to higher multiplets. We also demonstrate the utility of our approach toward the identification and annotation of higher order RNA multiplets, by investigating the occurrence contexts of selected examples in order to gain insights regarding their probable functional roles.
| S-EPMC6467009 | biostudies-literature
Project description:Inosine is found naturally in the anticodon loop of tRNA, is a product of adenosine deaminases that act on RNA, and can be used in oligonucleotide probes or to investigate the role of the exocyclic amino group of guanosine. Although the thermodynamics of I·U pairs in RNA have been systematically studied [Wright, D. J., Rice, J. L., Yanker, D. M., and Znosko, B. M. (2007) Biochemistry 46, 4625-4634], the thermodynamics of I·C pairs in RNA have not. Here, we have performed optical melting experiments on a series of RNA duplexes containing I·C pairs and compared their thermodynamics to the same duplexes containing A·C and G-C pairs. Nearest neighbor parameters for single I·C pairs adjacent to Watson-Crick pairs were derived. The derived nearest neighbor parameters are compared to those previously predicted blindly through a reweighting of energy-function collection with conformational ensemble sampling in Rosetta [Chou, F.-C., Kladwang, W., Kappel, K., and Das, R. (2016) Proc. Natl. Acad. Sci. U.S.A. 113, 8430-8435]. Scientists can use these nearest neighbor parameters to calculate the stability of ADAR products and to calculate the stability of an RNA duplex in which G-to-I substitution was used to determine the role of the exocyclic amino group of G.
| S-EPMC6294561 | biostudies-literature
Project description:Prediction of RNA tertiary structure from sequence is an important problem, but generating accurate structure models for even short sequences remains difficult. Predictions of RNA tertiary structure tend to be least accurate in loop regions, where non-canonical pairs are important for determining the details of structure. Non-canonical pairs can be predicted using a knowledge-based model of structure that scores nucleotide cyclic motifs, or NCMs. In this work, a partition function algorithm is introduced that allows the estimation of base pairing probabilities for both canonical and non-canonical interactions. Pairs that are predicted to be probable are more likely to be found in the true structure than pairs of lower probability. Pair probability estimates can be further improved by predicting the structure conserved across multiple homologous sequences using the TurboFold algorithm. These pairing probabilities, used in concert with prior knowledge of the canonical secondary structure, allow accurate inference of non-canonical pairs, an important step towards accurate prediction of the full tertiary structure. Software to predict non-canonical base pairs and pairing probabilities is now provided as part of the RNAstructure software package.
| S-EPMC5690697 | biostudies-literature
Project description:Adenosine to inosine (A?I) editing is the most common modification of double-stranded RNA (dsRNA). This change is mediated by adenosine deaminases acting on RNA (ADARs) enzymes with a preference of U>A>C>G for 5′ neighbor and G>C=A>U or G>C>U=A for 3′ neighbor. A?I editing occurs most frequently in the non-coding regions containing repetitive elements such as ALUs. It leads to disruption of RNA duplex structure, which prevents induction of innate immune response. We employed standard and biased molecular dynamics (MD) simulations to analyze the behavior of RNA duplexes with single and tandem inosine?uracil (I?U) base pairs in different sequence context. Our analysis showed that the I?U pairs induce changes in base pair and base pair step parameters and have different dynamics when compared with standard canonical base pairs. In particular, the first I?U pair from tandem I?U/I?U systems exhibited increased dynamics depending on its neighboring 5′ base. We discovered that UII sequence, which is frequently edited, has lower flexibility compared with other sequences (AII, GII, CII), hence it only modestly disrupts dsRNA. This might indicate that the UAA motifs in ALUs do not have to be sufficiently effective in preventing immune signaling.
| S-EPMC6070904 | biostudies-literature