本篇文章為大家展示了怎么使用IRESfinder分析環(huán)狀RNA的蛋白編碼潛能，內容簡明扼要并且容易理解，絕對能使你眼前一亮，通過這篇文章的詳細介紹希望你能有所收獲。

我們提供的服務有：網(wǎng)站設計制作、成都網(wǎng)站設計、微信公眾號開發(fā)、網(wǎng)站優(yōu)化、網(wǎng)站認證、通化縣ssl等。為上千余家企事業(yè)單位解決了網(wǎng)站和推廣的問題。提供周到的售前咨詢和貼心的售后服務，是有科學管理、有技術的通化縣網(wǎng)站制作公司

真核生物mRNA的5’ 帽子結構可以介導核糖體的結合，從而開啟翻譯過程，circRNA由于是一個閉合的環(huán)狀結構，缺失了5’端帽子結構，所以歸類為非編碼RNA的一種。

隨著環(huán)狀RNA研究的深入，有科學家發(fā)現(xiàn)部分環(huán)狀RNA可以編碼蛋白。除了常見的5’帽子結構介導核糖體結合外，還存在了一種特殊情況，在一些基因上存在一段長度在150-250bp的序列，這些序列能夠折疊成類似tRNA的結構，介導核糖體與RNA結合，起始蛋白質的翻譯。這樣的位點稱之為內部核糖體進入位點序列,Internal ribosome entry site, 簡稱IRES。

能夠編碼蛋白的環(huán)狀RNA上就是通過IRES來實現(xiàn)翻譯過程的，為例研究環(huán)狀RNA的蛋白編碼潛能，有學者開發(fā)了識別IRES的軟件，IRESfinder就是這樣一款軟件，用于識別真核生物的IRES位點。

利用實驗驗證過的583個人類的IRES序列，挑選了19個kmer用于區(qū)分IRES和非IRES的序列，在論文中，給出了測試數(shù)據(jù)集中二種序列的kmer分布

在上圖中，只包含了18個kmer,另外還有一個kmer是T，可以看到在IRES和非IRES序列之間，這些kmer的頻率分布是有差異的。

軟件采用Python進行開發(fā)，基本用法如下

python IRESfinder.py -f circRNA.fa -o IRES.out.xls

-f參數(shù)指定輸入的fasta格式的序列，-o參數(shù)指定輸出的結果文件。輸出文件的內容示意如下

ID Index Score
hsa_circ_0018046 IRES 0.817344005
hsa_circ_0039868 IRES 0.795083668
hsa_circ_0089160 IRES 0.53322605
hsa_circ_0048972 IRES 0.784080068
hsa_circ_0018658 IRES 0.745230164
hsa_circ_0067857 IRES 0.704497116
hsa_circ_0019137 IRES 0.742607966
hsa_circ_0063162 IRES 0.738372532
hsa_circ_0087609 IRES 0.793042932
hsa_circ_0006254 IRES 0.64587644

通過這個軟件，可以快速分析得到RNA上的IRES序列信息，不過軟件的結果中假陽性肯定是很高的，后續(xù)還要通過實驗手段來驗證。

上述內容就是怎么使用IRESfinder分析環(huán)狀RNA的蛋白編碼潛能，你們學到知識或技能了嗎？如果還想學到更多技能或者豐富自己的知識儲備，歡迎關注創(chuàng)新互聯(lián)行業(yè)資訊頻道。