Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP4873486B2 - 核酸配列の二次構造を予測する方法、核酸配列の二次構造の予測装置及び核酸配列の二次構造の予測プログラム - Google Patents
[go: Go Back, main page]

JP4873486B2 - 核酸配列の二次構造を予測する方法、核酸配列の二次構造の予測装置及び核酸配列の二次構造の予測プログラム - Google Patents

核酸配列の二次構造を予測する方法、核酸配列の二次構造の予測装置及び核酸配列の二次構造の予測プログラム Download PDF

Info

Publication number
JP4873486B2
JP4873486B2 JP2007134529A JP2007134529A JP4873486B2 JP 4873486 B2 JP4873486 B2 JP 4873486B2 JP 2007134529 A JP2007134529 A JP 2007134529A JP 2007134529 A JP2007134529 A JP 2007134529A JP 4873486 B2 JP4873486 B2 JP 4873486B2
Authority
JP
Japan
Prior art keywords
nucleic acid
acid sequence
secondary structure
feature vector
evaluated
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2007134529A
Other languages
English (en)
Other versions
JP2008283943A (ja
Inventor
穣 秋冨
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
NEC Solution Innovators Ltd
Original Assignee
NEC Solution Innovators Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by NEC Solution Innovators Ltd filed Critical NEC Solution Innovators Ltd
Priority to JP2007134529A priority Critical patent/JP4873486B2/ja
Priority to DE102008024282A priority patent/DE102008024282B4/de
Priority to US12/124,467 priority patent/US8370069B2/en
Publication of JP2008283943A publication Critical patent/JP2008283943A/ja
Application granted granted Critical
Publication of JP4873486B2 publication Critical patent/JP4873486B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B15/00ICT specially adapted for analysing two-dimensional [2D] or three-dimensional [3D] molecular structures, e.g. structural or functional relations or structure alignment
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B15/00ICT specially adapted for analysing two-dimensional [2D] or three-dimensional [3D] molecular structures, e.g. structural or functional relations or structure alignment
    • G16B15/10Nucleic acid folding

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Chemical & Material Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Biotechnology (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Theoretical Computer Science (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Biophysics (AREA)
  • Evolutionary Biology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biochemistry (AREA)
  • Molecular Biology (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Apparatus Associated With Microorganisms And Enzymes (AREA)
  • Information Retrieval, Db Structures And Fs Structures Therefor (AREA)

Description

本発明は、核酸配列の二次構造を予測する方法、並びにこの方法を実行する核酸配列の二次構造の予測装置及び核酸配列の二次構造の予測プログラムに関する。
DNA、RNAなどの核酸配列は、アデニン(A)、シトシン(C)、グアニン(G)又はチミン(T)若しくはウラシル(U)の4種類の塩基からなる分子である。この核酸配列の分子の機能は、その配列の二次構造に大きく依存することが知られ、この依存性は、エクソンなどのいわゆる構造遺伝子に代表される機能性核酸だけでなく、イントロンなどの非構造遺伝子や、セントラルドグマに準じて細胞内で内在的に転写され得ないアプタマーなどの各種の核酸配列においても、発揮される性質である。
核酸配列の二次構造を予測することは、非常に重要であるが、同様の配列を有していても、異なる構造を与えることもあり、異なる配列を有していても、同様又は類似の構造を与えることもあり、単に核酸配列の二次構造を予測した場合であっても、その構造が必ずしもその配列において取り得る最適な二次構造を与えているとは言えない場合も、数多く考えられる。実際、アプタマーと呼ばれる核酸配列など、人工的に合成される核酸配列に関しては、核酸配列の特定の重要な領域のみによって生物学的な機能を発揮する場合が多く、その特定の領域については構造の類似性が期待されるが、その他の領域に関しては必ずしも構造の類似性があるとは限らない。そのため、特定の核酸配列について得た各二次構造について、他の核酸配列の二次構造との間の類似性を評価することで、核酸配列の二次構造予測の精度を高める試みがなされている。
従来の核酸配列の二次構造の類似性比較の方法としては、進化的に保存されている核酸配列について、核酸配列を構成する塩基のアラインメントと核酸配列の二次構造の予測とを組み合わせたものによって評価するものが一般的である。また、核酸配列中の塩基対の形成確率などを独自に算出し、それを基にして予測構造を実際の構造に近い形に修正、再構築するような幾つかのアルゴリズムが報告されている(非特許文献1)。
だが、このような方法で核酸配列の二次構造の類似性を評価した場合には、基本的には構造の全体としての類似性を評価することになるため、全体的には類似でなくとも局所的に類似の構造を取り得る核酸配列の構造の類似性を的確に評価することは困難であった。
Gorodkinら著、Nucleic Acids Research、1997年、25巻、p.3724〜3732
本発明は、上記の問題に鑑みて提案されたものであって、核酸配列の二次構造の全体的な類似性のみならず、局所的な類似性を評価し得る、核酸配列の二次構造を予測する方法、並びにこの方法を実行する核酸配列の二次構造の予測装置及び予測プログラムを提供することを目的とする。
本発明による核酸配列の二次構造を予測する方法は、核酸配列の二次構造のうち、該二次構造の構造要素を抽出する工程と;前記構造要素の特徴ベクトルに基づいて、評価の対象である核酸配列のそれぞれの構造要素の類似性を評価する工程と;を有することを特徴とする。
本発明による核酸配列の二次構造の予測装置は、核酸配列の二次構造のうち、該二次構造の構造要素を抽出する構造要素抽出手段と、前記構造要素の特徴ベクトルに基づいて、評価の対象である核酸配列のそれぞれの構造要素の類似性を評価する局所構造類似性評価手段と、を有することを特徴とする。
本発明による核酸配列の二次構造の予測プログラムは、核酸配列の二次構造のうち、該二次構造の構造要素を抽出する工程と;
前記構造要素の特徴ベクトルに基づいて、評価の対象である核酸配列のそれぞれの構造要素の類似性を評価する工程と;
を実行することを特徴とする。
本発明によれば、局所的、全体的を問わず、任意のウィンドウサイズで配列の構造間の類似性の評価を行うことが可能である。
また、同じ機能を持つ核酸配列は一般的に構造が類似であるため、この発明に寄れば、複数の配列を入力に用いることにより、単独の配列からよりも更に正確な二次構造の予測が可能である。
さらに、局所的に類似な構造の局所的に類似な箇所が判明した場合、それがその核酸の機能にとって重要な領域である可能性が高い。したがって、この発明に寄れば、核酸配列中の機能的に重要な領域を予測することが可能である。
本発明は、核酸配列の二次構造の予測で得た、各構造の要素を抽出して、要素の特徴を示すベクトルとして表現し、対象となる核酸配列の二次構造と、被対象である核酸配列の二次構造とを、それぞれの構造に由来するベクトルの相関を計算することによって、核酸配列に含まれる各構造の種類の影響を最小限に抑えて、核酸配列の二次構造の類似性を検討するものである。
以下、本発明による核酸配列の二次構造の予測装置の概略図を参照しながら、本発明による核酸配列の二次構造を予測する方法、核酸配列の二次構造の予測装置、及び核酸配列の二次構造の予測プログラムについて詳細に説明する。
なお、本発明による核酸配列の二次構造の予測装置は、本発明による核酸配列の二次構造を予測する方法を実行する装置であり、本発明による核酸配列の二次構造の予測プログラムは、本発明による核酸配列の二次構造を予測する方法を実行するプログラムである。また、本発明において、「核酸配列」とは、DNA、RNAなど、種々の遺伝子の配列をいう。
図1は、本発明による核酸配列の二次構造の予測装置の概略図を示す。本発明による核酸配列の二次構造の予測装置は、キーボードなどの入力装置1と、プログラム制御により動作するデータ処理装置2と、情報を記憶する記憶装置3と、ディスプレイ装置や印刷装置などの出力装置4とを有する。
データ処理装置2は、構造要素抽出手段21と、局所構造類似性評価手段22とを備える。
構造要素抽出手段21は、まず、入力装置1から入力された、二次構造の予測の対象とする核酸配列と、その核酸配列の二次構造の構造候補の情報とを参照して、構造候補の各構造を、ステム構造と、ステム構造を構成する塩基以外の塩基で構成された周辺構造とに分割し、ステム構造及び周辺構造の構造要素を抽出する。構造要素抽出手段21は、次に、この構造要素を数値化して、特徴ベクトルを得る。特徴ベクトルには、数値化された構造要素の特徴がベクトルの形で並列に保持される。特徴ベクトルの集合として保持された核酸配列の構造情報は、構造要素記憶部31に記憶される。
局所構造類似性評価手段22は、構造要素記憶部31に記憶された、特徴ベクトルの集合として保持された核酸配列の構造候補の情報を取り出し、一の核酸配列の各特徴ベクトルと、他の核酸配列の各特徴ベクトルとの類似性を評価する。これにより、予測の対象とした核酸配列の核酸配列の構造の局所的な類似性を評価し得る。各核酸配列は、複数の構造候補情報、及び各構造候補ごとに複数の特徴ベクトルを保持しているため、核酸配列の構造の局所的な類似性の評価は、例えば各核酸配列の保持する特徴ベクトルの総当りによる類似性の評価によって決定される。
ここで、局所構造類似性評価手段22は、上記の通り構造要素記憶部31に記憶された特徴ベクトルを参照して、ひとつの構造候補から分割されたステム構造と周辺構造とに由来する特徴ベクトルを、それぞれ類似性の評価の対象として用いてもよく、ひとつのステム構造とこのステム構造の周囲に存在するひとつの周辺構造とに由来する各特徴ベクトルを組み合わせて、組み合わせて得た特徴ベクトルを、上記の類似性の評価の対象として用いてもよい。なお、特徴ベクトルを組み合わせる方法としては、ベクトルの加減乗除等の公知の算術方法を用いればよい。
各特徴ベクトルを組み合わせて用いる場合、核酸配列は、その構造候補の一つを基としてステム・ループ単位に分割され、ステム・ループの集合として取り扱われる。つまり、一つの核酸配列について複数の構造候補が存在する場合、一つの核酸配列について複数のステム・ループの集合が存在することとなる。また、単位となったステム・ループは、その構造要素に応じた特徴ベクトルの形で情報を保持する。
また、局所構造類似性評価手段22において採用され得る、特徴ベクトルの類似性の評価としては、ベクトル間の相関を評価し得る方法であれば、特に制約はない。なかでも、ベクトルの相関解析において採用されている内積の計算、サポートベクターマシーン(SVM)などによる各種パターン分類アルゴリズムの活用などの方法が挙げられる。
記憶装置3は、構造要素記憶部31を備える。
構造要素記憶部31は、構造要素抽出手段21によって抽出された、核酸配列の特徴ベクトルの情報を記憶する。ここで、特徴ベクトルに格納されるパラメータとしては、構造要素の特徴を示すものであれば、特に制約はなく、構造要素を構成する塩基の種類及び数、構造要素で取り得る自由エネルギーを始め、ステム構造やループ構造といった、構造の種類の別であってもよい。
(本発明による核酸配列の二次構造を予測する方法の各工程、並びに本発明による核酸配列の二次構造の予測装置及び核酸配列の二次構造の予測プログラムの動作等)
次に、図1の概略図、並びに図2乃至図4のフローチャートを参照して、本発明による核酸配列の二次構造を予測する方法の各工程、並びに本発明による核酸配列の二次構造の予測装置及び核酸配列の二次構造の予測プログラムの動作等について詳細に説明する。
なお、図2は、本発明による核酸配列の二次構造を予測する方法の各工程を示すフローチャートであり、図3は、構造要素を抽出する工程を示すフローチャートであり、図4は、各構造の特徴ベクトルに基づいて、構造類似性を評価する工程を示すフローチャートである。
本発明による核酸配列の二次構造の予測装置において、入力装置1から与えられた核酸配列、及びこの核酸配列の構造候補の情報は、構造要素抽出手段21に送られる(A1、A2、A3)。なお、核酸配列の二次構造の構造候補は、ステップA2のように、入力されてもよく、また、本発明による核酸配列の二次構造の予測装置に備えた、核酸配列の二次構造を予測し得る公知の二次構造予測手段から得た二次構造の構造候補を用いてもよい。
構造要素抽出手段21は、まず、最初の核酸配列の最初の構造候補の情報を参照して、ステム構造の単位に分割する(A31)。次に、構造要素抽出手段21は、構造候補から分割されたステム構造について、分析する(A32)。このステム構造の分析としては、このステム構造の種々の要素(本発明において、ステム構造の構造要素とも称する。)を抽出して、それぞれの構造要素を数値化することで、行う。このステム構造の構造要素としては、ステム構造の長さ、ステム構造を構成する塩基対の種類、塩基対の出現パターンが挙げられる。その後、構造要素抽出手段21は、ステム構造を構成する塩基以外の塩基で構成される周辺構造の種々の要素(本発明において、周辺構造の構造要素とも称し、上記のステム構造の構造要素とまとめて、単に構造要素とも称する。)について、ステップA32と同様に、周辺構造の構造要素を抽出してそれぞれの構造要素を数値化する分析を行う(A33)。この周辺構造としては、ステム構造を構成する塩基以外の塩基で構成される構造であれば、特に限定はなく、例えば、バルジループ構造、内部ループ構造、ヘアピンループ構造、多岐ループ構造、一本鎖末端などのループ構造や、シュードノット、Gカルテットなどの特殊な構造など、塩基配列の構造として公知の構造が挙げられる。また、この周辺構造の構造要素としては、周辺構造を構成する塩基の数、周辺構造を構成する塩基の種類、塩基の出現パターンなどが挙げられる。
このように、ステム構造及び周辺構造の各構造要素を数値化して得た値は、ステム構造及び周辺構造の特徴を示すものとして、特徴ベクトルとして保持される。
次に、構造要素抽出手段21は、最初の核酸配列に、未抽出の構造パターンがあるかどうかを判定する(A34)。未抽出の構造パターンがあると判定された場合には、その構造をセットして(A35)、上記のステップA31乃至A33までの各ステップを、最初の核酸配列に未抽出の構造パターンがなくなるまで、繰り返す。その後、構造要素抽出手段21は、特徴ベクトルを抽出していない核酸配列があるかどうかを判定し(A36)、未抽出の核酸配列があると判定された場合には、その配列をセットして(A37)、上記のステップA31乃至A35の各ステップを、未抽出の核酸配列がなくなるまで、繰り返す。核酸配列及びその構造候補からの特徴ベクトルの抽出が終了すると、それらの特徴ベクトルは構造要素記憶部31へ格納され(A38)、構造要素を抽出する工程は、終了する(A3)。
次に、本発明による核酸配列の二次構造の予測装置は、ステップA3の構造要素を抽出する工程が終了すると、ステップA3で得たステム構造及び周辺構造について得た各特徴ベクトルに基づいて、局所構造類似性評価手段22を用いて、ステップA4の構造候補の類似性の評価を行う。
局所構造類似性評価手段22は、まず、構造要素記憶部31から、評価の対象である、核酸配列の特徴ベクトルの情報を読み込む(A41)。次に、局所構造類似性評価手段22は、評価の被対象である、核酸配列の特徴ベクトルの情報を読み込む(A42)。その後、局所構造類似性評価手段22は、評価対象である特徴ベクトルと、被評価対象である特徴ベクトルとの間の類似性を評価して、評価対象である特徴ベクトルの候補評価値を得る(A43)。
ここで、各特徴ベクトルの類似性の評価の方法としては、二つのベクトルの内積の計算などのベクトルの相関解析が挙げられる。
次に、局所構造類似性評価手段22は、評価の被対象である特徴ベクトルのうち、類似性の評価の被対象とされていない特徴ベクトルがあるかどうかを判定する(A44)。被評価対象である特徴ベクトルのうち、未評価の特徴ベクトルがあると判定された場合には、局所構造類似性評価手段22は、その特徴ベクトルを被評価対象である特徴ベクトルとしてセットし(A45)、上記のステップA42乃至A43の各ステップを、被評価対象である特徴ベクトルのうち、未評価の特徴ベクトルがなくなるまで、繰り返す。このようにして、評価対象である特徴ベクトルについて、被評価対象である特徴ベクトルの全てに対する候補評価値が得られる。局所構造類似性評価手段22は、これらの候補評価値に基づいて、評価対象であるひとつの特徴ベクトルについての確定評価値を算出する(A46)。ここで、確定評価値としては、類似性を評価した各核酸配列について、それぞれ最高の類似性を示した特徴ベクトルとの候補評価値の平均や、一定以上の類似性があると評価された特徴ベクトルの個数が挙げられる。
次に、局所構造類似性評価手段22は、評価の対象である、核酸配列の特徴ベクトルが他にあるかどうかを判定する(A47)。評価対象の特徴ベクトルがあると判定された場合には、その特徴ベクトルをセットし(A48)、上記の通り、ステップA42乃至A46までの各ステップを、評価の対象である、核酸配列の特徴ベクトルがなくなるまで、繰り返す。
これら全ての特徴ベクトルについての類似性を示す確定評価値が得られることによって、ある核酸配列について、他の配列でも頻出する傾向にある特徴ベクトル、即ち局所的な構造の類似性が分かる。これらの評価結果の情報は、ディスプレイなどの出力装置4に出力される。
なお、上記の通り、本発明において、構造要素抽出手段21は、核酸配列の二次構造のうち、該二次構造の構造要素を抽出する機能を有し、局所構造類似性評価手段22は、構造要素の特徴ベクトルに基づいて、評価の対象である核酸配列のそれぞれの構造要素の類似性を評価する機能を有する。
まず、入力装置1より図5に示す、配列番号1及び2の配列、及びその構造が入力されたとすると、これは構造要素抽出手段21に送られる(A1、A2、A3)。
構造要素抽出手段21は、まず、構造1の構造情報を参照して、ステム構造の単位に分割する(A31)。構造1は、このステップにおいて、ステム1−1とステム1−2の二つのステム構造に分割される(A31)。それぞれのステム構造について、ステム、及びループの長さのみに着目した場合、ステム1−1のステムの長さは、3であり、ループの長さは4である。また、ステム1−2のステムの長さは、4であり、ループの長さは、4である。次に、構造要素抽出手段21は、これらの構造要素の特徴をそれぞれ、(3,4)、(4,4)と数値化し、これらを特徴ベクトルとして保持する(A32、A33)。未抽出の構造パターンとして、配列1には、構造1のほかに構造候補はなく(A34)、未抽出の配列として、配列1のほかに配列2があり(A36)、配列2には、構造2があるため、次に構造要素を抽出する対象として、配列2の構造2がセットされる(A37、A35)。次に、構造要素抽出手段21は、構造2について、構造1と同様にして、ステム2−1とステム2−2の二つのステム構造に分割し(A31)、それぞれのステム2−1及び2−2の構造要素の特徴を、構造1のステム構造と同様にして、それぞれ(4,4)、及び(2,4)と数値化し、これらを特徴ベクトルとして保持する(A32、A33)。未抽出の構造パターンとして、配列2には、構造2のほかに構造候補はなく(A34)、未抽出の配列として、配列2のほかに配列は存在しないため(A36)、構造要素抽出手段21は、ステム1−1乃至ステム2−2までの4つの特徴ベクトルを、構造要素記憶部31に格納し(A38)、構造要素を抽出する工程を終了する。
次に、局所構造類似性評価手段22は、構造要素記憶部31に記憶された4つの特徴ベクトルを読み込む(A41)。次に、局所構造類似性評価手段22は、最初の評価対象として、配列1の構造1のステム1−1の特徴ベクトルである(3,4)をセットし、最初の被評価対象の特徴ベクトルとして、配列2の構造2の一つ目のステムであるステム2−1の特徴ベクトルである(4,4)との類似性を評価する(A42、A43)。ここで、ベクトルの相関解析として、単純に二つのベクトルの内積を用いて類似性を評価する場合、この二つの特徴ベクトルの内積は、cosθ=0.9898となる(A43)。未評価の特徴ベクトルとして、配列2の構造2には、もう一つのステム構造であるステム2−2が存在するため(A44)、局所構造類似性評価手段22は、評価対象のステム構造であるステム1−1に対して、今度は被評価対象の特徴ベクトルとして、ステム2−2の特徴ベクトルである(2,4)との類似性を評価する(A45、A42)。ステム2−2についてもステム2−1と同様に類似性の評価を行うと、この二つのベクトルの内積は、cosθ=0.9839となる(A43)。
未評価の特徴ベクトルとして、配列2には他に特徴ベクトルはないので(A44)、局所構造類似性評価手段22は、評価対象である配列1の構造1のステム1−1についての候補評価値である「0.9898」及び「0.9839」に基づいて、このステム1−1の確定評価値を算出する(A46)。ここで、候補評価値のうち、最も大きな値を確定評価値とする場合、局所構造類似性評価手段22は、「0.9898」及び「0.9839」のうち、0.9898という値を算出して(A46)、この値を確定評価値として、ステム1−1についての類似性の評価を終了する。
次に、評価の対象である、核酸配列の特徴ベクトルとして、配列1には、ステム1−2という特徴ベクトルが存在するため(A47)、局所構造類似性評価手段22は、評価対象である、核酸配列の特徴ベクトルとして、ステム1−2の特徴ベクトルである(4,4)をセットし(A48)、上記のステップA42乃至A45と同様に、被評価対象である配列2の各ステムとの類似性の評価を行う。このようにして、ステム1−2と、ステム2−1及びステム2−2との類似性の評価の結果、ステム1−2の確定評価値として、特徴ベクトルの内積である1.0という値が得られる(A46)。
以下、配列2の構造2のステム2−1、及びステム2−2についても、同様の評価を行い、全ての特徴ベクトルの確定評価値が得られたところで、構造類似性の評価が終了する(A47)。この例の場合、配列1では、構造1のステム1−2について、配列2では構造2のステム2−1について、それぞれ最も大きな確定評価値である1.0が得られたので、局所構造類似性評価手段22は、それぞれの配列において、それぞれのステムが局所的な構造の類似性が高い構造であると予測する。
これらの評価結果は、必要に応じて出力装置4によってディスプレイなどに出力され(A4)、配列1及び2の予測が終了する。
以上、本発明の好適な実施の形態により本発明を説明した。ここでは特定の具体例を示して本発明を説明したが、特許請求の範囲に定義された本発明の広範な趣旨および範囲から逸脱することなく、これら具体例に様々な修正および変更を加えることができることは明らかである。すなわち、具体例の詳細および添付の図面により本発明が限定されるものと解釈してはならない。
本発明による核酸配列の二次構造の予測装置の概略図を示す。 本発明による核酸配列の二次構造を予測する方法の各工程を示すフローチャートである。 構造要素を抽出する工程を示すフローチャートである。 各構造の特徴ベクトルに基づいて、構造類似性を評価する工程を示すフローチャートである。 実施例1で用いた配列、及びその構造を示す。
符号の説明
1 入力装置
2 データ処理装置
3 記憶装置
4 出力装置
21 構造要素抽出手段
22 局所構造類似性評価手段
31 構造要素記憶部

Claims (31)

  1. 入力された評価対象の核酸配列の二次構造候補と、入力された被評価対象の核酸配列の二次構造候補とから、それぞれ、構造要素を抽出し、
    抽出した前記評価対象の各構造要素と、抽出した前記被評価対象の各構造要素とから、それぞれ特徴ベクトルを得る構造要素抽出工程と;
    記被評価対象の各構造要素の特徴ベクトルに対する、前記評価対象の各構造要素の特徴ベクトルの類似性を評価する局所構造類似性評価工程と;
    を有し、前記各工程が全てコンピュータ上で実行されることを特徴とする核酸配列の局所的な二次構造を予測する方法。
  2. 前記構造要素抽出工程は、
    前記二次構造候補構造単位分割し、前記構造単位の構造要素を抽出する工程を有することを特徴とする請求項1に記載の核酸配列の局所的な二次構造を予測する方法。
  3. 記核酸配列の二次構造候補を、前記構造単位としてステム構造と周辺構造とに分割することを特徴とする請求項に記載の核酸配列の局所的な二次構造を予測する方法。
  4. 前記構造要素抽出工程は、
    前記分割した前記ステム構造と前記周辺構造とから、
    前記ステム構造の構造要素と、前記周辺構造の構造要素とを抽出することを特徴とする請求項3に記載の核酸配列の局所的な二次構造を予測する方法。
  5. 前記構造要素抽出工程は、
    前記構造要素を数値化して、前記特徴ベクトルを得ることを特徴とする請求項1から4のいずれか一項に記載の核酸配列の局所的な二次構造を予測する方法。
  6. 前記特徴ベクトルが、前記ステム構造由来の特徴ベクトルと前記周辺構造由来の特徴ベクトルとを組み合わせた特徴ベクトルであることを特徴とする請求項から5のいずれか一項に記載の核酸配列の局所的な二次構造を予測する方法。
  7. 前記局所構造類似性評価工程において、
    前記類似性の評価を、ベクトルの相関解析により行うことを特徴とする請求項から6のいずれか一項に記載の核酸配列の局所的な二次構造を予測する方法。
  8. 前記局所構造類似性評価工程において
    前記被評価対象の各構造要素の特徴ベクトルに対する、前記評価対象の各構造要素の特徴ベクトルの類似性を数値化し、前記値を、前記評価対象の前記特徴ベクトルの候補評価値とし、
    前記候補評価値のうち、類似性が高い候補評価値を確定評価値とすることを特徴とする請求項1から7のいずれか一項に記載の核酸配列の局所的な二次構造を予測する方法。
  9. 前記被評価対象の核酸配列が、複数の特徴ベクトルを有する場合、
    前記局所構造類似性評価工程において、
    前記被評価対象の各特徴ベクトルに対する、前記評価対象の特徴ベクトルの類似性を数値化し、前記各値を、前記評価対象の特徴ベクトルの候補評価値とし、
    前記評価対象の特徴ベクトルの前記候補評価値のうち、類似性が高い候補評価値を確定評価値とすることを特徴とする請求項に記載の核酸配列の局所的な二次構造を予測する方法。
  10. 前記候補評価値のうち、最も高い類似性または一定以上の類似性を示す前記候補評価値を、前記確定評価値とすることを特徴とする請求項8または9に記載の核酸配列の局所的な二次構造を予測する方法。
  11. 前記確定評価値から、前記評価対象の核酸配列における、前記被評価対象の核酸配列に対して類似性が高い構造を予測することを特徴とする請求項から10のいずれか一項に記載の核酸配列の局所的な二次構造を予測する方法。
  12. 前記評価対象の核酸配列が、複数の特徴ベクトルを有する場合、
    前記局所構造類似性評価工程において、
    前記被評価対象の特徴ベクトルに対する前記評価対象の特徴ベクトルの類似性を数値化し、前記各値を、前記評価対象の特徴ベクトルの候補評価値とし、
    前記評価対象の特徴ベクトルの前記候補評価値のうち、類似性が高い候補評価値を確定評価値とすることを特徴とする請求項から11のいずれか一項に記載の核酸配列の局所的な二次構造を予測する方法。
  13. 前記核酸配列が、複数の二次構造候補を有する場合、
    前記構造要素抽出工程において、
    前記各二次構造候補から、それぞれ、前記構造要素を抽出することを特徴とする請求項1から12のいずれか一項に記載の核酸配列の局所的な二次構造を予測する方法。
  14. 前記構造要素抽出工程は、
    入力された前記評価対象の核酸配列と、入力された前記被評価対象の核酸配列とから、それぞれの二次構造候補を得る工程を含むことを特徴とする請求項1から13のいずれか一項に記載の核酸配列の局所的な二次構造を予測する方法。
  15. 前記評価対象の核酸配列が、核酸アプタマーの配列であることを特徴とする請求項1から14のいずれか一項に記載の核酸配列の局所的な二次構造を予測する方法。
  16. 評価対象の核酸配列の二次構造候補と、被評価対象の核酸配列の二次構造候補とを入力する入力手段と、
    前記評価対象の核酸配列の二次構造候補と、前記被評価対象の核酸配列の二次構造候補とから、それぞれ、構造要素を抽出し、抽出した前記評価対象の各構造要素と、抽出した前記被評価対象の各構造要素とから、それぞれ特徴ベクトルを得る構造要素抽出手段と、
    前記被評価対象の各構造要素の特徴ベクトルに対する、前記評価対象の各構造要素の特徴ベクトルの類似性を評価する局所構造類似性評価手段と、
    を有することを特徴とする核酸配列の局所的な二次構造の予測装置
  17. 前記構造要素抽出手段は、
    前記二次構造候補を構造単位に分割し、前記構造単位の構造要素を抽出することを特徴とする請求項16に記載の核酸配列の局所的な二次構造の予測装置
  18. 前記核酸配列の二次構造候補を、前記構造単位としてステム構造と周辺構造とに分割することを特徴とする請求項17に記載の核酸配列の局所的な二次構造の予測装置。
  19. 前記構造要素抽出手段は、
    前記分割した前記ステム構造と前記周辺構造とから、
    前記ステム構造の構造要素と、前記周辺構造の構造要素とを抽出することを特徴とする請求項18に記載の核酸配列の局所的な二次構造の予測装置。
  20. 前記構造要素抽出手段は、
    前記構造要素を数値化して、前記特徴ベクトルを得ることを特徴とする請求項16から19のいずれか一項に記載の核酸配列の局所的な二次構造の予測装置。
  21. 前記特徴ベクトルが、前記ステム構造由来の特徴ベクトルと前記周辺構造由来の特徴ベクトルとを組み合わせた特徴ベクトルであることを特徴とする請求項18から20のいずれか一項に記載の核酸配列の局所的な二次構造の予測装置。
  22. 前記局所構造類似性評価手段は、
    前記類似性の評価を、ベクトルの相関解析により行うことを特徴とする請求項16から21のいずれか一項に記載の核酸配列の局所的な二次構造の予測装置。
  23. 前記局所構造類似性評価手段は、
    前記被評価対象の各構造要素の特徴ベクトルに対する、前記評価対象の各構造要素の特徴ベクトルの類似性を数値化し、前記値を、前記評価対象の前記特徴ベクトルの候補評価値とし、
    前記候補評価値のうち、類似性が高いことを示す確定評価値を得ることを特徴とする請求項16から22のいずれか一項に記載の核酸配列の局所的な二次構造の予測装置。
  24. 前記被評価対象の核酸配列が、複数の特徴ベクトルを有する場合、
    前記局所構造類似性評価手段は、
    前記被評価対象の各特徴ベクトルに対する、前記評価対象の特徴ベクトルの類似性を数値化し、前記各値を、前記評価対象の特徴ベクトルの候補評価値とし、
    前記評価対象の特徴ベクトルの前記候補評価値のうち、類似性が高い候補評価値を確定評価値とすることを特徴とする請求項23に記載の核酸配列の局所的な二次構造の予測装置。
  25. 前記候補評価値のうち、最も高い類似性または一定以上の類似性を示す前記候補評価値を、前記確定評価値とすることを特徴とする請求項23または24に記載の核酸配列の局所的な二次構造の予測装置。
  26. 前記確定評価値から、前記評価対象の核酸配列における、前記被評価対象の核酸配列に対して類似性が高い構造を予測することを特徴とする請求項23から25のいずれか一項に記載の核酸配列の局所的な二次構造の予測装置。
  27. 前記評価対象の核酸配列が、複数の特徴ベクトルを有する場合、
    前記局所構造類似性評価手段は、
    前記被評価対象の特徴ベクトルに対する、前記評価対象の各特徴ベクトルの類似性を数値化し、前記各値を、前記評価対象の各特徴ベクトルの候補評価値とし、
    前記評価対象の特徴ベクトルの前記候補評価値のうち、類似性が高い候補評価値を確定評価値とすることを特徴とする請求項23から26のいずれか一項に記載の核酸配列の局所的な二次構造の予測装置。
  28. 前記核酸配列が、複数の二次構造候補を有する場合、
    前記構造要素抽出手段は、
    前記各二次構造候補から、それぞれ、前記構造要素を抽出することを特徴とする請求項16から27のいずれか一項に記載の核酸配列の局所的な二次構造の予測装置。
  29. 前記入力手段は、
    前記評価対象の核酸配列と、前記被評価対象の核酸配列とを入力する手段であり、
    前記構造要素抽出手段は、
    入力された前記評価対象の核酸配列と、入力された前記被評価対象の核酸配列とから、それぞれの二次構造候補を得ることを特徴とする請求項16から28のいずれか一項に記載の核酸配列の局所的な二次構造の予測装置。
  30. 前記評価対象の核酸配列が、核酸アプタマーの配列であることを特徴とする請求項16から29のいずれか一項に記載の核酸配列の局所的な二次構造の予測装置。
  31. 請求項1から15のいずれか一項に記載の核酸配列の局所的な二次構造を予測する方法を、コンピュータ上で実行可能なことを特徴とするプログラム。
JP2007134529A 2007-05-21 2007-05-21 核酸配列の二次構造を予測する方法、核酸配列の二次構造の予測装置及び核酸配列の二次構造の予測プログラム Active JP4873486B2 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2007134529A JP4873486B2 (ja) 2007-05-21 2007-05-21 核酸配列の二次構造を予測する方法、核酸配列の二次構造の予測装置及び核酸配列の二次構造の予測プログラム
DE102008024282A DE102008024282B4 (de) 2007-05-21 2008-05-20 Verfahren zur Vorhersage der Sekundärstruktur einer Nucleinsäuresequenz, Vorhersagevorrichtung für die Sekundärstruktur einer Nucleinsäuresequenz und Vorhersageprogramm zum Vorhersagen der Sekundärstruktur einer Nucleinsäuresequenz
US12/124,467 US8370069B2 (en) 2007-05-21 2008-05-21 Method for predicting secondary structure of nucleic acid sequence, a predictor for secondary structure of nucleic acid sequence and a predicting program for predicting secondary structure of nucleic acid sequence

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2007134529A JP4873486B2 (ja) 2007-05-21 2007-05-21 核酸配列の二次構造を予測する方法、核酸配列の二次構造の予測装置及び核酸配列の二次構造の予測プログラム

Publications (2)

Publication Number Publication Date
JP2008283943A JP2008283943A (ja) 2008-11-27
JP4873486B2 true JP4873486B2 (ja) 2012-02-08

Family

ID=39942315

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2007134529A Active JP4873486B2 (ja) 2007-05-21 2007-05-21 核酸配列の二次構造を予測する方法、核酸配列の二次構造の予測装置及び核酸配列の二次構造の予測プログラム

Country Status (3)

Country Link
US (1) US8370069B2 (ja)
JP (1) JP4873486B2 (ja)
DE (1) DE102008024282B4 (ja)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5561755B2 (ja) * 2009-04-01 2014-07-30 Necソリューションイノベータ株式会社 Selex法用のプライマーの設計方法、プライマーの製造方法、アプタマーの製造方法、プライマーの設計装置、プライマー設計用コンピュータプログラムおよび記録媒体
JP5419145B2 (ja) * 2009-07-24 2014-02-19 Necソフト株式会社 アプタマー分類装置、アプタマー分類方法、プログラムおよび記録媒体
WO2019079412A2 (en) * 2017-10-21 2019-04-25 Cummins Filtration Ip, Inc. FILTRATION SYSTEM COMPRISING A DOUBLE LIP SEAL
CN113936737B (zh) * 2021-10-14 2023-05-23 温州医科大学附属眼视光医院 一种基于rna基序向量比较rna结构的方法、家族聚类方法、评估变构效应的方法、功能注释的方法、系统和设备

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007134529A (ja) 2005-11-11 2007-05-31 Nidec Tosok Corp ボンディング装置

Also Published As

Publication number Publication date
US8370069B2 (en) 2013-02-05
US20080294351A1 (en) 2008-11-27
JP2008283943A (ja) 2008-11-27
DE102008024282B4 (de) 2011-12-22
DE102008024282A1 (de) 2008-12-11

Similar Documents

Publication Publication Date Title
Ben-Bassat et al. A deep neural network approach for learning intrinsic protein-RNA binding preferences
Zhang et al. iRSpot-ADPM: Identify recombination spots by incorporating the associated dinucleotide product model into Chou’s pseudo components
Kleftogiannis et al. Where we stand, where we are moving: surveying computational techniques for identifying miRNA genes and uncovering their regulatory role
Zhang et al. KD-KLNMF: Identification of lncRNAs subcellular localization with multiple features and nonnegative matrix factorization
Liu iEnhancer-PsedeKNC: Identification of enhancers and their subgroups based on Pseudo degenerate kmer nucleotide composition
JP4873486B2 (ja) 核酸配列の二次構造を予測する方法、核酸配列の二次構造の予測装置及び核酸配列の二次構造の予測プログラム
US8200441B2 (en) Method for identifying nucleotide sequence, method for acquiring secondary structure of nucleic acid molecule, apparatus for identifying nucleotide sequence, apparatus for acquiring secondary structure of nucleic acid molecule, program for identifying nucleotide sequence, and program for acquiring secondary structure of nucleic acid molecule
Lu et al. ERNIE-ac4C: a novel deep learning model for effectively predicting N4-acetylcytidine sites
Jiao et al. Prediction of interface residue based on the features of residue interaction network
Wei et al. Predicting transcription factor binding sites by a multi-modal representation learning method based on cross-attention network
Pashaei et al. A novel method for splice sites prediction using sequence component and hidden Markov model
Meng et al. Protein function prediction based on data fusion and functional interrelationship
Li et al. Understanding sequence conservation with deep learning
Tan et al. Deciphering RNA secondary structure prediction: A probabilistic k-rook matching perspective
Zhang et al. A heuristic cluster-based em algorithm for the planted (l, d) problem
Lihu et al. De novo motif prediction using the fireworks algorithm
Cascitti et al. RNACache: A scalable approach to rapid transcriptomic read mapping using locality sensitive hashing
JP4791737B2 (ja) 機能性核酸配列解析方法
CN111599412A (zh) 基于词向量与卷积神经网络的dna复制起始区域识别方法
Hassanzadeh et al. MotifMark: Finding regulatory motifs in DNA sequences
Cao et al. UFold: Fast and Accurate RNA Secondary Structure Prediction with Deep Learning
Kitaygorodsky et al. Predicting localized affinity of RNA binding proteins to transcripts with convolutional neural networks
Koo et al. ResidualBind: Uncovering Sequence-Structure Preferences of RNA-Binding Proteins with Deep Neural Networks
Pathak et al. EMS1: An elegant algorithm for edit distance based motif search
KR101373254B1 (ko) 도메인 예측 장치, 방법 및 이를 컴퓨터에서 실행하기 위한 프로그램을 기록하는 컴퓨터 판독 가능한 기록매체

Legal Events

Date Code Title Description
A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20090813

RD03 Notification of appointment of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7423

Effective date: 20090910

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20090915

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20091009

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20101022

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20101221

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20111028

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20111116

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20141202

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Ref document number: 4873486

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350