Theoretical analysis of noncanonical base pairing interactions in RNA molecules

Noncanonical base pairs in RNA have strong structural and functional implications but are currently not considered for secondary structure predictions. We present results of comparative ab initio studies of stabilities and interaction energies for the three standard and 24 selected unusual RNA base pairs reported in the literature. Hydrogen added models of isolated base pairs, with heavy atoms frozen in their ‘away from equilibrium’ geometries, built from coordinates extracted from NDB, were geometry optimized using HF/6-31G** basis set, both before and after unfreezing the heavy atoms. Interaction energies, including BSSE and deformation energy corrections, were calculated, compared with respective single point MP2 energies, and correlated with occurrence frequencies and with types and geometries of hydrogen bonding interactions. Systems having two or more N-H...O/N hydrogen bonds had reasonable interaction energies which correlated well with respective occurrence frequencies and highlighted the possibility of some of them playing important roles in improved secondary structure prediction methods. Several of the remaining base pairs with one N-H...O/N and/or one C-H...O/N interactions respectively, had poor interaction energies and negligible occurrences. High geometry variations on optimization of some of these were suggestive of their conformational switch like characteristics.     

Theoretical analysis of noncanonical base pairing interactions in RNA molecules

Noncanonical base pairs in RNA have strong structural and functional implications but are currently not considered for secondary structure predictions. We present results of comparative ab initio studies of stabilities and interaction energies for the three standard and 24 selected unusual RNA base pairs reported in the literature. Hydrogen added models of isolated base pairs, with heavy atoms frozen in their 'away from equilibrium' geometries, built from coordinates extracted from NDB, were geometry optimized using HF/6-31G** basis set, both before and after unfreezing the heavy atoms. Interaction energies, including BSSE and deformation energy corrections, were calculated, compared with respective single point MP2 energies, and correlated with occurrence frequencies and with types and geometries of hydrogen bonding interactions. Systems having two or more N-H...O/N hydrogen bonds had reasonable interaction energies which correlated well with respective occurrence frequencies and highlighted the possibility of some of them playing important roles in improved secondary structure prediction methods. Several of the remaining base pairs with one N-H...O/N and/or one C-H...O/N interactions respectively, had poor interaction energies and negligible occurrences. High geometry variations on optimization of some of these were suggestive of their conformational switch like characteristics.     

Quantum-chemical analysis of C-H...O and C-H...N interactions in RNA base pairs--H-bond versus anti-H-bond pattern

Geometries and interaction energies of unusual UU and AA base pairs with one standard hydrogen bond (H-bond) and additional C-H...O or C-H...N contacts have been determined by quantum-chemical methods taking into account electron correlation. Whereas the C-H bond length in the UU C-H...O contact increases upon complex formation (H-bond pattern), the C-H bond of the AA C-H....N interaction is shortened (anti-H-bond pattern). The same properties are found for model complexes between U or A and formaldehyde that have intermolecular C-H...acceptor contacts but no standard H-bonds. Both the C-H...acceptor H-bond and anti-H-bond interactions are attractive. A possible influence of the donor CH group charge distribution on the interaction pattern is discussed.     

Variable base pairing in a helix of eubacterial 5 S ribosomal RNA points to the existence of a conformational switch

A survey of 160 published sequences of eubacterial 5 S rRNAs shows that there exists structural variability in one of the helices of the generally accepted secondary structure model. Four structural variants are found, which differ with respect to the position and the number of bulges present. Most eubacterial 5 S RNAs fit into at least two of these conformations. A reaction scheme connecting the four observed conformations by changes in the base pairing scheme is proposed. Each of the known 5 S RNA sequences fits into conformations interconvertible by the proposed reactions.     

STRAL: progressive alignment of non-coding RNA using base pairing probability vectors in quadratic time

