Stacking properties of a highly hydrophobic dinucleotide sequence, N6, N6-dimethyladenylyl(3' leads to 5')N6, N6-dimethyladenosine, occurring in 16--18-S ribosomal RNA

Predicting coaxial helical stacking in RNA junctions

RNA junctions are important structural elements that form when three or more helices come together in space in the tertiary structures of RNA molecules. Determining their structural configuration is important for predicting RNA 3D structure. We introduce a computational method to predict, at the secondary structure level, the coaxial helical stacking arrangement in junctions, as well as classify the junction topology. Our approach uses a data mining approach known as random forests, which relies on a set of decision trees trained using length, sequence and other variables specified for any given junction. The resulting protocol predicts coaxial stacking within three- and four-way junctions with an accuracy of 81% and 77%, respectively; the accuracy increases to 83% and 87%, respectively, when knowledge from the junction family type is included. Coaxial stacking predictions for the five to ten-way junctions are less accurate (60%) due to sparse data available for training. Additionally, our application predicts the junction family with an accuracy of 85% for three-way junctions and 74% for four-way junctions. Comparisons with other methods, as well applications to unsolved RNAs, are also presented. The web server Junction-Explorer to predict junction topologies is freely available at: http://bioinformatics.njit.edu/junction.

Structural principles from large RNAs

Since the year 2000 a number of large RNA three-dimensional structures have been determined by X-ray crystallography. Structures composed of more than 100 nucleotide residues include the signal recognition particle RNA, group I intron, the GlmS ribozyme, RNAseP RNA, and ribosomal RNAs from Haloarcula morismortui, Escherichia coli, Thermus thermophilus, and Deinococcus radiodurans. These large RNAs are constructed from the same secondary and tertiary structural motifs identified in smaller RNAs but appear to have a larger organizational architecture. They are dominated by long continuous interhelical base stacking, tend to segregate into domains, and are planar in overall shape as opposed to their globular protein counterparts. These findings have consequences in RNA folding, intermolecular interaction, and packing, in addition to studies of design and engineering and structure prediction.

RNA-like polymer model: exact calculation on the Bethe lattice

We consider a lattice polymer model (random walk), in which the walk is allowed to visit lattice bonds at most twice. Such a model might have some relevance to describe statistical properties of RNA molecules. In order to mimic base pairing, we assign an attractive energy term to each doubly visited bond, and a further contribution to each pair of consecutive doubly visited bonds. The latter term is expected to mimic the stacking effect, whereas no effect of sequence, that is, of chemical specificity, is taken into account. The phase diagram is worked out exactly on a Bethe lattice, in a grand-canonical formulation. In the single molecule limit, the system undergoes two different phase transitions upon decreasing temperature: a Theta-like collapse from a swollen "coil" state to a "molten" state, with a low fraction of doubly visited bonds, and subsequently to a "paired" state, with empty or doubly visited bonds only. The stacking effect drives the latter transition from second to first order.

Quantitative analysis of nucleic acid three-dimensional structures

A new computer program to annotate DNA and RNA three-dimensional structures, MC-Annotate, is introduced. The goals of annotation are to efficiently extract and manipulate structural information, to simplify further structural analyses and searches, and to objectively represent structural knowledge. The input of MC-Annotate is a PDB formatted DNA or RNA three-dimensional structure. The output of MC-Annotate is composed of a structural graph that contains the annotations, and a series of HTML documents, one for each nucleotide conformation and base-base interaction present in the input structure. The atomic coordinates of all nucleotides and the homogeneous transformation matrices of all base-base interactions are stored in the structural graph. Symbolic classifications of nucleotide conformations, using sugar puckering modes and nitrogen base orientations around the glycosyl bond, and base-base interactions, using stacking and hydrogen bonding information, are introduced. Peculiarity factors of nucleotide conformations and base-base interactions are defined to indicate their marginalities with all other examples. The peculiarity factors allow us to identify irregular regions and possible stereochemical errors in 3-D structures without interactive visualization. The annotations attached to each nucleotide conformation include its class, its torsion angles, a distribution of the root-mean-square deviations with examples of the same class, the list of examples of the same class, and its peculiarity value. The annotations attached to each base-base interaction include its class, a distribution of distances with examples of the same class, the list of examples of the same class, and its peculiarity value. The distance between two homogeneous transformation matrices is evaluated using a new metric that distinguishes between the rotation and the translation of a transformation matrix in the context of nitrogen bases. MC-Annotate was used to build databases of nucleotide conformations and base-base interactions. It was applied to the ribosomal RNA fragment that binds to protein L11, which annotations revealed peculiar nucleotide conformations and base-base interactions in the regions where the RNA contacts the protein. The question of whether the current database of RNA three-dimensional structures is complete is addressed.

N 2-methylguanosine is iso-energetic with guanosine in RNA duplexes and GNRA tetraloops

