Unique tertiary and neighbor interactions determine conservation patterns of Cis Watson-Crick A/G base-pairs

X-ray, phylogenetic and quantum chemical analysis of molecular interactions and conservation patterns of cis Watson-Crick (W.C.) A/G base-pairs in 16S rRNA, 23S rRNA and other molecules was carried out. In these base-pairs, the A and G nucleotides interact with their W.C. edges with glycosidic bonds oriented cis relative to each other. The base-pair is stabilised by two hydrogen bonds, the C1'-C1' distance is enlarged and the G(N2) amino group is left unpaired. Quantum chemical calculations show that, in the absence of other interactions, the unpaired amino group is substantially non-planar due to its partial sp(3) pyramidalization, while the whole base-pair is internally propeller twisted and very flexible. The unique molecular properties of the cis W.C. A/G base-pairs make them distinct from other base-pairs. They occur mostly at the ends of canonical helices, where they serve as interfaces between the helix and other motifs. The cis W.C. A/G base-pairs play crucial roles in natural RNA structures with salient sequence conservation patterns. The key contribution to conservation is provided by the unpaired G(N2) amino group that is involved in a wide range of tertiary and neighbor contacts in the crystal structures. Many of them are oriented out of the plane of the guanine base and utilize the partial sp(3) pyramidalization of the G(N2). There is a lack of A/G to G/A covariation, which, except for the G(N2) position, would be entirely isosteric. On the contrary, there is a rather frequent occurrence of G/A to G/U covariation, as the G/U wobble base-pair has an unpaired amino group in the same position as the cis W.C. G/A base-pair. The cis W.C. A/G base-pairs are not conserved when there is no tertiary or neighbor interaction. Obtaining the proper picture of the interactions and phylogenetic patterns of the cis W.C. A/G base-pairs requires a detailed analysis of the relation between the molecular structures and the energetics of interactions at a level of single H-bonds and contacts.     

Structure of the pyrimidine-rich internal loop in the poliovirus 3'-UTR: the importance of maintaining pseudo-2-fold symmetry in RNA helices containing two adjacent non-canonical base-pairs

Formation of non-canonical base-pairs in RNA often plays a very important functional role. In addition they frequently serve as factors in stabilizing the secondary structure elements that provide the frame of large compact RNA structures. Here we describe the structure of an internal loop containing a 5'CU3'/5'UU3' non-canonical tandem base-pair motif, which is conserved within the 3'-UTR of poliovirus-like enteroviruses. Structural details reveal striking regularities of the local helix geometry, resulting from alternating geometrical adjustments, which are important for understanding and predicting stabilities and configurations of tandem non-canonical base-pairs. The C-U and U-U base-pairs severely contract the minor groove of the sugar-phosphate backbone, which might be important for protein recognition or binding to other RNA elements.     

Molecular dynamics simulation reveals conformational switching of water-mediated uracil-cytosine base-pairs in an RNA duplex

A 4 ns molecular dynamics simulation of an RNA duplex (r-GGACUUCGGUCC)(2 )in solution with Na+ and Cl- as counterions was performed. The X-ray structure of this duplex includes two water-mediated uracil-cytosine pairs. In contrast to the other base-pairs in the duplex the water-mediated pairs switch between different conformations. One conformation corresponds to the geometry of the water-mediated UC pairs in the duplex X-ray structure with water acting both as hydrogen-bond donor and acceptor. Another conformation is close to that of a water-mediated UC base-pair found in the X-ray structure of the 23 S rRNA sarcin/ricin domain. In this case the oxygen of the water molecule is linked to two-base donor sites. For a very short time also a direct UC base-pair and a further conformation that is similar to the one found in the RNA duplex structure but exhibits an increased H3(U)...N3(C) distance is observed. Water molecules with unusually long residence times are involved in the water-mediated conformations. These results indicate that the dynamic behaviour of the water-mediated UC base-pairs differs from that of the duplex Watson-Crick and non-canonical guanine-uracil pairs with two or three direct hydrogen bonds. The conformational variability and increased flexibility has to be taken into account when considering these base-pairs as RNA building blocks and as recognition motifs.     

AA.AG@helix.ends: A:A and A:G base-pairs at the ends of 16 S and 23 S rRNA helices.