MOTIVATION: Alignment of RNA has a wide range of applications, for example in phylogeny inference, consensus structure prediction and homology searches. Yet aligning structural or non-coding RNAs (ncRNAs) correctly is notoriously difficult as these RNA sequences may evolve by compensatory mutations, which maintain base pairing but destroy sequence homology. Ideally, alignment programs would take RNA structure into account. The Sankoff algorithm for the simultaneous solution of RNA structure prediction and RNA sequence alignment was proposed 20 years ago but suffers from its exponential complexity. A number of programs implement lightweight versions of the Sankoff algorithm by restricting its application to a limited type of structure and/or only pairwise alignment. Thus, despite recent advances, the proper alignment of multiple structural RNA sequences remains a problem. RESULTS: Here we present StrAl, a heuristic method for alignment of ncRNA that reduces sequence-structure alignment to a two-dimensional problem similar to standard multiple sequence alignment. The scoring function takes into account sequence similarity as well as up- and downstream pairing probability. To test the robustness of the algorithm and the performance of the program, we scored alignments produced by StrAl against a large set of published reference alignments. The quality of alignments predicted by StrAl is far better than that obtained by standard sequence alignment programs, especially when sequence homologies drop below approximately 65%; nevertheless StrAl's runtime is comparable to that of ClustalW.     

Aromatic interactions in tryptophan-containing peptides: crystal structures of model tryptophan peptides and phenylalanine analogs.

The crystal structures of the peptides, Boc-Leu-Trp-Val-OMe (1), Ac-Leu-Trp-Val-OMe (2a and 2b), Boc-Leu-Phe-Val-OMe (3), Ac-Leu-Phe-Val-OMe (4), and Boc-Ala-Aib-Leu-Trp-Val-OMe (5) have been determined by X-ray diffraction in order to explore the nature of interactions between aromatic rings, specifically the indole side chain of Trp residues. Peptide 1 adopts a type I beta-turn conformation stabilized by an intramolecular 4-->1 hydrogen bond. Molecules of 1 pack into helical columns stabilized by two intermolecular hydrogen bonds, Leu(1)NH...O(2)Trp(2) and IndoleNH...O(1)Leu(1). The superhelical columns further pack into the tetragonal space group P4(3) by means of a continuous network of indole-indole interactions. Peptide 2 crystallizes in two polymorphic forms, P2(1) (2a) and P2(1)2(1)2(1) (2b). In both forms, the peptide backbone is extended, with antiparallel beta-sheet association being observed in crystals. Extended strand conformations and antiparallel beta-sheet formation are also observed in the Phe-containing analogs, Boc-Leu-Phe-Val-OMe (3) and Ac-Leu-Phe-Val-OMe (4). Peptide 5 forms a short stretch of 3(10)-helix. Analysis of aromatic-aromatic and aromatic-amide interactions in the structures of peptides, 1, 2a, 2b are reported along with the examples of 14 Trp-containing peptides from the Cambridge Crystallographic Database. The results suggest that there is no dramatic preference for a preferred orientation of two proximal indole rings. In Trp-containing peptides specific orientations of the indole ring, with respect to the preceding and succeeding peptide units, appear to be preferred in beta-turns and extended structures.     

B cell development leads off with a base hit: dU:dG mismatches in class switching and hypermutation

The mechanisms underlying somatic hypermutation (SHM) and class switch recombination (CSR) have been the subject of much debate. Recent studies from the Neuberger and Honjo labs have lent insight into these distinct processes, and we discuss a new, comprehensive model for how AID, uracil DNA glycosylase (UNG) and the mismatch repair system function in both SHM and CSR.     

Structure of Musashi1 in a complex with target RNA: the role of aromatic stacking interactions

Mammalian Musashi1 (Msi1) is an RNA-binding protein that regulates the translation of target mRNAs, and participates in the maintenance of cell 'stemness' and tumorigenesis. Msi1 reportedly binds to the 3'-untranslated region of mRNA of Numb, which encodes Notch inhibitor, and impedes initiation of its translation by competing with eIF4G for PABP binding, resulting in triggering of Notch signaling. Here, the mechanism by which Msi1 recognizes the target RNA sequence using its Ribonucleoprotein (RNP)-type RNA-binding domains (RBDs), RBD1 and RBD2 has been revealed on identification of the minimal binding RNA for each RBD and determination of the three-dimensional structure of the RBD1:RNA complex. Unique interactions were found for the recognition of the target sequence by Msi1 RBD1: adenine is sandwiched by two phenylalanines and guanine is stacked on the tryptophan in the loop between β1 and α1. The minimal recognition sequences that we have defined for Msi1 RBD1 and RBD2 have actually been found in many Msi1 target mRNAs reported to date. The present study provides molecular clues for understanding the biology involving Musashi family proteins.     

A novel guanine-guanine base pairing: crystal structure of a complex between 7-methylguanosine and its iodide.