Modified nucleotides are resource-intensive alternatives to the four nucleotides that constitute the bulk of natural RNAs. Yet, even in cases where modifications are highly conserved, their functions are difficult to identify. One possible function might be to modulate the stability of RNA structures. To investigate this possibility for N 2-methylguanosine (m2G), which is present in a wide variety of RNAs, we have determined the thermodynamic consequences of substituting m2G for G in G-C Watson-Crick pairs and G@U wobble pairs within RNA duplexes. The m2G substitution is iso-energetic with G in all cases, except for aninternal m2G@U pair, where it has a modest (0.3 kcal/mol) stabilizing effect. We have also examined theconsequences of replacing G by m2G, and A by N 6, N 6-dimethyladenosine (m26A) in the helix 45 tetraloop of 16S rRNA, which would otherwise be a standard GNRA tetraloop. This loop is a conserved, hypermethylated region of the ribosome where methylation appears to modulate activity. m26A substitution destabilizes the tetraloop, presumably because it prevents the formation of the G@A sheared pair it would otherwise contain. m2G substitution has no effect on tetraloop stability. Together, these results suggest that m2G is equally stable as either the s-cis or s-trans rotamer. The lack of a significant effect on secondary structural stability in these systems suggests that m2G is introduced into naturally occurring RNAs for reasons other than modulation of duplex stability.

The structure of a methylated tetraloop in 16S ribosomal RNA

BACKGROUND

Ribosomal RNAs contain many modified nucleotides. The functions of these nucleotides are poorly understood and few of them are strongly conserved. The final stem loop in 16S-like rRNAs is an exception in both regards. In both prokaryotes and eukaryotes, the tetranucleotide loop that caps the 3'-terminal stem contains two N6, N6-dimethyladenosine residues. The sequence and pattern of methylation are conserved within the loop, and there is evidence that these methylated nucleotides play an important role in subunit association and the initiation of protein synthesis. Because of the integral role that helix 45 plays in ribosome function, it is important to know what consequences these methylated nucleotides have on its structure.

RESULTS

We have solved the solution structure of a 14-nucleotide analog of the terminal stem loop of bacterial 16S rRNA, which contains N2-methylguanosine as well as two N6,N6-dimethyladenosines.

CONCLUSIONS

The methylation of the 16S rRNA stem loop completely alters its conformation, which would otherwise be a GNRA tetraloop. It is likely that the conformation of this loop is crucial for its function, having implications for its interaction with ribosomal subunits and its role in the initiation of protein synthesis.

The conformation of a conserved stem-loop structure in ribosomal RNA

The RNA of small ribosomal subunits contains a conserved stem-loop structure near the 3' end. Characteristics for the hairpins are: (a) a nine-basepairs stem: (b) a conserved A-UU-G junction in the stem: (c) a conserved sequence Gm6(2)AM6(2)A sequence in the loop (except yeast mitochondria and mutants from bacteria). We are using UV-optics, micro-calorimetry and 500 MHz-NMR to investigate fragments of about 50 nucleotides cleaved from the 3' ends of small ribosomal subunit RNA's by bacteriocins. Our preliminary conclusions are: (1) Dimethylation of the adenines in the loop destabilizes the hairpin because of an increased stacking; (2) melting of the hairpin starts at the ends as well as in the middle at the A-UU-G junction; (3) basepair substitutions have an unexpectedly large effect on thermal stability.

Intramolecular stacking association and conformation properties of a 'cap' structure, m7G5'pppUm, and the related model compounds

The stacking equilibrium quotient of the m7G5'pppUm unit, which occurs as the 5'-terminal "cap" of certain eukaryotic mRNA's, was determined by temperature-dependent difference spectrophotometry as Kstack = 1.82 at 25 degrees and pH 5. In order to evaluate the contribution of different structural modifications to the net stabilization of the cap structures of mRNA, a variety of compounds related to m7G5'pppUm were synthesized and their stacking properties were studied by the same method and compared. The results are summarized as: (1) Introduction of a methyl group into N-7 of G residue results in an increase in base stacking. (2) Methylation at 2'-OH of U residue also stabilizes the stacked structure of G-containing dimers, but it does not influence stacking interaction in m7G-containing dimers. (3) The effect of different types of internucleotide linkages on the order of stacking tendencies is: N5'ppN' greater than N5'pppN' greater than NpN'. UV hypochromicity and CD spectral measurements of the relevant dimers were also conducted, and the hypochromicity values and CD spectra of dimers in their stacked conformation were estimated by making use of the determined Kstack values. The results indicate that, while 2'-O-methylation exerts very little effect on the stacked conformation of the dimers, methylation at N-7 and the nature of the internucleotide linkage strongly influence the stacked conformation, thereby forming unusual left-handed conformations in m7G5'pppU(m), m7G5'ppU(m), and G5'ppU(m).

Determination of microscopic ionization constants and microscopic stacking equilibrium quotients of adenylyl(3' → 5')adenosine

Using the basic ionization constants for a pair of isomers, m1ApA and Apm1A, and the measured values for the overlapping pK values of ApA, the microscopic ionization constants and intramolecular stacking quotients for the monoprotonated ApA were estimated. The results indicate that, in contrast to the case of GpG, ApA did not exhibit preferential protonation on either site of 3'- and 5'-linked nucleoside bases and neither enhanced nor diminished stacking was observed for ApA and ApA as compared to ApA.

