利用一級與三級結構進行兩個RNA的比對

Pairwise Alignment of RNA Using Primary and Tertiary Structures

施承廷(Shih, Cheng Ting)

指導教授 ： 盧錦隆

國立清華大學/電機資訊學院/資訊工程學系/碩士(2015年)

摘要

近年來，存放在PDB資料庫中的RNA 三級（3D）結構數量持續地增加中。由於RNA分子的三級結構在演化的速度上比一級序列來的慢，所以RNA的結構比對可以對RNA分子的功能與演化關係帶來更具意義的了解，而這是無法藉由只分析一級序列資訊所能夠偵側出來。在2010年，我們實驗室開發出一套名叫iPARTS的軟體工具，可以去比對兩個RNA三級結構並決定出他們之間的結構相似程度。iPARTS的基本步驟如下：首先利用二個假扭轉角(pseudo-torsion angles) η 與θ把PDB資料庫中RNA結構的核苷酸畫在二維的平面以得到一個Ramachandran-like圖。接著使用親合性互動式(affinity propagation)分群演算法對Ramachandran-like圖上的核苷酸進行分群，以得到一個含有23個核苷酸形構的結構字元集(structural alphabet)。再來利用這個結構字元集把RNA 3D結構轉成1D的結構字元序列。最後使用傳統的序列比對演算法去比較二條結構字元編碼的序列，以決定出他們的結構相似程度。事實上，自從2010年我們發表iPARTS之後，許多新的RNA 3D結構已經被存放在PDB資料庫中。除此之外，最近有一個研究指出核苷酸一級序列對RNA結構的比對可以提供有用的額外資訊。因此，在本研究中我們首先使用目前可以從PDB資料庫中取得的RNA 3D結構去更新目前的iPARTS。接著我們設計出兩個把RNA 一級序列的資訊結合到iPARTS中的方法，使得iPARTS在同時考慮RNA一級序列和三級結構的資訊下，可以更準確地去比對兩個RNA的3D結構。最後我們的實驗證實了新版本的iPARTS，我們稱之為iPARTS2，確實比其之前的版本在RNA結構比對與功能設定上有比較好的表現，除此之外，我們提出的兩個把RNA 一級序列的資訊結合到iPARTS2中的方法，在RNA功能設定上，又更進一步比只考慮三級結構的iPARTS2還要好。

關鍵字

生物資訊 ； 核醣核酸三級結構 ； 結構比對 ； 結構字元 ； 親合性互動式

並列摘要

Recently, the number of RNA tertiary (3D) structures deposited in the PDB database continues to grow. Since structures of molecules evolve more slowly than their sequences, the structural comparison of RNAs can bring more sig-nificant insights into their functions and evolutionary relationships that would not be detected by analyzing sequence information alone. In 2010, our labor-atory developed a tool, called iPARTS, that aligns two RNA 3D structures and determines their structural similarities. The basic steps of our iPARTS are as follows. First, a Ramachandran-like diagram of RNAs was derived by plotting nucleotides of RNA structures in the PDB database on a 2D axis us-ing their two pseudo-torsion angles η and θ. Then, affinity propagation clus-tering algorithm was applied to the η-θ plot to obtain a structural alphabet (SA) of 23 nucleotide conformations. Next, the SA was used to transform RNA 3D structures into 1D sequences of SA letters. Finally, classical se-quence alignment methods were utilized to compare two SA-encoded se-quences and determine their structural similarities. In fact, many new RNA 3D structures have been deposited in the PDB database since we published iPARTS in 2010. In addition, there is a recent study to report that nucleotide 1D sequences of RNAs can provide useful additional information for their structural alignment. In this study, therefore, we first update our current iPARTS by utilizing currently available RNA 3D structures in the PDB data-base and then design two methods to incorporate RNA 1D sequence infor-mation into our iPARTS such that it can more accurately align two RNA 3D structures based on both their sequences and structures. Finally, our experi-mental results demonstrate that the new version of iPARTS, named iPARTS2, indeed outperforms its previous version in terms of accuracies of RNA struc-ture alignment and functional assignment. In addition, both the methods we propose in this study to incorporate 1D sequence information into iPARTS2 further outperform iPARTS2 with considering only 3D structures of RNAs in terms of accuracy of functional assignment.

並列關鍵字

bioinformatics ； RNA tertiary structure ； structural alignment ； structural alphabet ； affinity propagation

參考文獻

[2] Kolodny,　R. and Linial,　N. (2004) Approximate protein structural alignment in poly-nomial time. Proceedings of the National Academy of Sciences of the United States of America, 101, 12201–12206.

[3] Dror, O., Nussinov, R. and Wolfson, H. J. (2005) ARTS: alignment of RNA tertiary structures. Bioinformatics, 21, 47–53.

[4] Dror, O., Nussinov, R. and Wolfson, H. J. (2006) The ARTS web server for aligning RNA tertiary structures. Nucleic Acids Research, 34, W412–W415.

[5] Ferr`e, F., Ponty, Y., Lorenz, W. A. and Clote, P. (2007) DIAL: a web server for the pairwise alignment of two RNA three-dimensional structures using nucleotide, dihedral angle and base-pairing similarities. Nucleic Acids Research, 35, W659–W668.

[6] Chang, Y. F., Huang, Y. L. and Lu, C. L. (2008) SARSA: a web tool for structural alignment of RNA using a structural alphabet. Nucleic Acids Research, 36, W19–W24.