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

Diversity of 5S rRNA genes within individual prokaryotic genomes

We examined intragenomic variation of paralogous 5S rRNA genes to evaluate the concept of ribosomal constraints. In a dataset containing 1161 genomes from 779 unique species, 96 species exhibited > 3% diversity. Twenty-seven species with > 10% diversity contained a total of 421 mismatches between all pairs of the most dissimilar copies of 5S rRNA genes. The large majority (401 of 421) of the diversified positions were conserved at the secondary structure level. The high diversity was associated with partial rRNA operon, split operon, or spacer length-related divergence. In total, these findings indicated that there are tight ribosomal constraints on paralogous 5S rRNA genes in a genome despite of the high degree of diversity at the primary structure level.

Prediction of functional tertiary interactions and intermolecular interfaces from primary sequence data

Given the availability of sequence information for many species, one can examine how the sequence of a gene varies among different organisms. This is accomplished by aligning the sequences and observing patterns of conservation, mutation and counter-mutation at different positions in the gene. Imbedded in these patterns is information on energetic coupling and macromolecular interactions, which can be deciphered by application of statistical algorithms. Here we report a robust approach for predicting interactions within (or between) any type of biopolymer, including proteins, RNAs and RNA-protein complexes. Rather than maximize the number of predictions, this approach is designed to detect a limited number of highly significant interactions, thereby providing accurate results from alignments that contain a modest number of sequences (20-60). The versatility and accuracy of the algorithm is demonstrated by the successful prediction of important intramolecular interactions within RNAs, modified RNAs, and proteins, as well as the prediction of RNA-protein and protein-protein interactions.

Evidence for a base triple in the free HIV-1 TAR RNA

We propose the existence of a novel base triple in the HIV-1 TAR hairpin. This triple is supported by covariation of loop residue 31 with residue 22, which is part of an unusual base pair with U40 below the 3-nucleotide bulge. A set of mutants was constructed to test the involvement of bases A22, U31, and U40 in a triple interaction. RNA structure probing, trans-activation assays, and structure modeling are consistent with the existence of this base triple in a bent conformation of the free TAR element. However, disruption of the base triple does not affect binding of a Tat-derived peptide. We therefore compared the structure of free and Tat-bound TAR RNA by footprinting and site-specific cross-linking analyses. These studies indicate that the Tat arginine-rich motif, in addition to its known binding site at the bulge, is in close contact with U31 in the TAR loop. Because binding of Tat to TAR is known to coincide with the formation of a base triple with residues U23, A27, and U38, we hypothesize that Tat binding and the associated straightening of TAR triggers the disruption of the (A22-U40)U31 triple.

Construction and analysis of base-paired regions of the 16S rRNA in the 30S ribosomal subunit determined by constraint satisfaction molecular modelling

Structure models for each of the secondary structure regions from the Escherichia coli 16S rRNA (58 separate elements) were constructed using a constraint satisfaction modelling program to determine which helices deviated from classic A-form geometry. Constraints for each rRNA element included the comparative secondary structure, H-bonding conformations predicted from patterns of base-pair covariation, tertiary interactions predicted from covariation analysis, chemical probing data, rRNA-rRNA crosslinking information, and coordinates from solved structures. Models for each element were built using the MC-SYM modelling algorithm and subsequently were subjected to energy minimization to correct unfavorable geometry. Approximately two-thirds of the structures that result from the input data are very similar to A-form geometry. In the remaining instances, the presence of internal loops and bulges, some sequences (and sequence covariants) and accessory information require deviation from A-form geometry. The structures of regions containing more complex base-pairing arrangements including the central pseudoknot, the 530 region, and the pseudoknot involving base-pairing between G570-U571/A865-C866 and G861-C862/G867-C868 were predicted by this approach. These molecular models provide insight into the connection between patterns of H-bonding, the presence of unpaired nucleotides, and the overall geometry of each element.

Novel design architecture for genetic stability of recombinant poliovirus: the manipulation of G/C contents and their distribution patterns increases the genetic stability of inserts in a poliovirus-based RPS-Vax vector system

Poliovirus has been studied as a live recombinant vaccine vector because of its attractive characteristics. The genetic instability, however, has hampered recombinant polioviruses (PVs) from being developed as an appropriate vaccine. A variety of different foreign inserts were cloned directly into our poliovirus Sabin 1-based RPS-Vax vector system, resulting in the production of recombinant PVs. The genetic stability of each recombinant PV was examined during 12 rounds of consecutive passage. It was found that the genetic stability of the recombinants was not well correlated with their insert size. Instead, elevated stability was frequently observed in recombinants with inserts of high G/C contents. Furthermore, a comparative study using different constructs of the human immunodeficiency virus env gene revealed that the internal deletion of the unstable insert was seemingly caused by the presence of the adjacent A/T-rich region. The instability of these inserts was completely remedied by (i) increasing the G/C contents and (ii) replacing the local A/T-rich region with the G/C-rich codon without a change of the amino acid. This means that stability is closely associated with the G/C content and the G/C distribution pattern. To see whether these findings can be applied to the design of genetically stable recombinant PV, we have reconstructed the heteromultimeric insert based on our design architecture, including the above-mentioned G/C rules and the template/ligation-free PCR protocol. The heteromultimeric insert was very unstable, as expected, but the manipulated insert with the same amino acid sequence showed complete genetic stability, not only in vitro, but also in vivo. Even though this guideline was established with our RPS-Vax vector system, to some extent, it can also be applied to other live viral vaccine vectors.

Predicting U-turns in ribosomal RNA with comparative sequence analysis

The U-turn is a well-known RNA motif characterized by a sharp reversal of the RNA backbone following a single-stranded uridine base. In experimentally determined U-turn motifs, the nucleotides 3' to the turn are frequently involved in tertiary interactions, rendering this motif particularly attractive in RNA modeling and functional studies. The U-turn signature is composed of an UNR sequence pattern flanked by a Y:Y, Y:A (Y=pyrimidine) or G:A base juxtaposition. We have identified 33 potential UNR-type U-turns and 25 related GNRA-type U-turns in a large set of aligned 16 S and 23 S rRNA sequences. U-turn candidates occur in hairpin loops (34 times) as well as in internal and multi-stem loops (24 times). These are classified into ten families based on loop type, sequence pattern (UNR or GNRA) and the nature of the closing base juxtaposition. In 13 cases, the bases on the 3' side of the turn, or on the immediate 5' side, are involved in tertiary covariations, making these sites strong candidates for tertiary interactions.

Calculating nucleic acid secondary structure

New results for calculating nucleic acid secondary structure by free energy minimization and phylogenetic comparisons have recently been reported. A complete set of DNA energy parameters is now available and the RNA parameters have been improved. Although databases of RNA secondary structures are still derived and expanded using computer-assisted, ad hoc comparative analysis, a number of new computer algorithms combine covariation analysis with energy methods.