Comparative sequence analysis as a tool for studying the secondary structure of mRNAs

Hidden messages in the nef gene of human immunodeficiency virus type 1 suggest a novel RNA secondary structure

The coexistence of multiple codes in the genome of human immunodeficiency virus type 1 (HIV-1) was analyzed. We explored factors constraining the variability of the virus genome primarily in relation to conserved RNA secondary structures overlapping coding sequences, and used a simple combination of algorithms for RNA secondary structure prediction based on the nearest-neighbor thermodynamic rules and a statistical approach. In our previous study, we applied this combination to a non- redundant data set of env nucleotide sequences, confirmed the conservative secondary structure of the rev-responsive element (RRE) and found a new RNA structure in the first conserved (C1) region of the env gene. In this study, we analyzed the variability of putative RNA secondary structures inside the nef gene of HIV-1 by applying these algorithms to a non-redundant data set of 104 nef sequences retrieved from the Los Alamos HIV database, and predicted the existence of a novel functional RNA secondary structure in the beta3/beta4 regions of nef. The predicted RNA fold in the beta3/beta4 region of nef appears in two forms with different loop sizes. The loop of the first fold consists of seven nucleotides (positions 494-500), with consensus UCAAGCU appearing in 79% of sequences. The other has a five-base loop (positions 495-499) with consensus CAAGC. The difference in size between these two loops may reflect the difference between respective counterparts in the hairpin recognition. This may also have an adaptive biological significance.

Pattern analysis of RNA secondary structure similarity and consensus of minimal-energy folding

We describe an automated procedure to search for consensus structures or substructures in a set of homologous or related RNA molecules. The procedure is based on the calculation of optimal and sub-optimal secondary structures using thermodynamic rules for base-pairing by energy-minimization. A linear representation of the secondary structures of the related RNAs is used so that they can be compared and classified using standard alignment and clusterings programs. We illustrate the method by means of two sets of homologous small RNAs, U2 and U3, and a set of alpha-globin mRNAs and show that biologically interesting consensus structures are obtained.

Chemical modification as a tool for analysis of messenger RNA secondary structure in ribonucleoprotein particles

Chemical modification of unpaired bases is demonstrated in this study to be a reliable method for determining the conformation of nucleotides in mRNA. The modified nucleotides are identified by primer extension using reverse transcriptase. We have used this procedure to compare the structure of limited regions of SV40 T-antigen mRNA in solution, in nonpolysome-bound cytoplasmic messenger ribonucleoprotein particles, and in nuclear ribonucleoprotein complexes. The results indicate that SV40 T-antigen mRNA adopts a specific structure both in solution and when complexed with cellular proteins. The structures adopted by the mRNA in solution and in native cellular protein particles are very similar.

Secondary structure model for mouse beta Maj globin mRNA derived from enzymatic digestion data, comparative sequence and computer analysis

A model for the secondary structure of mouse beta Maj globin messenger RNA is presented based on enzymatic digestion data, comparative sequence and computer analysis. Using 5'-32P-end-labeled beta globin mRNA as a substrate, single-stranded regions were determined with S1 and T1 nucleases and double-stranded regions with V1 ribonuclease from cobra venom. The structure data obtained for ca. 75% of the molecule was introduced into a computer algorithm which predicts secondary structures of minimum free energy consistent with the enzymatic data. Two prominent base paired regions independently derived by phylogenetic analysis were also present in the computer generated structure lending support for the model. An interesting feature of the model is the presence of long-range base pairing interactions which permit the beta globin mRNA to fold back on itself, thereby bringing the 5'- and 3'-noncoding regions within close proximity. This feature is consistent with data from other laboratories suggesting an interaction of the 5'- and 3'-domains in the mammalian globin mRNAs.