tutorial 11
DESCRIPTION
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/. - PowerPoint PPT PresentationTRANSCRIPT
![Page 1: Tutorial 11](https://reader035.vdocument.in/reader035/viewer/2022081603/56813ff3550346895dab0f3e/html5/thumbnails/1.jpg)
Tutorial 11
RNA Structure Prediction
![Page 2: Tutorial 11](https://reader035.vdocument.in/reader035/viewer/2022081603/56813ff3550346895dab0f3e/html5/thumbnails/2.jpg)
RNA Structure Prediction
• Rfam – RNA structures database
• RNAfold – RNA secondary structure prediction
• tRNAscan – tRNA prediction
• TargetScan – microRNA prediction
![Page 3: Tutorial 11](https://reader035.vdocument.in/reader035/viewer/2022081603/56813ff3550346895dab0f3e/html5/thumbnails/3.jpg)
RNA secondary structure types
![Page 4: Tutorial 11](https://reader035.vdocument.in/reader035/viewer/2022081603/56813ff3550346895dab0f3e/html5/thumbnails/4.jpg)
Rfam
• The Rfam database is a collection of RNA families, each represented by multiple sequence alignments and consensus secondary structures.
http://rfam.sanger.ac.uk/
![Page 5: Tutorial 11](https://reader035.vdocument.in/reader035/viewer/2022081603/56813ff3550346895dab0f3e/html5/thumbnails/5.jpg)
Rfam
Different search modes
http://rfam.sanger.ac.uk/
![Page 6: Tutorial 11](https://reader035.vdocument.in/reader035/viewer/2022081603/56813ff3550346895dab0f3e/html5/thumbnails/6.jpg)
Search Rfam
![Page 7: Tutorial 11](https://reader035.vdocument.in/reader035/viewer/2022081603/56813ff3550346895dab0f3e/html5/thumbnails/7.jpg)
The secondary structure
![Page 8: Tutorial 11](https://reader035.vdocument.in/reader035/viewer/2022081603/56813ff3550346895dab0f3e/html5/thumbnails/8.jpg)
:::::: free extremes((())) Stem<<<>>> Internal Stem______ Loop,,,,,, Internal loop
Structure representations
![Page 9: Tutorial 11](https://reader035.vdocument.in/reader035/viewer/2022081603/56813ff3550346895dab0f3e/html5/thumbnails/9.jpg)
RNA secondary structure prediction
GGGCUAUUAGCUCAGUUGGUUAGAGCGCACCCCUGAUAAGGGUGAGGUCGCUGAUUCGAAUUCAGCAUAGCCCA
![Page 10: Tutorial 11](https://reader035.vdocument.in/reader035/viewer/2022081603/56813ff3550346895dab0f3e/html5/thumbnails/10.jpg)
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.
![Page 11: Tutorial 11](https://reader035.vdocument.in/reader035/viewer/2022081603/56813ff3550346895dab0f3e/html5/thumbnails/11.jpg)
![Page 12: Tutorial 11](https://reader035.vdocument.in/reader035/viewer/2022081603/56813ff3550346895dab0f3e/html5/thumbnails/12.jpg)
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).
![Page 13: Tutorial 11](https://reader035.vdocument.in/reader035/viewer/2022081603/56813ff3550346895dab0f3e/html5/thumbnails/13.jpg)
RNAfold - input
RNA sequence
![Page 14: Tutorial 11](https://reader035.vdocument.in/reader035/viewer/2022081603/56813ff3550346895dab0f3e/html5/thumbnails/14.jpg)
Minimal free energy structure
Frequency of the structure
RNAfold - output
Sequence-Structure alignment
Best “average” structure
![Page 15: Tutorial 11](https://reader035.vdocument.in/reader035/viewer/2022081603/56813ff3550346895dab0f3e/html5/thumbnails/15.jpg)
Graphic representation
Coloring options
![Page 16: Tutorial 11](https://reader035.vdocument.in/reader035/viewer/2022081603/56813ff3550346895dab0f3e/html5/thumbnails/16.jpg)
Minimal free energy
Folding Probabiliy
Mountain Plot
• 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.
centroid
Entropy
![Page 17: Tutorial 11](https://reader035.vdocument.in/reader035/viewer/2022081603/56813ff3550346895dab0f3e/html5/thumbnails/17.jpg)
RNAfold structure representations
![Page 18: Tutorial 11](https://reader035.vdocument.in/reader035/viewer/2022081603/56813ff3550346895dab0f3e/html5/thumbnails/18.jpg)
RNAalifold - input
Alignment
![Page 19: Tutorial 11](https://reader035.vdocument.in/reader035/viewer/2022081603/56813ff3550346895dab0f3e/html5/thumbnails/19.jpg)
RNAalifold - output
![Page 20: Tutorial 11](https://reader035.vdocument.in/reader035/viewer/2022081603/56813ff3550346895dab0f3e/html5/thumbnails/20.jpg)
RNAinverse - input
![Page 21: Tutorial 11](https://reader035.vdocument.in/reader035/viewer/2022081603/56813ff3550346895dab0f3e/html5/thumbnails/21.jpg)
RNAinverse - output
![Page 22: Tutorial 11](https://reader035.vdocument.in/reader035/viewer/2022081603/56813ff3550346895dab0f3e/html5/thumbnails/22.jpg)
tRNAscan
• Detection of tRNA genes in raw genomic sequence orother types of sequence inputs.
http://lowelab.ucsc.edu/tRNAscan-SE/
![Page 23: Tutorial 11](https://reader035.vdocument.in/reader035/viewer/2022081603/56813ff3550346895dab0f3e/html5/thumbnails/23.jpg)
Sequence
tRNAscan - input
![Page 24: Tutorial 11](https://reader035.vdocument.in/reader035/viewer/2022081603/56813ff3550346895dab0f3e/html5/thumbnails/24.jpg)
tRNAscan - results
![Page 25: Tutorial 11](https://reader035.vdocument.in/reader035/viewer/2022081603/56813ff3550346895dab0f3e/html5/thumbnails/25.jpg)
http://trna.nagahama-i-bio.ac.jp/
![Page 26: Tutorial 11](https://reader035.vdocument.in/reader035/viewer/2022081603/56813ff3550346895dab0f3e/html5/thumbnails/26.jpg)
MicroRNA - reminder
![Page 27: Tutorial 11](https://reader035.vdocument.in/reader035/viewer/2022081603/56813ff3550346895dab0f3e/html5/thumbnails/27.jpg)
• 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
![Page 28: Tutorial 11](https://reader035.vdocument.in/reader035/viewer/2022081603/56813ff3550346895dab0f3e/html5/thumbnails/28.jpg)
![Page 29: Tutorial 11](https://reader035.vdocument.in/reader035/viewer/2022081603/56813ff3550346895dab0f3e/html5/thumbnails/29.jpg)
![Page 30: Tutorial 11](https://reader035.vdocument.in/reader035/viewer/2022081603/56813ff3550346895dab0f3e/html5/thumbnails/30.jpg)
![Page 31: Tutorial 11](https://reader035.vdocument.in/reader035/viewer/2022081603/56813ff3550346895dab0f3e/html5/thumbnails/31.jpg)
Mir 31 - broadly conserved* microRNA
* conserved across most vertebrates, usually to zebrafish
![Page 32: Tutorial 11](https://reader035.vdocument.in/reader035/viewer/2022081603/56813ff3550346895dab0f3e/html5/thumbnails/32.jpg)
Mir 136 - conserved* microRNA
* conserved across most mammals, but usually not beyond placental mammals