This study reveals that AA and AG oppositions occur frequently at the ends of helices in RNA crystal and NMR structures in the PDB database and in the 16 S and 23 S rRNA comparative structure models, with the G usually 3' to the helix for the AG oppositions. In addition, these oppositions are frequently base-paired and usually in the sheared conformation, although other conformations are present in NMR and crystal structures. These A:A and A:G base-pairs are present in a variety of structural environments, including GNRA tetraloops, E and E-like loops, interfaced between two helices that are coaxially stacked, tandem G:A base-pairs, U-turns, and adenosine platforms. Finally, given structural studies that reveal conformational rearrangements occurring in regions of the RNA with AA and AG oppositions at the ends of helices, we suggest that these conformationally unique helix extensions might be associated with functionally important structural rearrangements.     

Solution structure of a DNA double helix incorporating four consecutive non-Watson-Crick base-pairs

A series of DNA 21-mers containing a variety of the 4 x 4 internal loop sequence 5'-CAAG-3'/3'-ACGT-5' were studied using nuclear magnetic resonance (NMR) methodology and distance geometry (DG)/molecular dynamics (MD) approaches. Such oligomers exhibit excellent resolution in the NMR spectra and reveal many unusual NOEs (nuclear Overhauser effect) that allow for the detailed characterization of a DNA hairpin incorporating a track of four different non-Watson-Crick base-pairs in the stem. These include a wobble C.A base-pair, a sheared A.C base-pair, a sheared A.G base-pair, and a wobble G.T base-pair. Significantly different twisting angles were observed between the base-pairs in internal loop that results with excellent intra-strand and inter-strand base stacking within the four consecutive mismatches and the surrounding canonical base-pairs. This explains why it melts at 52 degrees C even though five out of ten base-pairs in the stem adopt non-Watson-Crick pairs. However, the 4 x 4 internal loop still fits into a B-DNA double helix very well without significant change in the backbone torsion angles; only zeta torsion angles between the tandem sheared base-pairs are changed to a great extent from the gauche(-) domain to the trans domain to accommodate the cross-strand base stacking in the internal loop. The observation that several consecutive non-canonical base-pairs can stably co-exist with Watson-Crick base-pairs greatly increases the limited repertoire of irregular DNA folds and reveals the possibility for unusual structural formation in the functionally important genomic regions that have potential to become single-stranded.     

The structure of the HIV-1 RRE high affinity rev binding site at 1.6 A resolution

The crystal structure of a 28 nt RNA fragment containing the human immunodeficiency virus type 1 (HIV-1) Rev response element high affinity binding site for Rev protein has been solved at 1.6 A resolution. The overall structure of the RRE helix is greatly distorted from A-form geometry by the presence of two purine-purine base-pairs and two single nucleotide bulges. G48 and G71 form a Hoogsteen-type asymmetric base-pair with G71 adopting a syn conformation. The non-canonical regions in the unliganded Rev response element molecule narrow the major groove width with respect to standard A-RNA. The Rev response element structure observed here represents a closed form of the Rev binding site and differs from conformations of the RNA observed previously by solution NMR studies.     

Identity and geometry of a base triple in 16S rRNA determined by comparative sequence analysis and molecular modeling

Comparative sequence analysis complements experimental methods for the determination of RNA three-dimensional structure. This approach is based on the concept that different sequences within the same gene family form similar higher-order structures. The large number of rRNA sequences with sufficient variation, along with improved covariation algorithms, are providing us with the opportunity to identify new base triples in 16S rRNA. The three-dimensional conformations for one of our strongest candidates involving U121 (C124:G237) and/or U121 (U125:A236) (Escherichia coli sequence and numbering) are analyzed here with different molecular modeling tools. Molecular modeling shows that U121 interacts with C124 in the U121 (C124:G237) base triple. This arrangement maintains isomorphic structures for the three most frequent sequence motifs (approximately 93% of known bacterial and archaeal sequences), is consistent with chemical reactivity of U121 in E. coli ribosomes, and is geometrically favorable. Further, the restricted set of observed canonical (GU, AU, GC) base-pair types at positions 124:237 and 125:236 is consistent with the fact that the canonical base-pair sets (for both base pairs) that are not observed in nature prevent the formation of the 121 (124:237) base triple. The analysis described here serves as a general scheme for the prediction of specific secondary and tertiary structure base pairing where there is a network of correlated base changes.     

Stacking geometry for non‐canonical G:U wobble base pair containing dinucleotide sequences in RNA: dispersion‐corrected DFT‐D study