7-Methylguanosine, one of the biologically important minor nucleosides, could be crystallized as a complex of its zwitterionic form and its iodide, and the crystal structure was determined by the X-ray diffraction method. The crystals belong to the triclinic space group P1 with the unit cell dimensions: a = 7.678(1), b = 18.094(3), and c = 5.711(1) A, alpha = 79.32(1), beta = 80.14(1) and gamma = 76.90(1) degrees. The structure was solved by the heavy atom method and refined by the least-squares method to give a final R index of 0.075. The novel reverse Watson-Crick type base pairing observed between a positively charged molecule and a deprotonated one indicates that the deprotonation at the N(1) position promoted by the alkylation at the N(7) position may interrupt the formation of the normal Watson-Crick type GC base pair. The conformations about the glycosidic bond and the sugar puckering are quite different between the two molecules: the former has anti and C(4')-exo,C(3')-endo and the latter syn and C(1')-exo-C(2')-endo.     

Computational evaluation of intermolecular interactions of a universal base 3-nitropyrrole in stacked dimers and DNA duplexes

The stacking interactions between a universal base of 3-nitropyrrole (3NP) and four canonical nucleobases were studied by means of ab initio molecular orbital calculations. The stabilities of the complexes are comparable to those of the stacked dimers of canonical bases reported previously. The detailed analysis of the interaction energies revealed the importance of the dipole-dipole interaction included in the Hartree-Fock terms to determine the geometry dependence of the stacking energies. It was also clarified that the dispersion energies included in the electron-correlation terms were essential to obtain adequate stabilities. The contribution of the nitro group was evaluated by the comparative studies of pyrrole and 3NP. The increased molecular dipole moment and surface are expected to account for the enhancement of the stability of the stacked dimers containing 3NP. The force field parameters required for calculation of the molecular mechanics of 3NP were obtained for 3NP on the basis of these molecular orbital calculations. The energy-minimized structures obtained by the molecular mechanics calculations of 3NP accorded with those obtained by the molecular orbital calculations described above. A DNA duplex structure containing 3NP-A, 3NP-T, or 3NP-C was calculated by use of these force field parameters. In the case of 3NP-A, the computationally calculated structure was in good agreement with that previously determined by use of (1)H-NMR except for the orientation of the nitro group.     

Altered structural fluctuations in duplex RNA versus DNA: a conformational switch involving base pair opening

DNA and RNA are known to have different structural properties. In the present study, molecular dynamics (MD) simulations on a series of RNA and DNA duplexes indicate differential structural flexibility for the two classes of oligonucleotides. In duplex RNA, multiple base pairs experienced local opening events into the major groove on the nanosecond time scale, while such events were not observed in the DNA simulations. Three factors are indicated to be responsible for the base opening events in RNA: solvent-base interactions, 2'OH(n)-O4'(n+1) intra-strand hydrogen bonding, and enhanced rigid body motion of RNA at the nucleoside level. Water molecules in the major groove of RNA contribute to initiation of base pair opening. Stabilization of the base pair open state is due to a 'conformational switch' comprised of 2'OH(n)-O4'(n+1) hydrogen bonding and a rigid body motion of the nucleoside moiety in RNA. This rigid body motion is associated with decreased flexibility of the glycosyl linkage and sugar moieties in A-form structures. The observed opening rates in RNA are consistent with the imino proton exchange experiments for AU base pairs, although not for GC base pairs, while structural and flexibility changes associated with the proposed conformational switch are consistent with survey data of RNA and DNA crystal structures. The possible relevance of base pair opening events in RNA to its many biological functions is discussed.     

Pyrimidine non-nucleoside analogs: A direct synthesis of a novel class of N-substituted amino and N-sulfonamide derivatives of pyrimidines

A convenient method for the regioselective synthesis of pyrimidine non-nucleoside analogs was developed. This study reports a novel and efficient method for the synthesis of a new type of N-substituted amino methylsulfanylpyrimidines and the corresponding pyrazolo[3,4-d]pyrimidines. This series of compounds was designed through the reaction of dimethyl N-cyanodithioiminocarbonate with 2-cyano-N′-(thiophen-2-yl-, furan-2-yl- and pyridin-4-ylmethylene)acetohydrazide and N′-(2-cyanoacetyl)arylsulfonohydrazides. The scope and limitation of the method are demonstrated. The antibacterial and antifungal activities of the synthesized compounds were also evaluated.     

