380 likes | 679 Views
Tutorial 11. RNA Structure Prediction. RNA Structure Prediction. Rfam – RNA structures database RNAfold – RNA secondary structure prediction tRNAscan – tRNA prediction TargetScan – microRNA prediction. RNA secondary structure types. Rfam. http://rfam.sanger.ac.uk/.
E N D
Tutorial 11 RNA Structure Prediction
RNA Structure Prediction • Rfam – RNA structures database • RNAfold – RNA secondary structure prediction • tRNAscan – tRNA prediction • TargetScan – microRNA prediction
Rfam http://rfam.sanger.ac.uk/ • The Rfam database is a collection of RNA families, each represented by multiple sequence alignments andconsensus secondary structures.
Rfam http://rfam.sanger.ac.uk/ Different search modes
Structure representations :::::: free extremes ((())) Stem <<<>>> Internal Stem ______ Loop ,,,,,, Internal loop
RNA secondary structure prediction GGGCUAUUAGCUCAGUUGGUUAGAGCGCACCCCUGAUAAGGGUGAGGUCGCUGAUUCGAAUUCAGCAUAGCCCA
RNA structure prediction by Vienna RNA package RNAfold server minimum free energy structures and base pair probabilities from single RNA or DNA sequences. RNAalifold server consensus secondary structures from an alignment of several related RNA or DNA sequences. You need to upload an alignment. RNAinverse server design RNA sequences for any desired target secondary structure.
RNAfold • Gives best stabilized structure (structure with minimal free energy (MFE)) • uses a dynamic programming algorithm that exploits base pairing and thermodynamic probabilities in order to predict the most likely structures of an RNA molecule (partition function algorithm).
RNAfold - input RNA sequence
RNAfold - output Minimal free energy structure Sequence-Structure alignment Frequency of the structure Best “average” structure
Graphic representation Coloring options
Minimal free energy Mountain Plot centroid • A mountain plot represents a secondary structure in a plot of height versus position. • The height m(k) is given by the number of base pairs enclosing the base at position k. • Loops correspond to plateaus and stems correspond to slopes. • The closer the two curves, the better defined the structure. Folding Probabiliy Entropy
RNAalifold - input Alignment
tRNAscan http://lowelab.ucsc.edu/tRNAscan-SE/ • Detection of tRNA genes in raw genomic sequence or • other types of sequence inputs.
tRNAscan - input Sequence
Search for predicted microRNA targets in mammals (/worm/fly) 3’ UTRs. • Find conserved 8mer and 7mer sites that match the seed region of each miRNA. • Predictions are ranked based on the predicted efficacy of targeting as calculated using the context+ scores of the sites
Mir 31 - broadly conserved* microRNA * conserved across most vertebrates, usually to zebrafish
Mir 136 - conserved* microRNA * conserved across most mammals, but usually not beyond placental mammals