Emergence of thousands of crystal structures of noncoding RNA molecules indicates its structural and functional diversity. RNA function is based upon a large variety of structural elements which are specifically assembled in the folded molecules. Along with the canonical Watson‐Crick base pairs, different orientations of the bases to form hydrogen‐bonded non‐canonical base pairs have also been observed in the available RNA structures. Frequencies of occurrences of different non‐canonical base pairs in RNA indicate their important role to maintain overall structure and functions of RNA. There are several reports on geometry and energetic stabilities of these non‐canonical base pairs. However, their stacking geometry and stacking stability with the neighboring base pairs are not well studied. Among the different non‐canonical base pairs, the G:U wobble base pair (G:U W:WC) is most frequently observed in the RNA double helices. Using quantum chemical method and available experimental data set we have studied the stacking geometry of G:U W:WC base pair containing dinucleotide sequences in roll‐slide parameters hyperspace for different values of twist. This study indicates that the G:U W:WC base pair can stack well with the canonical base pairs giving rise to large interaction energy. The overall preferred stacking geometry in terms of roll, twist and slide for the eleven possible dinucleotide sequences is seen to be quite dependent on their sequences. © 2015 Wiley Periodicals, Inc. Biopolymers 103: 328–338, 2015.     

Database proton NMR chemical shifts for RNA signal assignment and validation

The Biological Magnetic Resonance Data Bank contains NMR chemical shift depositions for 132 RNAs and RNA-containing complexes. We have analyzed the ¹H NMR chemical shifts reported for non-exchangeable protons of residues that reside within A-form helical regions of these RNAs. The analysis focused on the central base pair within a stretch of three adjacent base pairs (BP triplets), and included both Watson–Crick (WC; G:C, A:U) and G:U wobble pairs. Chemical shift values were included for all 4³ possible WC-BP triplets, as well as 137 additional triplets that contain one or more G:U wobbles. Sequence-dependent chemical shift correlations were identified, including correlations involving terminating base pairs within the triplets and canonical and non-canonical structures adjacent to the BP triplets (i.e. bulges, loops, WC and non-WC BPs), despite the fact that the NMR data were obtained under different conditions of pH, buffer, ionic strength, and temperature. A computer program (RNAShifts) was developed that enables convenient comparison of RNA ¹H NMR assignments with database predictions, which should facilitate future signal assignment/validation efforts and enable rapid identification of non-canonical RNA structures and RNA-ligand/protein interaction sites.     

Genetic conversion of G.C base-pairs to A.U base-pairs in a transfer RNA

The cloverleaf stem segments of the suppressor gene of bacteriophage T4 tRNA(Gln) contain ten G.C and ten A.U base-pairs. To gain a better appreciation of the G.C base-pair requirement, we isolated multiple mutants of this suppressor gene in which base-pairs of G.C were replaced by A.U. One active suppressor gene contained only A.U base-pairs on the anticodon stem, indicating that G.C base-pairs in this region of tRNA(Gln) are not essential for function. In contrast, replacement was not possible at two base-pairs on the D stem and at one base-pair on the T stem.     

Inferring the conformation of RNA base pairs and triples from patterns of sequence variation.

The success of comparative analysis in resolving RNA secondary structure and numerous tertiary interactions relies on the presence of base covariations. Although the majority of base covariations in aligned sequences is associated to Watson-Crick base pairs, many involve non-canonical or restricted base pair exchanges (e.g. only G:C/A:U), reflecting more specific structural constraints. We have developed a computer program that determines potential base pairing conformations for a given set of paired nucleotides in a sequence alignment. This program (ISOPAIR) assumes that the base pair conformation is maintained through sequence variation without significantly affecting the path of the sugar-phosphate backbone. ISOPAIR identifies such 'isomorphic' structures for any set of input base pair or base triple sequences. The program was applied to base pairs and triples with known structures and sequence exchanges. In several instances, isomorphic structures were correctly identified with ISOPAIR. Thus, ISOPAIR is useful when assessing non-canonical base pair conformations in comparative analysis. ISOPAIR applications are limited to those cases where unusual base pair exchanges indeed reflect a non-canonical conformation.     

Stacking geometry between two sheared Watson-Crick basepairs: Computational chemistry and bioinformatics based prediction