Crystal structure of E. coli YhbY: a representative of a novel class of RNA binding proteins

E. coli YhbY belongs to a conserved family of hypothetical proteins represented in eubacteria, archaea, and plants (Pfam code UPF0044). Three maize proteins harboring UPF0044-like domains are required for chloroplast group II intron splicing, and bioinformatic data suggest a role for prokaryotic UPF0044 members in translation. The crystal structure of YhbY has been determined. YhbY has a fold similar to that of the C-terminal domain of translation initiation factor 3 (IF3C), which binds to 16S rRNA in the 30S ribosome. Modeling studies indicate that the same surface is highly basic in all members of UPF0044, suggesting a conserved RNA binding surface. Taken together, the evidence suggests that members of UPF0044 constitute a previously unrecognized class of RNA binding domain.     

Identification and characterization of O-acetylpeptidoglycan esterase: a novel enzyme discovered in Neisseria gonorrhoeae

Modification of the bacterial cell wall heteropolymer peptidoglycan by addition of an acetyl group to the C-6 hydroxyl group of N-acetylmuramoyl residues is known to inhibit the activity of muramidases (lysozymes) of innate immune systems. The O-acetylation of peptidoglycan also precludes the action of intrinsic lytic transglycosylases, enzymes that require a free C-6 hydroxyl group to generate their 1,6-anhydromuropeptide products. This class of autolysins is ubiquitous in peptidoglycan-synthesizing bacteria as they are responsible for insertion of pores and flagella, spore formation, and the general metabolism of peptidoglycan. We recently discovered a cluster of genes in the Neisseria gonorrhoeae chromosome that are proposed to participate in peptidoglycan O-acetylation (Weadge, J. T., Pfeffer, J. M., and Clarke, A. J. (2005) BMC Microb. 5, 49). In the current study, we demonstrate that one of these genes, ape1 functions as an O-acetylpeptidoglycan esterase. The ape1 gene was cloned and overexpressed in Escherichia coli as a fusion protein with a hexa-histidine tag. The expressed protein was purified to apparent homogeneity and assayed for activity as an esterase using three different assays involving high-performance liquid chromatography and chromogenic detection methods which measured the release of ester-linked acetate from a variety of polymer and soluble substrates. These assays demonstrated that Ape1 has a higher specific activity on O-acetylated peptidoglycan compared to O-acetylated xylan. Consequently, Ape1 represents the first enzyme characterized as an O-acetylpeptidoglycan esterase. The physicochemical and kinetic parameters of Ape1 were determined using soluble chromogenic substrates for convenience. Thus, its pH optima for stability and activity were observed to be 6.0 and 6.2, respectively, while its optimum temperature for activity was 55 degrees C. Two forms of truncated Ape1 are generated in E. coli, one lacked the complete predicted N-terminal signal sequence, while the second involved a proteolytic cleavage within this signal sequence. The smaller truncated form was localized predominantly to the periplasm, whereas the larger form was mainly associated with the outer membrane, and to a lesser extent, the cytoplasmic membrane, sites expected for the maintenance of peptidoglycan.     

Legacy effects of low N diets on a vernal folivore: deprivation leads to gluttony

We tested for legacy effects of low-N diets offered to newly emerged second-instar spruce budworm (Choristoneura fumiferana) larvae for a duration of either one or two full instars on their growth and nutritional physiology in the sixth instar. The experimental design evaluated the effects of initial diet, final diet, and sex on energy consumption, assimilation, retention, and growth rates. Legacy effects were apparent after two instars of low-N diet exposure and were manifested by elevated ( approximately 10%) consumption rates (RCRs) coupled to elevated ( approximately 10%) growth rates (RGRs) and elevated ( approximately 3%) body energy densities, i.e., heightened fat deposition. However, initial dietary N levels had no legacy main effects on food assimilation efficiencies (ADs), and gross (ECI) and net (ECD) food conversion efficiencies. RCR and AD were dependent on an initial x final diet interaction (i.e., nonlinear legacy effects). RGR depended on an initial diet x sex interaction but not on an initial x final diet interaction. Therefore, the legacy effects of low-N initial diets on RGR and body energy density were simply additive to final diet effects. Final diet universally affected all indices and interacted with sex. Low-N final diets increased RCR ( approximately 41%) and decreased AD (14-18%) but unexpectedly increased ECD (21-24%) and RGR ( approximately 36%). Females generally had higher performance than males on the low-N diets but often only matched males on the high-N diets. Low-N initial diets extended larval development times ( approximately 7-26%) and lowered growth rates (6-24%) to the sixth instar, depending on duration of diet exposure, but did not affect total growth achieved by the start of the sixth instar.     