A carbon-13 nuclear magnetic resonance study of the 3'-terminus of 16S ribosomal RNA of Escherichia coli specifically labeled with carbon-13 in the methylgroups of the m6(2)Am6(2)A sequence

30S ribosomes were isolated from a kasugamycin resistant mutant of E. coli that lacks methylgroups on two adjacent adenines in 16S ribosomal RNA. These ribosomes were methylated in vitro with a purified methylating enzyme and 5-S-adenosyl-(13C-methyl)-L-methionine chloride ((13C-methyl)-SAM) as methyldonor. After in situ cleavage of the 16S ribosomal RNA by the bacteriocin cloacin DF13, the 49 nucleotide fragment from the 3'-end of the RNA was isolated. The carbon-13 nuclear magnetic resonance spectra of the fragment at various temperatures were compared with those of 6-N-dimethyladenosine (m6(2)A) and 6-N-dimethyladenylyl-(3' leads to 5')-6-N-dimethyladenosine (m6(2)Am6(2)A). The data show that the two methylated adenines, which are part of a four membered hairpin loop, show a strong tendency to be stacked in analogy to the dinucleotide m6(2)Am6(2).

Influence of the 2'-hydroxyl group and of 6-N-methylation on the conformation of adenine dinucleoside monophosphates in solution. A nuclear magnetic resonance and circular dichroism study

Proton NMR studies at 360 MHz are reported on the adenine dinucleoside monophosphates N6-dimethyladenyly(3'-5')-N6-dimethyladenosine (m(6)(2)Apm(6)(2)A), ApA, rApdA, dAprA and on the methyl phosphate esters of the monomers m(6)(2)Ap, pm(6)(2)A, Ap and pA. Complete 1H-NMR spectral assignments are given. The dimers were also investigated by means of circular dichroism to obtain accurate thermodynamic parameters of the stacking equilibrium. With the aid of the thermodynamic data NMR coupling constants are extrapolated to values appropriate to the stacked conformers. A modernized version of pseudorotation analysis is used to delineate the conformational behaviour of the ribose and 2'-deoxyribose rings. It is shown that the unmethylated dimers can be arranged in two groups (dApdA/dAprA vs ApA/rApdA) according to their melting temperatures. ApA and the fully N6-methylated dimer m(6)(2)Apm(6)(2)A prefer to adopt the classical right-handed N-N stacked conformation. Both dimers with a 2'-deoxyribose ring at the 5'-OH end (dApdA and dAprA) behave similarly and occur in solution as a 75:25 mixture of S-S and S-N stacked states. The fully stacked hybrid dimer rApdA displays an unexpectedly large amount of S conformers (greater than 40%) in both sugar rings. This finding is rationalized by the postulation of a right-handed helical S-S stacked state on the basis of NMR and circular dichroic data.

Destabilization of secondary structure in 16S ribosomal RNA by dimethylation of two adjacent adenosines

Fragments comprising the 49 nucleotides from the 3'-end have been purified from 16S ribosomal RNA of wild-type Escherichia coli and from a kasugamycin-resistant mutant that specifically lacks dimethylation of two adjacent adenines near the 3'-terminus. These fragments, obtained after treatment of ribosomes in vitro with the bacteriocin cloacin DF13, were used to study the effect of the methyl groups on the temperature dependent unfolding of double-stranded regions. Both fragments contain at least 3 independent melting transitions, of which the one with the highest Tm corresponds with the unfolding of a nine-basepair long central hairpin. Dimethylation of the adenines in the loop of this hairpin lowers the melting temperature (Tm) by approximately 2 degrees C at 0.2 M NaCl and by about 5 degrees C at 0.15 M NaCl. It is suggested that m6(2)Am6(2)A is more antagonistic to loop formation that ApA and that the function of the methyl groups is to help to destabilize the 3'-terminal hairpin in 16S rRNA in order to facilitate intermolecular interactions.

Adenosine dimethylation of 16S ribosomal RNA: effect of the methylgroups on local conformational stability as deduced from electrophoretic mobility of RNA fragments in denaturing polyacrylamide gels

The electrophoretic mobility of RNA fragments derived from the 3'-end of 16S rRNA on slabs of polyacrylamide gel in the presence of urea is strongly influenced by dimethylation of the N6-aminogroup of two adjacent adenosines. This is not due to the presence of the methylgroups per se, but must be ascribed to an effect of methylation on long range intramolecular interactions at these denaturing conditions. When it is assumed that the electrophoretic mobilities of the RNA fragments in the polyacrylamide matrix are determined by the conformational state(s) of the fragments, dimethylation of the adenosines leads in the smaller fragments to a less compact average conformation and in the larger fragments to a more compact average conformation. An effort is made to comprehend the effects of adenosine dimethylation in terms of secondary structure based on nucleotide sequence.