BackgroundMolecular modeling of RNA double helices is possible using most probable values of basepair parameters obtained from crystal structure database. The A:A w:wC non-canonical basepair, involving Watson-Crick edges of two Adenines in cis orientation, appears quite frequently in database. Bimodal distribution of its Shear, due to two different H-bonding schemes, introduces the confusion in assigning most the probable value. Its effect is pronounced when the A:A w:wC basepair stacks on Sheared wobble G:U W:WC basepairs.MethodsWe employed molecular dynamics simulations of three possible double helices with GAG, UAG and GAU sequence motifs at their centers and quantum chemical calculation for non-canonical A:A w:wC basepair stacked on G:U W:WC basepair.ResultsWe noticed stable structures of GAG motif with specifically negative Shear of the A:A basepair but stabilities of the other motifs were not found with A:A w:wC basepairing. Hybrid DFT-D and MP2 stacking energy analyses on dinucleotide step sequences, A:A w:wC::G:U W:WC and A:A w:wC::U:G W:WC reveal that viable orientation of A:A::G:U prefers one of the H-bonding modes with negative Shear, supported by crystal structure database. The A:A::U:G dinucleotide, however, prefers structure with only positive Shear.ConclusionsThe quantum chemical calculations explain why MD simulations of GAG sequence motif only appear stable. In the cases of the GAU and UAG motifs “tug of war” situation between positive and negative Shears of A:A w:wC basepair induces conformational plasticity.General significanceWe have projected comprehensive reason behind the promiscuous nature of A:A w:wC basepair which brings occasional structural plasticity.     

Identification and assignment of base pairs in three helical stems of Torulopsis utilis ribosomal 5S RNA and its RNase T1 and RNase T2 cleaved fragments via 500-MHz proton homonuclear overhauser enhancements

Imino proton resonances in the downfield region (10-14 ppm) of the 500-MHz 1H NMR spectrum of Torulopsis utilis 5S RNA are identified (A X U, G X C, or G X U) and assigned to base pairs in helices I, IV, and V via analysis of homonuclear Overhauser enhancements (NOE) from intact T. utilis 5S RNA, its RNase T1 and RNase T2 digested fragments, and a second yeast (Saccharomyces cerevisiae) 5S RNA whose nucleotide sequence differs at only six residues from that of T. utilis 5S RNA. The near-identical chemical shifts and NOE behavior of most of the common peaks from these four RNAs strongly suggest that helices I, IV, and V retain the same conformation after RNase digestion and that both T. utilis and S. cerevisiae 5S RNAs share a common secondary and tertiary structure. Of the four G X U base pairs identified in the intact 5S RNA, two are assigned to the terminal stem (helix I) and the other two to helices IV and V. Seven of the nine base pairs of the terminal stem have been assigned. Our experimental demonstration of a G X U base pair in helix V supports the 5S RNA secondary structural model of Luehrsen and Fox [Luehrsen, K. R., & Fox, G.E. (1981) Proc. Natl. Acad. Sci. U.S.A. 78, 2150-2154]. Finally, the base-pair proton peak assigned to the terminal G X U in helix V of the RNase T2 cleaved fragment is shifted downfield from that in the intact 5S RNA, suggesting that helices I and V may be coaxial in intact T. utilis 5S RNA.     

Interaction of 23S ribosomal RNA helices 89 and 91 of Escherichia coli contributes to the activity of IF2 but is insignificant for elongation factors functioning

The non-canonical base-pair C2475/G2529 joins helices 89 and 91 of the 23S rRNA in the large subunit of E. coli ribosomes. These nucleotides are located at the "crossroads" between the peptidyl transferase center, the sarcin-ricin loop and the GTPase-associated center. We probed the functional role of nucleotides C2475/G2529 by the mutations C2475G, C2475G/G2529C and deltaA2471/U2479 of 23S rRNA. All these mutations had no influence on the elongation factors activity but had different effects on the cell growth, 23S rRNA conformation and translation initiation. C2475G/G2529C and C2475G mutations led to more or less substantial decrease in IF2.GDPNP binding to the ribosomes, and IF2-assisted initiation complex formation. Ribosome-dependent GTPase activity of IF2 was enhanced by both C2475G/G2529C and C2475G mutations. Mutation deltaA2471/U2479 has no influence on IF2.GDPNP binding to the ribosome, but reduces IF2-dependent formation of initiation complex and the ribosome-dependent GTPase activity. Thus, the contact between helices 89 and 91 is important for efficient IF2 functioning in translation initiation.     

New information content in RNA base pairing deduced from quantitative analysis of high-resolution structures

Non-canonical base pairs play important roles in organizing the complex three-dimensional folding of RNA. Here, we outline methodology developed both to analyze the spatial patterns of interacting base pairs in known RNA structures and to reconstruct models from the collective experimental information. We focus attention on the structural context and deformability of the seven pairing patterns found in greatest abundance in the helical segments in a set of well-resolved crystal structures, including (i–ii) the canonical A·U and G·C Watson–Crick base pairs, (iii) the G·U wobble pair, (iv) the sheared G·A pair, (v) the A·U Hoogsteen pair, (vi) the U·U wobble pair, and (vii) the G·A Watson–Crick-like pair. The non-canonical pairs stand out from the canonical associations in terms of apparent deformability, spanning a broader range of conformational states as measured by the six rigid-body parameters used to describe the spatial arrangements of the interacting bases, the root-mean-square deviations of the base-pair atoms, and the fluctuations in hydrogen-bonding geometry. The deformabilties, the modes of base-pair deformation, and the preferred sites of occurrence depend on sequence. We also characterize the positioning and overlap of the base pairs with respect to the base pairs that stack immediately above and below them in double-helical fragments. We incorporate the observed positions of the bases, base pairs, and intervening phosphorus atoms in models to predict the effects of the non-canonical interactions on overall helical structure.     

