A computer method for finding common base paired helices in aligned sequences: application to the analysis of random...

Area of Science:

• Bioinformatics
• Computational Biology
• Molecular Biology

Background:

• Identifying conserved secondary structures in RNA is crucial for understanding RNA function and evolution.
• Existing methods may be limited in analyzing large datasets of nucleotide sequences.

Purpose of the Study:

• To introduce a novel computer program for identifying conserved secondary structures in aligned single-stranded RNA sequences.
• To assess the program's utility in analyzing large biological sequences and comparing conserved helices in real RNA with random sequences.

Main Methods:

• The program utilizes hash tables to detect and categorize common base-paired helices at identical positions across multiple sequences.
• It analyzes conserved helices, including those with compensating base changes.
• The program was tested on random RNA sequences with compositions similar to Escherichia coli 16S ribosomal RNA.

Main Results:

• The program successfully identified conserved helices, including those with compensating base changes.
• Analysis of random sequences showed significantly fewer conserved helices compared to actual ribosomal RNA sequences.
• Conserved helices in 16S ribosomal RNA were found to be longer on average than those in comparable random sequences.

Conclusions:

• The developed program is effective for analyzing conserved secondary structures in large nucleotide sequence datasets.
• The findings suggest that conserved RNA secondary structures, particularly in ribosomal RNA, are not easily explained by random sequence generation.
• The program has potential applications in analyzing long non-ribosomal nucleotide sequences.