Comparison of the base pairing properties of a series of nitroazole nucleobase analogs in the oligodeoxyribonucleotide sequence 5'-d(CGCXAATTYGCG)-3'.

The nucleoside analogs 1-(2'-deoxy-beta-D-ribofuranosyl)- 3-nitropyrrole (9), 1-(2'-deoxy-beta-D-ribofuranosyl)-4-nitropyrazole (10), 1-(2'-deoxy-beta-D-ribofuranosyl)-4-nitroimidazole (11) and 1-(2'-deoxy-beta-D-ribofuranosyl)-5-nitroindole (21) were incorporated into the oligonucleotide 5'-d(CGCXAATTYGCG)-3'in the fourth position from the 5'-end. Procedures for synthesis of two of the nitroazole nucleosides, 10 and 11, were developed for this study. Each of the nitroazoles was converted into a 3'-phosphoramidite for oligonucleotide synthesis by conventional automated protocols. Four oligonucleotides were synthesized for each modified nucleoside in order to obtain duplexes in which each of the four natural bases was placed opposite (position 9) the nitroazole. In order to assess the role of the nitro group on base stacking interaction, sequences were also synthesized in which the fourth base was 1-(2'-deoxy-beta-D-ribofuranosyl)pyrazole. Corresponding sequences containing an abasic site, as well as sequences containing inosine, were synthesized for comparison. Thermal melting studies yielded T m values and thermodynamic parameters. Each nucleoside analog displayed a unique pattern of base pairing preferences. The least discriminating analog was 3-nitropyrrole, for which T m values differed by 5 degrees C and Delta G 25 degrees C ranged from -6.1 to -6.5 kcal/mol. 5-Nitroindole gave duplexes with significantly higher thermal stability, with Tm values varying from 35.0 to 46.5 degrees C and -Delta G 25 degrees C ranging from 7.7 to 8.5 kcal/mol. Deoxyinosine (22), a natural analog which has found extensive use as a universal nucleoside, is far less non-discriminating than any of the nitroazole derivatives. Tm values ranged from 35.4 degrees C when paired with G to 62.3 degrees C when paired with C. The significance of the nitro substituent was determined by comparison of the base pairing properties of a simple azole nucleoside, 1-(2'-deoxy-beta-D-ribofuranosyl)pyrazole (12). The pyrazole-containing sequences melt at 10-20 degrees C lower than the corresponding nitropyrazole-containing sequences. On average, the pyrazole-containing sequences were equivalent in stability (average Delta G = -4.8 kcal/mol) to the sequences containing an abasic site (average Delta G = -4.7 kcal/mol).     

NCIR: a database of non-canonical interactions in known RNA structures

The secondary and tertiary structure of an RNA molecule typically includes a number of non-canonical base–base interactions. The known occurrences of these interactions are tabulated in the NCIR database, which can be accessed from http://prion.bchs.uh.edu/bp_type/. The number of examples is now over 1400, which is an increase of >700% since the database was first published. This dramatic increase reflects the addition of data from the recently published crystal structures of the 50S (2.4 Å) and 30S (3.0 Å) ribosomal subunits. In addition, non-canonical interactions observed in published crystal and NMR structures of tRNAs, group I introns, ribozymes, RNA aptamers and synthetic oligonucleotides are included. Properties associated with these interactions, such as sequence context, sugar pucker conformation, glycosidic angle conformation, melting temperature, chemical shift and free energy, are also reported when available. Out of the 29 anticipated pairs with at least two hydrogen bonds, 28 have been observed to date. In addition, several novel examples, not generally predicted, have also been encountered, bringing the total of such pairs to 36. Added to this list are a variety of single, bifurcated, triple and quadruple interactions. The most common non-canonical pairs are the sheared GA, GA imino, AU reverse Hoogsteen, and the GU and AC wobble pairs. The most frequent triple interaction connects N3 of an A with the amino of a G that is also involved in a standard Watson–Crick pair.