Solvent-accessible surfaces of nucleic acids

Static solvent-accessible surface areas were calculated for DNA and RNA double helices of varied conformation, composition and sequence, for the single helix of poly(rC), and for a transfer RNA. The results show that for DNA and RNA double helices, two thirds of the water-accessible surface area become buried on double helix formation; phosphate oxygens retain near maximal exposure while the bases are 80% buried. Transfer RNA exposes slightly less surface per residue than does double-helical RNA, despite the presence of several additional “modified” groups, all of which are exposed significantly.When a probe corresponding to a single water molecule is used, both the total and atom type exposures are very similar for A-DNA and B-DNA, although marked differences appear in the major and minor groove exposures between the two conformations. For a given base-pair, the accessible surface area buried upon double-helical stacking is nearly constant (within 5%) for different sequences of neighboring base-pairs.For probes larger than single water molecules, there exist considerable differences in the total and atom type exposures of A-DNA and B-DNA. Conformational transitions between the A-DNA and B-DNA helical forms can thus be related to differences in the accessible areas for “structured” water, or a secondary hydration shell, rather than to interactions with individual water molecules of the primary hydration shell. The base-composition dependence of DNA helical conformation can be explained in terms of the opposing effects of thymine methyl groups of A · T base-pairs and the amino groups of G · C base-pairs upon the solvent within the grooves.The area calculations show that primarily the major groove of B-DNA and the minor groove of A-DNA have sufficient accessible surface area to be recognized by a probe size corresponding to the side-chains of amino acids.     

RNA Sequence Determinants for Aminoglycoside Binding to an A-site rRNA Model Oligonucleotide

The codon-anticodon interaction on the ribosome occurs in the A site of the 30 S subunit. Aminoglycoside antibiotics, which bind to ribosomal RNA in the A site, cause misreading of the genetic code and inhibit translocation. Biochemical studies and nuclear magnetic resonance spectroscopy were used to characterize the interaction between the aminoglycoside antibiotic paromomycin and a small model oligonucle otide that mimics the A site ofEscherichia coli16 S ribosomal RNA. Upon chemical modification, the RNA oligonucleotide exhibits an accessibility pattern similar to that of 16 S rRNA in the 30 S subunit. In addition, the oligonucleotide binds specifically aminoglycoside antibiotics. The anti biotic binding site forms an asymmetric internal loop, caused by non-canonical base-pairs. Nucleotides that are important for binding of paromomycin were identified by performing quantitative footprinting on oligonucleotide sequence variants and include the C1407·G1494 base-pair, and A·U base-pair at positions 1410/1490, and nucleotides A1408, A1493 and U1495. The asymmetry of the internal loop, which requires the presence of a nucleotide in position 1492, is also crucial for antibiotic binding. Introduction into the oligonucleotide of base changes that are known to confer aminoglycoside resistance in 16 S rRNA result in weaker binding of paromomycin to the oligonucleotide. Oligonucleotides homologous to eukaryotic rRNA sequences show reduced binding of paromomycin, suggesting a physical origin for the species-specific action of aminoglycosides.     

Contribution of opening and bending dynamics to specific recognition of DNA damage

Guanine-uracil (G.U) wobble base-pairs are a detrimental lesion in DNA. Previous investigations have shown that such wobble base-pairs are more prone to base-opening than the normal G.C base-pairs. To investigate the sequence-dependence of base-pair opening we have performed 5ns molecular dynamics simulations on G.U wobble base-pairs in two different sequence contexts, TGT/AUA and CGC/GUG. Furthermore, we have investigated the effect of replacing the guanine base in each sequence with a fluorescent guanine analogue, 6-methylisoxanthopterin (6MI). Our results indicate that each sequence opens spontaneously towards the major groove in the course of the simulations. The TGT/AUA sequence has a greater proportion of structures in the open state than the CGC/GUG sequence. Incorporation of 6MI yields wobble base-pairs that open more readily than their guanine counterparts. In order of increasing open population, the sequences are ordered as CGC相似文献