JP2019515651A - 活性化合物によって誘発される特定の効果を選択するための、rna編集に基づくアルゴリズム及びインビトロの方法 - Google Patents

活性化合物によって誘発される特定の効果を選択するための、rna編集に基づくアルゴリズム及びインビトロの方法 Download PDF

Info

Publication number
JP2019515651A
JP2019515651A JP2018548117A JP2018548117A JP2019515651A JP 2019515651 A JP2019515651 A JP 2019515651A JP 2018548117 A JP2018548117 A JP 2018548117A JP 2018548117 A JP2018548117 A JP 2018548117A JP 2019515651 A JP2019515651 A JP 2019515651A
Authority
JP
Japan
Prior art keywords
editing
isoforms
algorithm
risk
rna
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.)
Granted
Application number
JP2018548117A
Other languages
English (en)
Other versions
JP7227006B2 (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.)
Alcediag SAS
Original Assignee
Alcediag SAS
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 Alcediag SAS filed Critical Alcediag SAS
Publication of JP2019515651A publication Critical patent/JP2019515651A/ja
Priority to JP2022084981A priority Critical patent/JP2022134134A/ja
Application granted granted Critical
Publication of JP7227006B2 publication Critical patent/JP7227006B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6876Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
    • C12Q1/6883Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for diseases caused by alterations of genetic material
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6869Methods for sequencing
    • C12Q1/6874Methods for sequencing involving nucleic acid arrays, e.g. sequencing by hybridisation
    • 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 or three-dimensional 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 or three-dimensional molecular structures, e.g. structural or functional relations or structure alignment
    • G16B15/30Drug targeting using structural data; Docking or binding prediction
    • 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
    • G16B20/00ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
    • 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
    • G16B30/00ICT specially adapted for sequence analysis involving nucleotides or amino acids
    • 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
    • G16B50/00ICT programming tools or database systems specially adapted for bioinformatics
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2600/00Oligonucleotides characterized by their use
    • C12Q2600/106Pharmacogenomics, i.e. genetic variability in individual responses to drugs and drug metabolism
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2600/00Oligonucleotides characterized by their use
    • C12Q2600/136Screening for pharmacological compounds
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2600/00Oligonucleotides characterized by their use
    • C12Q2600/142Toxicological screening, e.g. expression profiles which identify toxicity
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2600/00Oligonucleotides characterized by their use
    • C12Q2600/158Expression markers

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Biotechnology (AREA)
  • Biophysics (AREA)
  • General Health & Medical Sciences (AREA)
  • Organic Chemistry (AREA)
  • Analytical Chemistry (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Theoretical Computer Science (AREA)
  • Medical Informatics (AREA)
  • Evolutionary Biology (AREA)
  • Wood Science & Technology (AREA)
  • Zoology (AREA)
  • Genetics & Genomics (AREA)
  • Molecular Biology (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Biochemistry (AREA)
  • General Engineering & Computer Science (AREA)
  • Immunology (AREA)
  • Microbiology (AREA)
  • Pathology (AREA)
  • Medicinal Chemistry (AREA)
  • Pharmacology & Pharmacy (AREA)
  • Databases & Information Systems (AREA)
  • Bioethics (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
  • Investigating Or Analysing Biological Materials (AREA)
  • Enzymes And Modification Thereof (AREA)
  • Saccharide Compounds (AREA)

Abstract

本発明は、薬物又は化合物が患者において特定の効果を誘発する確率をインビトロで予測するためのアルゴリズム、及びRNAのA−to−I編集を示す少なくとも1つの標的を使用する、前記アルゴリズムの使用方法に関する。本発明は、この方法を実施するためのキットにも関する。

Description

本発明は、薬物又は化合物が患者において特定の効果を誘発する確率をインビトロで予測するためのアルゴリズム、及びRNAのA−to−I編集を示す少なくとも1つの標的を使用する、前記アルゴリズムの使用方法に関する。本発明は、この方法を実施するためのキットにも関する。
精神障害は世界中の医療制度においてますます多くを占めるようになっている(1)。それらは西洋社会における一般的な障害であり、5人のうち1人は生涯に少なくとも1回は罹患する。精神障害は、脳回路、神経伝達、及び神経可塑性に影響を与える、分子経路の撹乱によって引き起こされる。最近の研究は、メチル化、アセチル化、及び脱アミノ化等の、DNA及びRNAに対するエピジェネティック修飾の変化が、例えば大うつ病、双極性障害、及び統合失調症に関連することを示している(2、3)。最近の研究はまた、RNA上のアデノシン脱アミノ化(RNAのA−to−I編集)を触媒する編集酵素の重要性を明らかにした。この特異的機序は、高度に保存された神経伝達物質及びシナプス関連因子を本質的にコードする遺伝子の機能を直接調節することが示されている(4〜7)。重要なことに、このRNA編集機構及び同種のADAR酵素(RNAに作用するアデノシンデアミナーゼ(Adenosine Deaminases Acting on RNA))の健康及び疾患における役割が、精神障害に罹患している患者の脳におけるその調節解除の証拠の蓄積により、更に明らかになってきた(8、9)。ADARは、二本鎖pre−mRNAステムループに作用して、アデノシン残基を優先的、特異的に脱アミノ化する。コード配列中に存在する残基の脱アミノ化は、異なる薬理学的特性を有する受容体の変異体(例えば、セロトニン2c受容体、グルタミン酸受容体)を産生するアミノ酸置換をもたらす(10)。
脳におけるセロトニン生態学の異常は、うつ病及び/又は自殺行動の根底にある特徴的な形質であることが提唱されている(11〜13)。自殺犠牲者の死後脳組織を分析することにより、本発明者らは、5−HT2CRの薬理学的特性を大きく損なうことが知られている、セロトニン受容体2C(5HT2cR)pre−mRNAに対するRNA編集活性の明確な変化を観察した(10、14)。興味深いことに、自殺犠牲者のヒト皮質における5−HT2cR mRNA編集プロファイルにおけるこれらの変化は、部分的に、SH−SY5Y細胞で観察されるインターフェロン誘発性変化と重なる。本発明者らは、うつ病/自殺患者に関連する5−HT2cRの「RNA編集のサイン」を特徴づけるための特異的バイオマーカを正確に指摘した。
異なる治療クラスに属するいくつかの薬物が、重度の精神医学的有害作用、特にうつ病及び自殺念慮を誘発する可能性があると報告されている(15〜18)。今日では、このような分子を同定するための認可された試験はなく、食品医薬品局(FDA)は全治療クラスに関する一般的なアラートのみを出すことができる。
このように、薬物又は候補薬物が有害な副作用を誘発するリスクを、高い正診率及び高い識別力で判定することができるインビトロ試験を提供する必要がある。
(WHO) WHO.Preventing suicide:A global imperative.2014. The Official Journal of the International Academy for Suicide Research.2010;14(4):291-310. The Official Journal of the Society for Neuroscience.2002 Dec 15;22(24):10529-32. Molecular psychiatry.2001 Jul;6(4):373-9. eLife.2015;4. Rna.2016 Feb;22(2):290-302. Nature communications.2014;5:4726. Brain research Brain research reviews.1998 May;26(2-3):217-29. Brain research Molecular brain research.2004 Apr 29;124(1):70-8. Molecular psychiatry.2008 Nov;13(11):1001-10. Official Publication of the American College of Neuropsychopharmacology.2001 May;24(5):467-77. European psychiatry:the journal of the Association of European Psychiatrists.2010 Jun;25(5):268-71. Archives of general psychiatry.2000 Aug;57(8):729-38. Translational psychiatry.2016; Minor revision. Lancet.2007 Nov 17;370(9600):1706-13. Psychiatria Danubina.2010 Mar;22(1):79-84. Revista brasileira de psiquiatria.2009 Jun;31(2):145-53. Bmj.2010;341:c5812.
本発明者らは、慎重に選択された細胞株(SH−SY5Y)を用いて薬物誘発性精神医学的副作用を予測する、これまでに設計された革新的インビトロアッセイを検証した。本発明者らは、市場で承認された260種にわたる化合物をスクリーニングし、5−HT2cRの編集の薬物誘発性変化を調べた。化合物は、自殺念慮を誘発する可能性が公知の、(FDA警告ラベル及び/又は多数の症例報告を有する)広範囲の治療クラス(抗うつ薬、抗精神病薬、抗肥満薬、抗ウイルス薬、抗炎症薬、抗真菌薬、抗てんかん薬、気分安定薬など)又は誘発しない(精神医学的副作用が報告されていない)広範囲の治療クラスを選択した。このデータを用いて、高い特異度及び敏感度で「リスクあり」化合物を同定した。
第1の態様において、本発明は、患者において化合物、特に薬物が特定の効果を誘発する確率をインビトロで予測するためのアルゴリズムであって、
a)−pre−mRNAがADAR酵素(RNAに作用するアデノシンデアミナーゼ(Adenosine Deaminases Acting on RNA))の基質であるRNAの、前記ADARの作用により異なるアイソフォーム/又は部位の生成がもたらされるA−to−I編集を示す少なくとも1つの標的を選択する工程、
−前記少なくとも1つの標的及び少なくとも前記ADAR酵素を内因的に発現する少なくとも1つの細胞株を選択する工程、
−前記細胞株の細胞が処理された場合、前記(1又は複数の)標的アイソフォーム/又は(1又は複数の)編集部位の相対的割合を用量依存的に変化させることができる陽性対照化合物を選択する工程、
−前記特定の効果を誘発するリスクスコアを注釈付けされた、一定の薬物又は化合物比で構成された分子の集合を選択する工程、
b)前記細胞株の細胞を、陰性対照及び前記陽性対照と併せて、前記分子の集合の個々の単一分子で処理する工程、
c)前記標的のRNA編集レベルの割合を、その編集アイソフォーム/又は部位のそれぞれについて、及び前記集合の分子それぞれについて得るために、その集合の分子で処理された個々の試料の前記少なくとも1つの標的RNA編集プロファイルを解析する工程、
d)−i)単変量解析統計法により、個々のアイソフォーム/又は編集部位について、分子が前記特定の効果を誘発するリスクを識別するその正診率及びその能力を評価する工程;及び/又は
−ii)多変量解析統計法により、アイソフォーム/又は編集部位の各組み合わせについて、分子が前記特定の効果を誘発するリスクを識別するその正診率及びその能力を評価する工程、及び
−iii)最良の識別能を示す組み合わせを選択する工程、
e)アイソフォーム/又は編集部位の前記選択された組み合わせを用いてアルゴリズムを構築し、このようにして得られた前記アルゴリズムを、薬物、化合物、又は分子が患者において前記特定の効果を誘発する確率を予測するために使用する工程、
を含む方法によって得られるアルゴリズムに関する。
化合物とは、本明細書では、特にヒト若しくは動物患者、又は植物において活性であり得る鉱物、化学的、又は生物学的化合物を示すことが意図される。
本明細書において、用語「患者」には、植物も含まれる。
用語「アルゴリズム」には、統計モデル(Cartモデル等)も含まれる。
好ましい実施形態では、本発明による前記アルゴリズムにおいて、前記特定の1又は複数の効果は副作用であり、好ましくは有害な副作用又は所望の副作用から選択され、好ましくは有害な副作用である。
好ましい実施形態では、RNAのA−to−I編集を示す前記標的は、5−HT2cR、PDE8A(ホスホジエステラーゼ8A)、GRIA2(グルタミン酸受容体2)、GRIA3、GRIA4、GRIK1、GRIK2、GRIN2C、GRM4、GRM6、FLNB(フィラミンB)、5−HT2A、GABRA3(GABAα3)、FLNA、CYFIP2からなる群から選択される。
好ましい実施形態では、前記特定の効果、好ましくは副作用、より好ましくは所望又は有害な副作用は、心臓血管、アレルギー、CNS、特に精神医学、皮膚病学、内分泌学、胃腸病学、血液学、感染症、代謝、神経筋、腫瘍学、炎症性、及び肥満の有害副作用を含む群から選択されている。
より好ましくは、精神医学的有害副作用である。
好ましい実施形態では、本発明のアルゴリズムによる前記細胞株の細胞は、前記標的及び(1又は複数の)ADARを内因的に発現する細胞株由来である。
より好ましくは、前記細胞株は、
−前記標的を内因的に発現し、ヒト皮質で観察されるものと同様のADAR酵素の発現を定常状態で示すことができるヒト又は動物の細胞株、
−神経芽細胞腫細胞株、好ましくはヒト細胞株、
−陰性対照またはビヒクル対照に対して正規化された場合に、陽性対照が少なくとも4、好ましくは少なくとも5又は6の誘発倍率でADAR1a発現を誘発した神経芽細胞腫細胞株、並びに
−ヒトSH−SY5Y細胞株、
からなる群から選択される。
好ましい実施形態では、本発明のアルゴリズムの工程b)において、前記細胞株の細胞は、12時間〜72時間、より好ましくは48時間±4時間の間、被験の分子又は対照で処理され、48時間が最も好ましい。
好ましい実施形態では、本発明によるアルゴリズムにおいて、前記陽性対照は、インターフェロンアルファ、又は100IU/mlでインターフェロンのRNA編集プロファイル曲線を再現することができる化合物である(例えば、図6に示される)。SH−SY5Yヒト神経芽細胞腫細胞株は、5−HT2cR mRNAを内因的に発現し、ヒト皮質インターフェロンアルファにおいて観察されるものと同様のADAR酵素発現を定常状態で示すため、これが使用された。
好ましい実施形態では、本発明によるアルゴリズム又はモデルにおいて、工程c)は、個々の分子及び個々の編集アイソフォーム又は編集部位について、前記標的のRNA編集レベルの平均/中央値の相対的割合を得るために、ビヒクル処理対照細胞と比較した、前記細胞株中の個々のアイソフォーム又は部位についてのRNA編集の基底レベルを決定する工程を含む。
好ましくは、前記ビヒクル処理対照細胞はDMSO処理対照細胞である。
好ましい実施形態では、本発明によるアルゴリズム又はモデルにおいて、前記方法は、化合物、特に薬物が前記特定の1又は複数の効果、好ましくは副作用、好ましくは有害又は所望の副作用から選択される副作用、好ましくは有害副作用を誘発する確率を、リスクなし若しくは低リスク、又は高リスク、好ましくはリスクなし又は高リスクによりインビトロで予測する方法である。
特定の好ましい実施形態では、本発明によるアルゴリズム又はモデルにおいて、前記分子の集合は、副作用、好ましくは有害又は所望の副作用から選択される副作用、好ましくは有害副作用である前記特定の1又は複数の効果の誘発に対して高リスク及び非常に低リスク、好ましくはリスクなしのスコアを注釈付けされた、平衡化された分子比で構成される。
「平衡化された分子比」とは、前記有害な副作用を誘発するリスクなし若しくは低リスク、又は高リスクであることが公知であり、特に、心臓血管、アレルギー、CNS、特に精神医学、皮膚病学、内分泌学、胃腸病学、血液学、感染症、代謝、神経筋、腫瘍学、炎症性、及び肥満の治療のクラスの群から選択される、少なくとも3つ、好ましくは少なくとも4つ又は5つの異なる治療クラスを提示する、前記所望の有害な副作用について十分に注釈付けされた分子の集合を示すことが意図される。
好ましくは、前記少なくとも3、4、5、6、7、又は8つの異なる治療クラスのそれぞれに含まれる分子の数は、前記集合の分子の合計の少なくとも10%に相当する。
より好ましい実施形態では、所望の特定の1又は複数の効果、好ましくは副作用、好ましくは有害又は所望の副作用から選択される副作用、好ましくは有害な副作用のクラスを提示する治療クラスは、集合の全分子の20%超、好ましくは25%、30%、又は35%を含む。
好ましい実施形態では、本発明によるアルゴリズムの工程c)において、前記分子の集合は、好ましくは前記集合の個々の分子について異なる濃度で同時に解析される。
好ましい実施形態では、本発明によるアルゴリズムにおいて、工程1)d)i)は、個々のアイソフォームもしくは部位、またはそれらの組み合わせについて、:
−前記特定の1又は複数の効果、好ましくは副作用、好ましくは有害又は所望の副作用から選択される副作用、好ましくは1又は複数の有害副作用について、少なくとも60%、好ましくは70%、好ましくは80%以上の敏感度(Se%)、並びに少なくとも60%、好ましくは70%、好ましくは80%以上の特異度(Sp%)、の最適閾値;
−真の存在[真陽性/(真陽性+偽陽性)]及び真の不在[真陰性/(真陰性+偽陰性)]の割合を評価するための陽性的中率(PPV、%)及び陰性的中率(NPV、%);
を計算する工程を含み、前記方法は、前記(1又は複数の)アイソフォーム/若しくは(1又は複数の)部位又はそれらの組み合わせの選択の全体性能の決定を可能にする。
好ましい実施形態では、本発明によるアルゴリズムまたはモデルにおいて、工程c)において、RNA編集プロファイルは、以下の工程:
−好ましくは2ステップPCR法を用いて(1又は複数の)標的の((1又は複数の)編集部位を含む)目的の配列断片を選択的にシーケンシングする、NGSライブラリーの調製を含むNGS法(次世代シーケンシング);
−得られた全てのNGSライブラリーのシーケンシング;
並びに任意選択で、
−前記シーケンシングデータのバイオインフォマティクス解析
を含み、
前記バイオインフォマティクス解析が、好ましくは以下の工程:
−配列のプリアライメント処理と品質管理
−参照配列に対するアライメント;及び
−編集レベルのコーリング、
を含む方法によって実行され、標的の編集プロファイルを得る。
好ましい実施形態では、本発明によるアルゴリズムにおいて、工程d)i)及びd)ii)並びに工程e)において、前記アルゴリズムの取得を可能にする前記統計的方法は、
−mROCプログラム、特に、AUC(曲線下面積)ROCを最大化し、それぞれの組み合わせについての方程式が新しい仮想マーカZとして以下:
Z=a1(アイソフォーム1)+a2(アイソフォーム2)+...ai(アイソフォームi)+...an(アイソフォームn)
(式中、a1は計算された係数であり、(アイソフォームi)はアイソフォームの標的の個々のRNA編集レベルの相対的割合である)
のように提供され、使用され得る、線形結合を特定するためのプログラム;
及び/又は、
−単変量及び多変量解析に適用され、異なる(1又は複数の)アイソフォーム/若しくは(1又は複数の)編集部位の値で分子の相対的リスクを推定するロジスティック回帰モデル;
及び/又は、
−(1又は複数の)アイソフォーム/若しくは(1又は複数の)編集部位の組み合わせを評価するために適用されるCART(分類ツリー及び回帰ツリー(Classification And Regression Trees))アプローチ;
及び/又は、
−アイソフォーム/若しくは編集部位の組み合わせを評価し、特に、編集アイソフォーム/若しくは部位の重要度をランク付けし、最良のアイソフォーム/若しくは編集部位を組み合わせて、分子の「相対的リスク」を分類するために適用されるランダムフォレスト(RF)アプローチ、
及び/又は任意選択で、
−以下の群:
−サポートベクターマシン(SVM)アプローチ;
−人工ニューラルネットワーク(ANN)アプローチ;
−ベイジアンネットワークアプローチ;
−wKNN(重み付けk近傍)アプローチ;
−部分最小二乗法−判別分析(PLS−DA);並びに、
−線形及び二次判別分析(LDA/QDA);
から選択する、分子の「相対的リスク」についてアイソフォーム/若しくは編集部位の組み合わせを評価するために適用される多変量解析
を含む方法によって実行される。
好ましい実施形態では、本発明によるアルゴリズムにおいて、
−前記少なくとも1つの標的は5−HT2cRであり、
−前記有害副作用は精神医学的有害副作用であり、
−細胞株はヒトSH−SY5Y神経芽細胞腫細胞株であり、
−陽性対照はインターフェロンアルファであり、
並びに、
−試験薬物による前記精神医学的有害副作用の誘発が低リスクであるか、または高リスクであるかを識別することができる部位の組み合わせは、以下の5−HT2cRの部位:
A、B、C、D、及びE、
好ましくは前記部位の少なくとも3、4、又は5つの組み合わせ、
で構成される群から選択される少なくとも2、3、4、又は5つの単一部位の少なくとも1つの組み合わせを含むか、
−あるいは試験薬物による前記精神医学的有害副作用の誘発が低リスクであるか、又は高リスクであるかを識別することができるアイソフォームの組み合わせは、以下の5−HT2cRのアイソフォーム:
A、B、AB、ABC、AC、C、D、AD、AE、ACD、AEC、ABCD、及びNE
好ましくは前記アイソフォームの少なくとも5、6、又は7つの組み合わせ、
で構成される群から選択される少なくとも2、3、4、5、6、7、8、9、10、11、12、又は13の、単一アイソフォームの少なくとも1つの組み合わせを含み、
並びに任意選択で、
前記アルゴリズム又はモデルの取得を可能にする前記統計的方法は、
−mROCプログラム、ランダムフォレストアプローチ及び/又はCartアルゴリズム
を含む方法により実行される。
第2の態様では、本発明は、薬物、化合物、又は分子が、患者において特定の効果、好ましくは副作用、より好ましくは有害又は所望の副作用を誘発する確率又はリスクを予測するインビトロの方法であって、pre−mRNAがADAR酵素の基質であるRNAの、前記ADARの作用により異なるアイソフォームまたは編集部位の生成がもたらされるA−to−I編集を示す標的を使用し、
A)前記標的のRNA編集レベルの割合を、その編集アイソフォームのそれぞれについて得るために、前記薬物又は化合物又は分子で処理された試料の標的RNA編集プロファイルを解析し、
前記標的RNA編集プロファイルが、前記特定の効果について得られた、請求項1〜15のいずれか1項に記載のアルゴリズム又はモデルにおける分子集合の分子について得られたものとして得られる工程、
B)請求項1〜15のいずれか1項に記載の、前記標的及び前記特定の効果について得られたアルゴリズム又はモデルを用いて、前記薬物若しくは化合物について得られた最終値又は適用されたアルゴリズム若しくはモデルを計算する工程;並びに、
C)工程B)で得られた結果を考慮して、前記薬物又は化合物が、患者において前記特定の効果を誘発するリスクがあるかどうかを、特に低リスク対高リスクで判定する工程、
を含む方法に関する。
別の実施形態では、薬物、化合物、又は分子が、患者において特定の効果を誘発する確率又はリスクを予測する、本発明による前記インビトロの方法は、pre−mRNAがADAR酵素の基質であるRNAの、前記ADARの作用により異なるアイソフォーム又は部位の生成がもたらされるA−to−I編集を示す少なくとも2、3、又は4つの標的の組み合わせを使用し、
A)前記標的のそれぞれについてRNA編集レベルの割合を、その編集アイソフォーム又は部位のそれぞれについて得るために、前記薬物又は化合物または分子で処理された試料の前記標的組み合わせの個々の標的RNA編集プロファイルを解析し、
前記標的RNA編集プロファイルの前記それぞれが、前記特定の効果について得られた、請求項1〜15のいずれか1項に記載のアルゴリズム又はモデルにおける分子集合の分子について得られたものとして得られる工程、
B)請求項1〜15のいずれか1項に記載の、前記標的及び前記特定の効果について得られたアルゴリズム又はモデルを用いて、前記薬物若しくは化合物について得られた最終値又は適用されたアルゴリズム若しくはモデルを計算する工程;並びに、
C)工程B)で得られた結果を考慮して、前記薬物又は化合物が、患者において前記特定の効果を誘発するリスクがあるかどうかを、特にリスクなし、若しくは低リスク対高リスクで判定する工程;
を含む。
別の好ましい実施形態では、pre−mRNAがADAR酵素の基質であるRNAのA−to−I編集を示す少なくとも2、3、又は4つの標的の前記組み合わせ、RNAのA−to−I編集を示す標的は、5−HT2cR、PDE8A(ホスホジエステラーゼ8A)、GRIA2(グルタミン酸受容体2)、GRIA3、GRIA4、GRIK1、GRIK2、GRIN2C、GRM4、GRM6、FLNB(フィラミンB)、5−HT2A、GABRA3(GABAα3)、FLNA、CYFIP2からなる群から選択される標的の組み合わせから選択される。
第3の態様では、本発明は、化合物、好ましくは薬物が、患者において前記特定の1又は複数の効果、好ましくは副作用、好ましくは有害又は所望の副作用から選択される副作用、好ましくは1又は複数の有害副作用を誘発するリスクがあるか、特に低リスク対高リスクで判定するためのキットであって、
1)解析により、患者において前記有害副作用を誘発するリスクが判定される最終値を前記試験薬物について得るために、本発明によるアルゴリズムを使用するための説明書、又は化合物若しくは好ましくは薬物が、患者において前記特定の1又は複数の効果、好ましくは副作用、好ましくは有害若しくは所望の副作用から選択される副作用、好ましくは有害副作用を誘発する確率またはリスクを予測する、本発明による方法を適用するための説明書であり、任意選択でROC曲線又はCart決定木を含む説明書;並びに、
2)1)の前記説明書の前記アルゴリズム又は前記モデルを判定するために使用される分子集合の個々の分子について編集RNAプロファイルを得るために必要な試薬に準じた、前記試験薬物について得られる編集RNAプロファイルを判定するための試薬;
を含むキットに関する。
好ましい実施形態において、前記試薬は、本発明のアルゴリズム又はモデルに2ステップPCRを使用する場合、NGSライブラリー調製のためのこの方法に必要なプライマーのセットを含む。
より好ましい態様において、前記試薬は、前記標的の少なくとも1つについて、又は少なくとも2、3、若しくは4つの標的の組み合わせについて、請求項1から17に記載のRNA編集プロファイルを得るために使用されるオリゴヌクレオチド配列を含む。
より好ましい実施形態では、前記試薬は、NGSライブラリー調製のための2ステップPCRに必要なプライマーセットの1つ又は組み合わせを含み、前記少なくとも1つの標的又は前記標的の組み合わせは、5−HT2cR、PDE8A(ホスホジエステラーゼ8A)、GRIA2(グルタミン酸受容体2)、GRIA3、GRIA4、GRIK1、GRIK2、GRIN2C、GRM4、GRM6、FLNB(フィラミンB)、5−HT2A、GABRA3(GABAα3)、FLNA、CYFIP2からなる群から選択される標的から選択される。
別のより好ましい実施形態では、前記試薬は、以下からなる群から選択されるプライマーセットの1つ又は組み合わせを含む:
Figure 2019515651
以下の実施例、図面、及び凡例は、本発明を実施及び使用することができるようにするために、当業者に完全な説明を提供するために選択された。これらの実施例は、本発明者がその発明として考える範囲を限定する意図のものではなく、これ以下の実験のみが実施されたことを示す意図のものでもない。
本発明の他の特徴及び利点は、本明細書の以下に凡例が示される、実施例及び図面の説明の残りの部分において明らかになると思われる。
インターフェロンアルファ誘発性RNA編集(用量応答)の図である。SH−SY5Yヒト神経芽細胞腫細胞株における(IFNα)5−HT2cR mRNA編集「プロファイル」を示す。FNαによる48時間の処理後の、インターフェロンアルファ(IFNα)の効果の用量応答解析である。5−HT2cR mRNAの相対的割合を、NGSベースのシーケンシングによって解析した。プロファイルは、ビヒクル処理対照細胞における5−HT2cR mRNA編集の相対的割合を、IFNα処理細胞において測定した5−HT2cR mRNA編集の相対的割合から差し引くことによって得た。 インビトロアッセイで試験した全260の化合物の治療クラス分類の円グラフである。中枢神経系(CNS)作用化合物の更なるサブクラス分類を、図のB部分に示す。 選択された分子の試験中に適用された実験構成及びアプローチの概略図である。全260の化合物は、5つの独立した生物学的反復実験において試験している。それぞれの個々の細胞培養プレートを10分子、ビヒクル対照(DMSO)及び100IU/mlのインターフェロンアルファで処理した。5つの独立した生物学的反復実験を試験し、NGSベースのRNA編集定量法によって同一の方法で処理された正確に1620個の試料を作製した。 個々のウェルにおけるADAR1a mRNAの発現の図である。分子で48時間処理したSH−SY5Y細胞におけるADAR1a発現の定量的PCR(qPCR)分析。ADAR1a mRNAの発現レベルは、分子、ビヒクル(DMSO)又はIFNαで48時間処理した後、各試料において定量した。単一の生物学的反復実験(n=1)を示す。予想通り、IFNαで処理した各ウェルは、ADAR1a発現の増加を示した(A〜J)。注目すべきことに、分子165もまた、この分子への曝露後にADAR1a mRNA発現レベルの強い増加を示した。 全てのビヒクル対照及びIFNα処理(100UI/ml)SH−SY5Y細胞の生データ(A)である。全150のビヒクル対照(DMSO)及びIFN処理ウェルの全体解析。(A)この表は、全ての5−HT2cR mRNA編集アイソフォームの全基本的統計的特徴を表示している。全実験中に得られたビヒクル及びIFNα処理条件(n=150)を解析にプールし、5−HT2cRに対するIFN誘発性RNA編集変化の標準測定値を生成した。(B)ヒストグラムは、IFN処理によって最も有意に影響を受けた5HT2cR編集アイソフォームを示す。ビヒクル処理(DMSO)及びIFNα処理ウェルについての、全ての5−HT2cR mRNA編集アイソフォームについて平均、中央値、標準偏差、及び変動係数(CV、パーセンテージとして表される)を示す。 プロファイル曲線−RNA編集曲線IFN100である。5HT2cR mRNA編集プロファイルは、ビヒクル処理対照細胞における5−HT2cR mRNA編集の相対的割合を、IFNα処理細胞において測定した5−HT2cR mRNA編集の相対的割合から差し引くことによって得た。平均及び中央値を示し、エラーバーは平均の標準誤差(sem、n=150)を表す。 個々の分子による48時間の処理後に得られた5HT2cR mRNA編集プロファイルの実例である。例は、4種の「リスクあり」化合物(アリピプラゾール(Aririprazole)、セルトラリン、イソトレチノイン、及びタラナバント)(A)並びに4種の「低リスク」分子(リチウム、ケタミン、オンダンセトロン、及びリバビリン)(B)のセットについて示している。参照のIFN(黒色)を各グラフに示す。平均値を示し、エラーバーは平均の標準誤差(sem、n=5)を表す。 低リスク分子と高リスク分子とを識別するための最も代表的な5HT2cR mRNA編集アイソフォームの診断可能性の実例である。箱ひげ図は、5数要約(最小値(smallest observation)、下側四分位(Q1)、中央値(Q2)、上側四分位(Q3)、および最大値(largest observation))を用いて数値データのグループを図式的に示す便利な方法である。箱ひげ図は、基礎となる統計分布を前提とせずに集団間の差異を表示するのに有用である。p値にはウィルコクソン順位和検定を用いた。記号*はp値≦0.05を示し、**はp値≦0.01を示し、***はp値≦0,001を示す。 図15の13のアイソフォームの群から選択された2つのアイソフォームの組み合わせを分子データセット(n=143、低リスク対高リスク分子)に対して使用した、受信者動作特性(ROC)曲線の実例である。決定則:Z=0.121×ACD−0.142×NE。 図15の13個のアイソフォームの群から選択された3つのアイソフォームの組み合わせを分子データセット(n=143、低リスク対高リスク分子)に対して使用した、受信者動作特性(ROC)曲線の実例である。決定則:Z=−0.1449×C+0.569×AE−0.1548×NE。 図15の13個のアイソフォームの群から選択された4つのアイソフォームの組み合わせを分子データセット(n=143、低リスク対高リスク分子)に対して使用した、受信者動作特性(ROC)曲線の実例である。決定則:Z=0.0235xAB+0.1567×ACD+0.3880×AEC−0.1355×NE。 図15の13個のアイソフォームの群から選択された5つのアイソフォームの組み合わせを分子データセット(n=143、低リスク対高リスク分子)に対して使用した、受信者動作特性(ROC)曲線の実例である。決定則:Z=0.016xAB−0.0563×ABC+0.183×ACD+0.386xAEC−0.1428×NE。 図15の13個のアイソフォームの群から選択された6つのアイソフォームの組み合わせを分子データセット(n=143、低リスク対高リスク分子)に対して使用した、受信者動作特性(ROC)曲線の実例である。決定則:Z=0.0157×AB−0.0557×ABC+0.0187×D+0.1817×ACD+0.3883×AEC−0.1426×NE。 図15の13個のアイソフォームの群から選択された7つのアイソフォームの組み合わせを分子データセット(n=143、低リスク対高リスク分子)に対して使用した、受信者動作特性(ROC)曲線の実例である。決定則:Z=−0.0505×B+0.0224×AB+0.001×D+0.163×ACD+0.389×AEC−0.1402×ABCD−0.1385×NE。 13個のアイソフォームの組み合わせを分子データセット(n=143、低リスク対高リスク分子)に対して使用した、受信者動作特性(ROC)曲線の実例である。決定則:0.2035×A+0.1283×B+0.1979×AB+0.1147×ABC+0.1860×AC+0.04331×C+0.1884×D+0.1259×AD+0.7739×AE+0.4295×ACD+0.4775×AEC−0.0415×ABCD+0.0245×NE。 図14の7つのアイソフォームの組み合わせを分子データセット(n=143、低リスク対高リスク分子)に対して使用した、ランダムフォレスト(RF)アルゴリズムの受信者動作特性(ROC)曲線の実例である。すべてのデータセットのROC曲線を黒線で表し、試験データセットのROC曲線を点線で表す(A)。RFモデルにおけるアイソフォームの重要度(重み)(B)(C)。 RFアプローチによる診断能の例である。図15の13のアイソフォームの組み合わせを分子データセット(n=143、低リスク対高リスク分子)に対して使用した、ランダムフォレスト(RF)アルゴリズムの受信者動作特性(ROC)曲線の実例である。全てのデータセットのROC曲線を黒線で表し、試験データセットのROC曲線を点線で表す(A)。RFモデルにおけるアイソフォームの重要度(重み)(B)(C)。 更なる標的:GRIA2(A)、FLNB(B)、及びPDE8A(C)によって測定されたRNA編集活性の定量を示す図である。全ての場合において、IFN処理は、未編集(NE)mRNAの減少によって示されるように、編集されたアイソフォームの相対的割合の増加を誘発した。 LN18(A)及びLN229(B)神経芽細胞腫細胞株(HTR2C)の図である。LN18細胞(A)及びLN229細胞(B)において、5HT2cR mRNA編集プロファイルを、ビヒクル処理対照細胞における5−HT2cR mRNA編集の相対的割合を、IFNα処理細胞において測定した5−HT2cR mRNA編集の相対的割合から差し引くことによって得た。5HT2cR mRNAの平均mRNA編集プロファイルを示している。 予測y CARTアルゴリズムを示す図である。6つのアイソフォームを分子データセット(n=143、低リスク対高リスク分子)に対して使用するCARTアルゴリズムの代表的な決定ツリー及び診断能の実例である。 患者における特定の効果の誘発が低リスク、またはリスクなしである2つの化合物について得られたRNA編集プロファイルを示す図である。例として、ビヒクル対照処理細胞と比較した、リドカイン(A)及びオンダンセトロン(B)で得られたRNA編集プロファイルを提供する。患者における特定の効果の誘発が高リスクである、レセルピン(C)及びフルオキセチン(D)など2つの化合物についてRNA編集プロファイルを得た。 アリピプラゾール(A)、インターフェロン(IFN)(B)、及びレセルピン(C)によって観察された、HTR2Cに対するRNA編集変化の経時的解析を示す図である。 3つの異なる化合物:クロザピン(A)、セルトラリン(B)、及びケタミン(C)によるSH−SY5Y細胞の処理後のRNA編集プロファイルの用量依存的変化を示す図である。
実施例1:材料及び方法
薬物誘発性精神医学的有害副作用に関するデータベースの作成
正確に10mMでDMSOに溶解した1280種の小分子の集合を含む化学ライブラリーをPrestwick Chemicalsから購入した。ライブラリーに含まれる全ての小分子は、100%認可された薬物(FDA、EMA、及び他の機関)であり、最大の薬物類似性を示し、それらの高い化学的及び薬理学的多様性並びにヒトにおけるそれらの公知のバイオアベイラビリティに関して選択されている。化学ライブラリー(Prestwick Chemicals)の購入時に、各単一分子の標的、治療クラス/効果、特許及びADMETに関する詳細な情報を含む、高度に注釈を付けられたデータベースが提供された。本発明者らは、安全情報や症例報告(FDA Medwatch、EMEAなど)を、定期的に更新するデータベースに問い合わせることで、ヒトに処方された場合のこれらの薬物の自殺及びうつ病に関連した有害副作用について出された報告書を検索した。次に、質問の結果をまとめ、化学物質ライブラリーに含まれる各薬物にリスクスコアを帰属させた。自殺及び/又はうつ病に関連する有害副作用が報告された症例の数、薬物の処方の程度、WHOによる必須医薬品のリストに掲載されている、等の様々なパラメータを考慮に入れて、薬物による有害な精神医学的副作用(うつ病及び/又は自殺関連有害副作用)の誘発可能性のリスクを定量化するために、スコア付けシステムを確立した。本発明者らは、有害な精神医学的副作用を誘発するリスクに関する特定の情報を含む包括的なデータベースを取得した。
細胞培養
SH−SY5Yヒト神経芽細胞腫細胞株は、5−HT2cR mRNAを内因的に発現し、ヒト皮質で観察されるものと同様のADAR1酵素発現の定常状態を示すため、これを使用した(Cavarec et al.2013,Weissmann et al.2016 Translational Psychiatry, Patent TOXADAR)。SH−SY5Yヒト神経芽細胞腫細胞株は、Sigma Aldrichから購入した。細胞を、標準的な条件で、37℃において5%COの加湿雰囲気(humified atmosphere)中で通常通りに培養した。血清中にしばしば存在するセロトニンによる5−HT2cR mRNA発現の脱感作及び下方制御のために、透析胎児ウシ血清(FBS Science Tech 参照番号FB−1280D/500)が非透析より好ましかった(Saucier et al.1998)。実験の経過中、細胞を継代数P8からP22までの間培養した。細胞を12ウェル細胞培養プレートに播種する前に、トリプシン処理された細胞懸濁液をKovaslide(Kova International)(透明なプラスチックで作られた、血球計数グリッドの付いた使い捨て顕微鏡用スライド)のチャンバーに2つの独立した添加を行うことによって細胞数を推定した。両方のチャンバーを2人の実験技術者によって計数し、4つの独立した計数結果の平均を、更に12ウェル細胞培養プレートの細胞数およびプレーティングの計算に使用した。
薬理学的処理及び細胞溶解
受領後、Prestwick化学ライブラリー全体を個別のチューブに移し、コード化し、アリコートに分け、その後の使用まで−80℃で保存した。社内で生成された薬物誘発性の精神医学的有害副作用データベースから、高リスクおよび非常に低リスクのスコアを注釈付けされた、平衡化された薬物比で構成された260の分子を選択した。薬物をコード化し、実験構成全体を通して分子をランダムに処理するよう注意を払った。陰性対照(ビヒクルDMSO)及び陽性対照(インターフェロンアルファ)と共に、各実験において260分子全てを同時に解析した。各12ウェル細胞培養プレートにおいて、試験分子のための10個の空の位置を残して、陰性対照及び陽性対照を添加した。順に、単一の反復実験はそれぞれ、27個の12ウェル培養プレート(参照)で構成された。各被験分子について5つの独立した生物学的反復実験(n=5)を作製するように、実験を全く同じように5回繰り返した。実験の過程にわたって、合計1620個、即ち27(ウェルプレートの数)×12(プレート当たりのウェル数)×5(反復実験数)の試料が作製された。予備実験により、10μMでSH−SY5Y細胞に対して致死的であった7分子を同定することができた。これらの分子については、毒性の低下が検出されるまで濃度を適応させ、低下させた。実験前に、分子及び対照の全ての希釈物を調製し、ラックに配置した。324ウェル(27×12ウェル)全ての細胞密度、形態、生存率、及び汚染は、処理前に顕微鏡によって制御した。更に、各ウェルの画像をCanon EOS700デジタルカメラを用いて撮影した。細胞を分子で処理してから正確に48時間後、各ウェルの画像を、規定されたパラメータを用いてデジタルカメラで撮影した。増殖培地を慎重に除去した後、1%β−メルカプトエタノールを含有するRLT溶解緩衝液(Qiagen)350μlを、細胞の完全な化学的溶解のために添加した。12ウェルプレートを積み重ね、RNA抽出まで冷凍庫に保存した。
全RNA抽出、品質管理、及び逆転写
全RNA抽出は、製造業者のガイドライン(Qiagen)に従って行った。RNeasy Mini Kitは、シリカメンブレンRNeasyスピンカラムを用いて細胞から高品質のRNAを迅速に精製する。全ての細胞溶解物を、完全自動化試料調製QIAcubeを用いて抽出した。抽出は、適切なプロトコルを用いて、1回の操作につき12の試料のバッチ(1つの完全な12ウェルプレート)において標準的な手順を用いて処理した。試料調製及びRNA抽出中、RNAseによるRNA分解を避けるための標準的な注意を払った。全ての抽出RNA試料をlabChipGx(Perkin Elmer)により解析して、全RNAを定量及び定性した。更にQubitによる蛍光ベースの定量化を、LabChipGxデータを検証するために実施した。各個別の試料について、RNA品質スコア(RQSスコア)を決定した(1620試料の平均RQSスコア=9.6/10)。次に、試料を正規化し、精製RNAの逆転写を、最終反応容積20μl中1μgのRNA材料から開始して、Takaraキット(PrimeScript RT Takara 参照番号RR037A)を用いて行った。cDNA合成をPeqstar 96xサーモサイクラーで42℃において15分間行い、反応混合物を更なる使用まで4℃において保存した。
定量的PCR(qPCR)による相対的mRNA発現
cDNA合成後、試料を4℃において保存し、その後、LC480システム(Roche)においてqPCRによりADAR1a mRNA発現の解析を行った。標準曲線法を用いて、qPCRデータを定量した。ADAR1aのmRNA発現は、インターフェロンアルファ処理(IFNα)によって誘発されることが公知である。想定されたように、IFNαで48時間処理された全ての試料は、ADAR1a発現の増加を示し、6〜7の間の遺伝子発現の誘発倍率を示した。更に、レセルピン処理はまた、ADAR1a mRNAレベルを一貫して増加させた。
NGSライブラリーの調製
NGSライブラリー調製のために、2ステップPCR法を使用して、前述の、本発明者らなどによってRNA編集に供されたことが確認された5−HT2cRのエキソンVを選択的にシーケンシングした。検証されたPCRプライマーを使用して、目的の領域をPCRによって増幅した。PCR増幅のために、Q5Hot Start High Fidelity酵素(New England Biolabs)を製造業者のガイドライン(参照番号M0494S)に従って使用した。最適化されたPCRプロトコルを用いて、Peqstar 96xサーモサイクラーでPCR反応を行った。PCR後、全ての試料をLabChipGx(Perkin Elmer)により解析し、PCR産物の量及び質の両方を評価した。アンプリコンの純度を決定し、蛍光ベースのQubit法を用いて定量を行った。品質管理後、磁気ビーズ(Mokascience製のHigh Prep PCR MAGbioシステム)を用いて96種類のPCR反応物(マイクロプレート)を精製した。精製後、Qubitシステムを用いてDNAを定量し、精製収率を計算した。次に、Q5 Hot start High fidelity PCR酵素(New England Biolabs)およびIllumina 96 Indexesキット(Nextera XTインデックスキット;Illumina)を用いたPCR増幅により、試料を個別にインデックス付けした。PCR後、試料をライブラリーにプールし、Magbio PCRクリーンアップシステムを用いて精製した。ライブラリーを変性させ、MiSeqプラットフォームでFASTQのみをシーケンシングするためのIlluminaのガイドラインに従ってシーケンシングカートリッジに負荷した。決定された量の5HT2cRアイソフォームを含むプラスミドのプールを各ライブラリーに含めて、各シーケンシング操作におけるシーケンシングの質及びエラーを制御した。更に、標準RNAプールをライブラリーに組み込んで、実験の過程で異なるシーケンシングフローセル間の可変性を決定した。全1620の試料をシーケンシングするためには、18個のMiSeq ReagentキットV3が必要であった(Illumina)。全てのNGSライブラリーを14pMでシーケンシングし、10%Phix(PhiX Control V3)をスパイクして、ライブラリーの多様性を導入した。
実施例2:シーケンシングデータのバイオインフォマティクス解析
1.Fastq配列のプリアライメント処理および品質管理
シーケンシングデータを、Miseqシーケンサ(Illumina)からfastqファイルとしてダウンロードした。シーケンシング品質を評価するために、FastQCソフトウェアバージョン0.11.5を使用して、各生のfastqファイルの初期品質調査を実行した。アダプター配列の除去、それらのサイズ及び品質スコアによる配列のフィルタリング(全ての短いリード(<50nt)及び平均QC<30のリードを除去)からなる予備処理ステップを行った。次に、配列のアライメントの質を向上させ、改善するために、Illumina NGSデータ用のフレキシブルリードトリミングツールを使用した(トリムモマティック・プログラム バージョン0.35)。予備処理ステップを実行後、クリーニングした各fastqファイルの更なる品質管理を実施し、その後、更なる配列処理を実行した。
2.参照配列に対するアライメント
処理されたリードのアライメントを、bowtie2バージョン2.2.5をエンドツーエンドのsensitiveモードで使用して実行した。アライメントはヒトゲノム配列(UCSC hg38)の最新の注釈に対して行われ、多重アライメント領域のリード、アライメント品質の低いリード(Q<40)又は挿入/欠失(INDEL)を含むリードを更なる解析から取り除いた。ファイルアライメントのフィルタリングはSAMtoolsソフトウェアバージョン1.2を使用して実行した。SAMtoolsソフトウェアバージョン1.2は、ソート、マージ、インデックス作成、及び位置ごとのフォーマットでのアライメント生成を含む、SAMフォーマットでアライメントを扱うための様々なユーティリティを提供する。
3.編集レベルのコーリング
次に、SAMtools mpileupを使用して、複数の試料から同時に得られたアライメント結果のデータをパイルアップした。各ゲノム位置における異なるATGCヌクレオチドの数(「塩基数」)を数えるために社内スクリプトを実行した。従って、各ゲノム位置について、自家製スクリプトは「G」を有するリードのパーセンテージを計算する、[「G」リードの数/(「G」リードの数+「A」リードの数)×100]。「G」リードのパーセンテージ>0.5のゲノム位置「A」参照は、このスクリプトによって自動的に検出され、「A−to−I編集部位」とみなされる。最後の段階は、前述の「A−to−I編集部位」の全ての可能な組み合わせのパーセンテージを計算し、標的の編集プロファイルを得ることであった。
4.標的のベースラインと分子編集プロファイルとの比較
本発明者らは、広範囲の分子セット(n=260)の5HT2cR RNA編集プロファイルを解析した。分子を一緒に比較するために、本発明者らは、最初の工程において、ビヒクル処理(DMSO)対照細胞と比較して、SH−SY5Yヒト神経芽細胞腫細胞株の各アイソフォーム/又は部位について本発明者らの標的のRNA編集の基底レベルを決定した。このために、本発明者らは、実施例により示された150にわたるビヒクル独立実験(反復実験)からの5HT2cRのRNA編集レベルの平均を計算した。第2に、社内のスクリプトは、基準対照(CTRL)に対する分子の各反復実験(n=5)の偏差を自動的に計算した。
最後に、個々の分子及び個々の編集アイソフォーム/又は部位について、標的のRNA編集レベルの平均/中央値の相対的割合を得た。
実施例3:統計解析
全ての統計と数値は、統計オープンソースソフトウェア「R/Bioconductor」(19、20)を使用して計算した。RNA編集の値は、通常平均±平均の標準誤差(SEM)として表される。非パラメトリックウィルコクソン順位和検定とウェルチのt検定を用いて差異解析を行った。複数の試験方法では、個々の編集アイソフォーム(例として:5つの編集部位(A、B、C、E、D)からの5HT2cRの未編集アイソフォーム(Ne)を含む32のRNA編集アイソフォーム))のp値を調整して、偽発見率(FDR)を制御することが重要である。Benjamini及びHochberg(BH)手順(21)を「multtest package」を用いて全ての統計的検定に適用し、調整した0.05未満のp値を統計学的に有意とみなした。編集レベルの相対的割合は正規分布していたため、正規化を適用しなかった。全てのデータ分布は、個々の有意なアイソフォームについての中央値及び棒グラフ(barplot)又は箱ひげ図として示す。有意なアイソフォームによる、SH−SY5Yヒト神経芽細胞腫細胞株における5HT2cRのRNA編集レベルを提示する編集プロファイル曲線も各分子について示す。全ての分子群についてアイソフォーム相関を同定するためにピアソン相関検定を適用した。
5HT2cR編集アイソフォーム診断能は、「高リスク分子」群を検出する能力を表す敏感度及び「リスクなし分子又は低リスク分子」群を検出する能力を表す特異度を特徴とすることができる。診断テストの評価結果は、これらの2つの十分に定義された群を比較する2×2の分割表に要約することができる。カットオフを固定することにより、試験の結果に従って2つの群を陽性又は陰性のいずれかに分類することができた。特定のアイソフォームが与えられた場合、a「高リスク」群の中で陽性の試験結果を有する多数の分子(「真陽性」:TP)、及びb「低リスク」群の中で分子の間で、陽性の試験結果を有する分子(「真陰性」:TN)を同定することができる。同じように、c「高リスク」群の中で陰性の試験結果を有する分子(「偽陽性」:FP)、及びd「低リスク」群の中で陰性の試験結果を有する分子(「偽陰性」:FN)が観察される。敏感度はTP/(TP+FN)と定義され、これを本明細書では「真陽性率」と呼ぶ。特異度は、TN/(TN+FP)として定義され、これを本明細書では「真陰性率」と呼ぶ。
各5HT2cR編集アイソフォームの正診率及びその識別力を、受信者動作特性(ROC)解析を用いて評価した。ROC曲線は、様々な値に対する試験の敏感度(Se)と特異度(Sp)との間の相反関係をグラフで視覚化したものである。
更に、全ての5HT2cR編集アイソフォームを互いに組み合わせて、mROCプログラム[Comput.Methods Programs Biomed.2001; 66:199−207]、ロジスティック回帰(22)、並びに2つの教師付き学習アルゴリズム、CART(23)、及びランダムフォレスト(24)のいくつかのアプローチを用いて、敏感度(sensibility)及び特異度の増加可能性を評価した。
mROCはAUC(曲線下面積)ROC(27)を最大にする線形結合(25、26)を特定するための専用プログラムである。それぞれの組み合わせのための方程式が提供され、以下のように新しい仮想マーカZとして使用できる:
Z=a×アイソフォーム1+b×アイソフォーム2+c×アイソフォーム3
(式中、a、b、cは計算された係数であり、アイソフォーム1、2、3はアイソフォームの標的の個々のRNA編集レベルの相対的割合である)。
2、3、又は4つの標的の組み合わせを互いに組み合わせて、例えばmROCプログラム又はロジスティック回帰のような多変量アプローチを用いて敏感度(sensibility)及び特異度の増加可能性を評価することができる。個々の組み合わせのための方程式を計算することができ、次のように新しい仮想マーカZnとして使用することができる:
Zn=n1×標的1+n2×標的2+n3×標的3、
(式中、n1、n2、n3...は、計算された係数であり、標的1、2、3は、例えば、標的のレベルと相関する値である)。
異なるアイソフォーム又は部位の値における分子の相対的リスクを推定するために、単変量及び多変量解析にロジスティック回帰モデルを適用した。本発明者らは、アイソフォームを連続変数(データ非掲載)及びカテゴリー変数(三分位値をカットポイントとして使用する)の両方として解析した。最後のケースでは、オッズ比(OR)及びその95%信頼区間を計算する。罰則付きロジスティック回帰(LASSO、リッジ又はElastic−Netアプローチ)も連続変数に適用した。これらの方法の場合、パッケージ:Rソフトウェアバージョン3.2.3のglmnetバージョン2.0−3を使用する。
CART(分類ツリー及び回帰ツリー(Classification And Regression Trees))アプローチも、アイソフォームの組み合わせを評価するために適用した。この決定ツリーアプローチは、ユーザにとって容易に理解できる階層グラフによって表される、分類則のセットを生成することを可能にする。ツリーの各ノードで決定が行われる。慣例により、左側の枝は、対象の質問に対する肯定的な応答に対応し、右側の枝は、対象の質問に対する否定的な応答に対応する。次に、分類手順は、「IF−THEN」則のセットとして翻訳することができる(例については図20を参照)。
以前と同様にランダムフォレスト(RF)アプローチを適用し、アイソフォームの組み合わせを評価した。この方法は、分散が制御された決定ツリーの集合を構築するために、Breimanの「バギング」アイデアと機能のランダムな選択を組み合わせる。従って、ランダムフォレストは、編集アイソフォームの重要度をランク付けし、最良のアイソフォームを組み合わせて分子の「相対的リスク」を分類するために使用することができる(図16及び図17を参照)。
CARTとランダムフォレストは教師付き学習方法である。これらの方法は、モデルを構築するために使用されるトレーニングセットと、それを検証するためのテストセットの使用を必要とする。本発明者らのデータセットを分割し、データセットの2/3を学習フェーズに使用し、1/3を検証フェーズに使用する。この分割はランダム化されており、各試料の様々な規制(statutes)の初期割合を尊重している。これらの分類器の予測誤差を推定するために、10倍交差検定法を使用して、オーバーフィット問題を回避するために10回繰り返した。これらのアプローチでは、Rソフトウェアバージョン3.2.3の「rpartパッケージ4.1−10」と「ランダムフォレストパッケージ4.6−12」を使用した。
別の多変量解析を使用して、分子の「相対的リスク」について5HT2cR編集アイソフォームの組み合わせを評価することができる:
−サポートベクターマシン(SVM)アプローチ(28);
−人工ニューラルネットワーク(ANN)アプローチ(29);
−ベイジアンネットワークアプローチ(30);
−wKNN(重み付けk近傍)アプローチ(31);
−部分最小二乗法−判別分析(PLS−DA)(32);
−線形および二次判別分析(LDA/QDA)(33);等。
実施例4:結果
SH−SY5Y細胞株の検証
実験の前に、ヒト神経芽細胞腫細胞株(SH−SY5Y)を漸増用量のインターフェロンで処理し、5HT2cRのRNA編集をNGSベースのアプローチを用いて測定した。予想通り、5HT2cRアイソフォームの相対的割合は変化しており、特に用量依存的に増加可能であり(図1)、この特定の培養細胞株における前述のIFN誘発性応答が確認される。IFNプロファイルは、全く異なる解析法を用いて以前に得られたデータと厳密に一致した(34、35)。
実験手順
一旦細胞株が安定した増殖特性を示し、それに応じてIFN処理に応答した時点で、260分子のスクリーニングを調製した。社内で定められた基準に基づいて、リスクスコアは、化学ライブラリー中の個々の1280分子に帰属させた。実用上の理由から、260分子を選択して、独自のインビトロアッセイで更に試験した。分子の選択手順の間、化学ライブラリー(図2)に含まれる図2で同定された主要な治療クラスのすべて、好ましくは少なくとも3、4、5、6、又は7つの部分を含むように注意を払った。260分子のうち112は抗けいれん剤、抗うつ剤等の中枢神経系障害のための処方薬である(図2B)。全ての分子を移し、処理前に適切な管にアリコートに分けた。260分子のスクリーニングのために選択された実験構成は、10種の分子、ビヒクル対照(DMSO)及び100IU/mlインターフェロンアルファで個別に処理された26個のウェルプレート(12ウェルのプレート)からなり、各細胞培養プレートについて陽性対照及び陰性対照を設けた。追加の細胞培養プレートを使用して、追加の対照ウェルを加えた。各分子を、3週間の間隔で5回の生物学的反復実験において試験した(図3)。厳密に48時間の処理後、細胞を適切な溶解緩衝液で溶解し、更なる処理まで−20℃で保存した。全てのRNA抽出はQiacube自動RNA抽出を用いて行い、プレートは個別に処理した(1回の抽出につき12試料のバッチ)。
相対的ADAR1a mRNA発現
RNA抽出後、cDNAを合成し、384ウェルマイクロプレート中でLC480ライトサイクラー(Roche)においてADAR1a発現を評価した。このようにして、同じバッチの全ての試料を単一のqPCR操作で分析することが可能であった。応答の堅牢性を反映して、12ウェルプレートそれぞれの全てのIFN処理細胞について、ADAR1aのインターフェロン依存性誘発が観察された。興味深いことに、分子165もADAR1a mRNA発現を誘発した(図4A〜4I、プレート17)。この応答は、全ての生物学的反復実験(n=5)で見ることができた。以前にSH−SY5Y細胞で観察されたように、IFNは、ビヒクル対照に対して正規化した場合、6.6の誘発倍率でADAR1a発現を誘発した(表1)。9.31%の変動係数は、生物学的現象の再現性をはっきりと示している。
表1:SH−SY5Y細胞におけるIFN処理後のADAR1a mRNA発現の基本統計学的特徴。平均誘発倍率(DMSO処理対照細胞と比較)標準偏差、中央値、及びCV(パーセンテージとして表される)。
Figure 2019515651
A)5HT2cR編集アイソフォームの単変量解析
SH−SY5Y細胞におけるIFN RNA編集アイソフォームと対照との比較
cDNA合成工程の後、5HT2cRのエキソンVを標的とする2ステップPCRアプローチを適用してNGSライブラリーを構築し、全ての試料中の個々の5HT2cR mRNAの相対的割合を正確に定量した。全てのビヒクル対照及びIFN処理ウェル(n=150)の平均値を図5Aに示し、ヒストグラムとして表す。ビヒクル対照とIFN処理条件との間で、アイソフォームの相対的割合に明確な差異を観察することができる(図5B)。これらのデータは、前述のRNA編集プロファイル(図7A〜7B)を生成するRNA編集プロファイルとして表わされ、前述のプロファイル(図1及びCavarecらを参照されたい)と非常によく一致していた。
例として、SH−SY5Y細胞株における、IFNの存在下での5−HT2cR RNA編集アイソフォーム(n=150)のレベルと、ビヒクル対照(ビヒクル、n=150)のレベルとを比較した場合、5−HT2cRのAC、ABC、AB、A、AE、ACE、D、ABCD、ABE、C、B、BC、及びABCEのRNA編集レベルが有意に変化した(図5A〜5B及び図6)。5−HT2cRの未編集アイソフォーム(Ne)のレベルは、IFN分子とビヒクル対照(基底0)との比較のために最も有意である。更に、本発明者らは、AC、ABC、AB、A、AE、AEC、ABCD、ABE、C、B、BC、及びABECの5−HT2cR RNA編集のレベルの増加、並びに5−HT2cR RNA編集のD及び未編集(Ne)アイソフォームのレベルの減少を観察した。これらの結果は、全体として、SH−SY5Y細胞においてIFNの存在下で5−HT2cRに対するRNA編集活性が増加することを示唆している。
表2:IFN分子(n=150)と対照(n=150)とを比較した場合の5−HT2cR RNA編集レベルの差異解析
Figure 2019515651
SH−SY5Y細胞における高リスク分子対低リスク分子のRNA編集アイソフォームレベルの比較
例えば、低リスク分子(n=82)と高リスク分子(n=61)とを比較すると、5−HT2cRアイソフォームの単一編集又は未編集(Ne)レベルを有意に変化させることができる(図7A〜7B)。5−HT2cRアイソフォームのRNA編集レベル、個々のアイソフォームの曲線下面積(AUC)の受信者動作特性(ROC)解析に基づいて、低リスク又は高リスクの分子の識別が可能であった(表3)。
表3:低リスク分子(n=82)と高リスク分子(n=61)とを比較した場合の、単一編集アイソフォームの識別能
Figure 2019515651
各アイソフォームの正診率及びその識別力を、受信者動作特性(ROC)解析を用いて評価した。ROC曲線は、様々な値に対する試験の敏感度(Se)と特異度(Sp)との間の相反関係をグラフで視覚化したものである。AUCは、その信頼区間(CI)を有する曲線下面積を意味する。ROC曲線は、単一マーカの敏感度(Se%)及び特異度(Sp%)の最適閾値を計算することによる、分子の相対リスクの予測モデルに基づく。単一RNA編集アイソフォームの陽性的中率(PPV、%)及び陰性的中率(NPV、%)を計算して、「自殺副作用群」における高リスク分子の真の存在の割合[真陽性/(真陽性+偽陽性)及び真の不在の割合[真陰性/(真陰性+偽陰性)]を評価した。
B)5−HT2cR編集アイソフォームの多変量解析
低リスク分子と高リスク分子とを比較する場合、mROC(multiple Receiving−Operating−Characteristic)アプローチによる複数マーカ解析が有意にAUCを改善した。以下の組み合わせの13のアイソフォームの群から選択される、例えば2、3、4、5、6、7、又は13のアイソフォームに関連するアイソフォームの組み合わせ:A+B+AB+ABC+AC+C+D+AD+AE+ACD+AEC+ABCD+NE(本発明の方法によって得られた組み合わせ)は、Cavarecら(2013)で得られたものより高い敏感度及び特異度によって報告されるように、高リスク分子において自殺の副作用の高いリスク予測値を有する。本発明の組み合わせで特定された2、3、4、5、6、7、8、9、0、11、12、及び13のアイソフォームを組み合わせた統計解析により、一連の決定則を生成し、図9〜図15及び以下の対応する表4〜9(低リスク分子対高リスク分子)に示すように、各組み合わせについて新しい仮想マーカ(Z)を計算した。
マルチアイソフォームパネルの正診率及びその識別力を、受信者動作特性(ROC)解析を使用して評価した。ROC曲線は、様々な値に対する試験の敏感度(Se)と特異度(Sp)との間の相反関係をグラフで視覚化したものである。AUCは、その信頼区間(CI)を有する曲線下面積を意味する。ROC曲線は、マルチアイソフォームパネルの敏感度(Se%)及び特異度(Sp%)の最適閾値を計算することによる、毒性の高リスクの予測モデルに基づく。組み合わせマーカの陽性的中率(PPV、%)及び陰性的中率(NPV、%)を計算して、自殺/うつ病誘導有害副作用の高リスク分子の真の存在の割合[真陽性/(真陽性+偽陽性)及び真の不在の割合[真陰性/(真陰性+偽陰性)]を評価した。
表4:2つのアイソフォーム(低リスク対高リスク分子)による多変量解析を用いた5−HT2cR編集アイソフォームの性能
Figure 2019515651
決定則:
RD1:Z=0.121×ACD−0.142×NE
表5:3つのアイソフォーム(低リスク対高リスク分子)による多変量解析を用いた5−HT2cR編集アイソフォームの性能
Figure 2019515651
決定則:
RD1:Z=−0.1449×C+0.569×AE−0.1548×NE
表6:4つのアイソフォーム(低リスク対高リスク分子)による多変量解析を用いた5−HT2cR編集アイソフォームの性能
Figure 2019515651
決定則:
RD1:Z=0.0235×AB+0.1567×ACD+0.3880×AEC−0.1355×NE
表7:5つのアイソフォーム(低リスク対高リスク分子)による多変量解析を用いた5−HT2cR編集アイソフォームの性能
Figure 2019515651
決定則:
RD1:Z=0.016×AB−0.0563×ABC+0.183×ACD+0.386×AEC−0.1428×NE
表8:6つのアイソフォーム(低リスク対高リスク分子)による多変量解析を用いた5−HT2cR編集アイソフォームの性能
Figure 2019515651
決定則:
RD1:Z=0.0157×AB−0.0557×ABC+0.0187×D+0.1817×ACD+0.3883×AEC−0.1426×NE
表9:7つのアイソフォーム(低リスク対高リスク分子)による多変量解析を用いた5−HT2cR編集アイソフォームの性能
Figure 2019515651
決定則:
RD1:Z=−0.0505×B+0.0224×AB+0.001×D+0.163×ACD+0.389×AEC−0.1402×ABCD−0.1385×NE
表10:13のアイソフォーム(低リスク対高リスク分子)による多変量解析を用いた5−HT2cR編集アイソフォームの性能
Figure 2019515651
Z=0.2035×A+0.1283×B+0.1979×AB+0.1147×ABC+0.1860×AC+0.04331×C+0.1884×D+0.1259×AD+0.7739×AE+0.4295×ACD+0.4775×AEC−0.0415×ABCD+0.0245×NE。
C)決定ツリーアプローチ:多変量解析
「分類ツリー及び回帰ツリー(Classification And Regression Trees)」を表すCARTアルゴリズムは、決定ツリーアプローチである。これらのツリーは、ユーザにとって容易に理解できる階層グラフとして表される、分類則のセットを構築するのに役立つ。ツリーは、内部ノード(決定ノード)、エッジ、及び末端リーフからなる。これらのノードは、テスト及びこのノードからのエッジのラベルと一致する、該テストに対する応答候補によってラベル付けされる。決定ツリーがバイナリである場合、慣例により、左エッジはテストに対する肯定応答に対応し、右エッジは否定応答に対応する。得られた分類の手順は、決定則に関して即座に翻訳される。
決定ツリーは、教師付き分類の一般的で効率的な方法である。この方法は、モデルを構築するためのトレーニングセットと、それを検証するためのテストセットの使用を必要とする。データセットを構築するにあたっては、「曖昧でない」分子(n=143)のリストを分割し、データセットの90%を学習フェーズ(n=93の薬物)に使用し、10%をテストフェーズ(50の薬物)に使用する。この分割はランダム化されており、各分子の様々な規制(statutes)の初期割合を尊重している。更に、
例として、6個のRNA編集アイソフォームを「IFNプロファイル」に組み合わせ、CART法を使用して意思決定モデルを構築した(図20)。
表11:いくつかのアイソフォーム又はC13の組み合わせから選択された5個のアイソフォーム、X1、X2、X3、X4、及びX5の組み合わせによる、分子のデータセット(低リスク対高リスク分子)に対する、CARTアルゴリズムを用いた5−HT2cR編集アイソフォーム診断能
Figure 2019515651
データテストにおける5−HT2cRの5つのRNA編集アイソフォームを用いたCARTモデルの診断能は、低リスク分子対高リスク分子の識別に関しても非常に興味深いと思われる。
D)ランダムフォレストアプローチ:多変量解析
ランダムフォレストは、教師付き分類の一般的且つ効率的な方法である。この方法は、モデルを構築するためのトレーニングセットと、それを検証するためのテストセットの使用を必要とする。データセットを構築するにあたっては、「曖昧でない」分子(n=143)のリストを分割し、データセットの65%を学習フェーズ(n=93の薬物)に使用し、35%をテストフェーズ(50の薬物)に使用する。この分割はランダム化されており、各分子の様々な規制(statutes)の初期割合を尊重している。更に、IFNによる学習データセットの重み付けを行い、「IFNプロファイル」を有する薬物及び「基底0プロファイル」を有する薬物の分離力を改善した。従って、本発明者らはランダムに採取された12のIFN分子と8個の対照(基底0)を学習セット(n=113)に加えた。
例として、本発明者らは、意思決定のモデルを構築するためのランダムフォレスト(RF)アルゴリズムを用いて、「IFNプロファイル」の7及び13の代表的RNA編集アイソフォーム(C7のRD1、表9、及びC13のRD1、表10を参照)を組み合わせた(RFモデルのパラメータ:mtryStart=1、stepFactor=2、ntree=500、improve=0.01;RFモデルのOut Of Bag(OOB)推定値=0.21)(図16A〜C及び17A〜C))。
表12:分子データセット(低リスク分子対高リスク分子)に対して7個のアイソフォームを用いた、ランダムフォレストアルゴリズムを使用した分割表
Figure 2019515651
表13:分子データセット(低リスク対高リスク分子)に対して13個のアイソフォームを用いた、ランダムフォレストアルゴリズムを使用した5−HT2cR編集アイソフォーム診断能
Figure 2019515651
5−HT2cRの7又は13のRNA編集アイソフォームを用いたRFモデルの診断能は、低リスク分子対高リスク分子の識別に関して非常に興味深く、敏感度、特異度、及び正診率は90%(C7の場合)及び95%(C13の場合)より高く、Cavarecら(2013)に開示されたものよりも有意に高い。
実施例5:標的多様化
SH−SY5Y細胞における5HT2cR mRNA編集を更に補うために、本発明者らは更なるADAR基質(GRIA2、FLNB、PDE8A、GRIK2、及びGABRA3)を解析した。興味深いことに、IFN処理は、研究した3つの標的全てについて、RNA編集アイソフォームの相対的割合を変化させた(図11)。従って、バイオマーカを追加することにより、試験の診断能が更に高まることが予見可能である。
実施例6:様々な標的のNGSベースの解析によって得られた、化合物特異的RNA編集プロファイル
化合物特異的RNA編集プロファイルを、標的GABRA3、GRIA2、GRIK2、及びHTR2CのNGSベースの解析によって得た(図21A〜21B参照)。図21A及び21Bにおいて、ヒストグラムは、ビヒクル対照処理細胞と比較した、表示の化合物で処理したヒトSH−SY5Y神経芽細胞腫細胞株の各特異的部位において定量されたRNA編集レベルの相対的割合を示す。
正の値(%)は、ビヒクル処理細胞と比較した、化合物によって誘発される特異的部位におけるRNA編集の増加を示す。反対に、負の値(%)は、ビヒクル処理細胞と比較した、化合物による処理の結果としての特異的部位におけるRNA編集の減少を示す。
患者における特定の効果の誘発が低リスク又はリスクなしである2つの化合物について、RNA編集プロファイルを得た(図21A、21B参照)。例として、ビヒクル対照処理細胞と比較した、リドカイン(A)及びオンダンセトロン(B)で得られたRNA編集プロファイルを提供する。
レセルピン(図21C参照)及びフルオキセチン(図21D参照)のような、患者における特定の効果の誘発が高リスクである2つの化合物について、RNA編集プロファイルを得た。
実施例7:RNA編集の経時的解析
HTR2Cに対するアリピプラゾール、インターフェロン(IFN)、及びレセルピンによるRNA編集の変化の経時的解析を観察した(図22A〜22C参照)。SH−SY5Y細胞を3個全ての化合物で処理すると、RNA編集プロファイルの時間依存的変化が生じた。これは、未編集HTR2Cのそれぞれの相対的割合が経時的な減少を示すことによって、明確に示されている。興味深いことに、処理によって誘発された変化の特異度は、アリピプラゾール(図22A参照)及びインターフェロン(図22B参照)又はレセルピン(図22C参照)の間で得られたプロファイルの違いによって示される。最も好ましいアルゴリズムを適用して、各化合物のリスクスコアを、研究の各時点において決定した(prob(アルゴリズム))。インターフェロンとレセルピンのリスクスコアは全ての時点で高かったが、アリピプラゾール処理は24時間以降にリスクがあることが確実に特定された(下記表14参照)。
表14:アリピプラゾール、インターフェロン、及びレセルピン処理後の、いくつかの時点におけるリスクスコアのレベル
Figure 2019515651
実施例8:SH−SY5Y細胞を異なる化合物で処理した後の、RNA編集プロファイルの用量依存的変化
クロザピン、セルトラリン、及びケタミンの3種の異なる化合物でSH−SY5Y細胞を処理した後、RNA編集プロファイルの用量依存的変化が得られた(図23A〜23C参照)。
RNA編集プロファイルは、ビヒクル処理したSH−SY5Y細胞と比較した、HTR2C RNA編集のそれぞれの相対的割合を表す。
参考文献
1. (WHO) WHO.Preventing suicide:A global imperative.2014.
2. Labonte B, Turecki G.The epigenetics of suicide:explaining the biological effects of early life environmental adversity.Archives of suicide research :official journal of the International Academy for Suicide Research.2010;14(4):291-310.PubMed PMID:21082447.
3. Gurevich I, Englander MT, Adlersberg M, Siegal NB, Schmauss C.Modulation of serotonin 2C receptor editing by sustained changes in serotonergic neurotransmission.The Journal of neuroscience:the official journal of the Society for Neuroscience.2002 Dec 15;22(24):10529-32.PubMed PMID:12486144.
4. Sodhi MS, Burnet PW, Makoff AJ, Kerwin RW, Harrison PJ.RNA editing of the 5-HT(2C) receptor is reduced in schizophrenia.Molecular psychiatry.2001 Jul;6(4):373-9.PubMed PMID:11443520.
5. Alon S, Garrett SC, Levanon EY, Olson S, Graveley BR, Rosenthal JJ, et al.The majority of transcripts in the squid nervous system are extensively recoded by A-to-I RNA editing.eLife.2015;4.PubMed PMID:25569156.Pubmed Central PMCID:4384741.
6. Khermesh K, D'Erchia AM, Barak M, Annese A, Wachtel C, Levanon EY, et al.Reduced levels of protein recoding by A-to-I RNA editing in Alzheimer′s disease.Rna.2016 Feb;22(2):290-302.PubMed PMID:26655226.Pubmed Central PMCID:4712678.
7. Porath HT, Carmi S, Levanon EY.A genome-wide map of hyper-edited RNA reveals numerous new sites.Nature communications.2014;5:4726.PubMed PMID:25158696.Pubmed Central PMCID:4365171.
8. Seeburg PH, Higuchi M, Sprengel R.RNA editing of brain glutamate receptor channels:mechanism and physiology.Brain research Brain research reviews.1998 May;26(2-3):217-29.PubMed PMID:9651532.
9. Yang W, Wang Q, Kanes SJ, Murray JM, Nishikura K.Altered RNA editing of serotonin 5-HT2C receptor induced by interferon:implications for depression associated with cytokine therapy.Brain research Molecular brain research.2004 Apr 29;124(1):70-8.PubMed PMID:15093687.
10. Dracheva S, Patel N, Woo DA, Marcus SM, Siever LJ, Haroutunian V.Increased serotonin 2C receptor mRNA editing:a possible risk factor for suicide.Molecular psychiatry.2008 Nov;13(11):1001-10.PubMed PMID:17848916.
11. Mann JJ, Brent DA, Arango V.The neurobiology and genetics of suicide and attempted suicide:a focus on the serotonergic system.Neuropsychopharmacology:official publication of the American College of Neuropsychopharmacology.2001 May;24(5):467-77.PubMed PMID:11282247.
12. Mann JJ, Currier DM.Stress, genetics and epigenetic effects on the neurobiology of suicidal behavior and depression.European psychiatry:the journal of the Association of European Psychiatrists.2010 Jun;25(5):268-71.PubMed PMID:20451357.Pubmed Central PMCID:2896004.
13. Mann JJ, Huang YY, Underwood MD, Kassir SA, Oppenheim S, Kelly TM, et al.A serotonin transporter gene promoter polymorphism (5-HTTLPR) and prefrontal cortical binding in major depression and suicide.Archives of general psychiatry.2000 Aug;57(8):729-38.PubMed PMID:10920459.
14. Dinah Weissmann, Laurent Vincent, Mark D.Underwood, Laurent Cavarec, Nicolas Salvetat, Siem van der Laan, et al.REGION SPECIFIC ALTERATIONS OF RNA EDITING OF SEROTONIN 2C RECEPTOR IN CORTEX OF SUICIDES WITH MAJOR DEPRESSION.Translational psychiatry.2016; Minor revision.
15. Christensen R, Kristensen PK, Bartels EM, Bliddal H, Astrup A.Efficacy and safety of the weight-loss drug rimonabant:a meta-analysis of randomised trials.Lancet.2007 Nov 17;370(9600):1706-13.PubMed PMID:18022033.
16. Mihanovic M, Restek-Petrovic B, Bodor D, Molnar S, Oreskovic A, Presecki P.Suicidality and side effects of antidepressants and antipsychotics.Psychiatria Danubina.2010 Mar;22(1):79-84.PubMed PMID:20305596.
17. Moreira FA, Crippa JA.The psychiatric side-effects of rimonabant.Revista brasileira de psiquiatria.2009 Jun;31(2):145-53.PubMed PMID:19578688.
18. Sundstrom A, Alfredsson L, Sjolin-Forsberg G, Gerden B, Bergman U, Jokinen J.Association of suicide attempts with acne and treatment with isotretinoin:retrospective Swedish cohort study.Bmj.2010;341:c5812.PubMed PMID:21071484.Pubmed Central PMCID:2978759.
19. Ge Y, S D, Speed TP.Resampling-based multiple testing for microarray data hypothesis.Sociedad de estadistica e investigacion operativa test.2003;12:1-77.
20. Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, et al.Bioconductor:open software development for computational biology and bioinformatics.Genome biology.2004;5(10):R80.PubMed PMID:15461798.Pubmed Central PMCID:545600.
21. Y B, Y H.Controlling the false discovery rate:a practical and powerful approach to multiple testing.J Roy Statist Soc Ser.1995; B 57, 1:289-300.
22. D.G.K, L.L.K, E.M.M.Applied regression analysis and other multivariate methods.PWS-KENT Publishing Company, Boston.1988.
23. L.Breiman, J.Friedman, Olshen R, Stone.C.CART:Classification and Regression Trees.Wadsworth International.1984.
24. Breiman L.Random Forests.Machine Learning.2001;45 ((1)):5-32.
25. Su JQ, Liu JS.Linear combinations of multiple diagnostic markers.Journal of the American Statistical Association; 1993 (88):1350-5.
26. Wang H.A note on iterative marginal optimization:a simple algorithm for maximum rank correlation estimation.Computational Statistics and Data Analysis 2007 (51):2803-12.
27. Staack A, Badendieck S, Schnorr D, Loening SA, Jung K.Combined determination of plasma MMP2, MMP9, and TIMP1 improves the non-invasive detection of transitional cell carcinoma of the bladder.BMC urology.2006;6:19.PubMed PMID:16901349.Pubmed Central PMCID:1560390.
28. Cortes C, Vapnik.V.Support-Vector Networks.Machine Learning.1995;20.
29. Baxt WG.Application of artificial neural networks to clinical medicine.Lancet.1995 Oct 28;346(8983):1135-8.PubMed PMID:7475607.
30. N.F, al.e.Using Bayesian networks to analyze expression data.J Comput Biol.2000;7((3-4)):601-20.
31. K.H, K.P.S.Weighted k-Nearest-Neighbor Techniques and Ordinal Classification.Discussion Paper 399, SFB 386, Ludwig-Maximilians University Munich 2004.
32. Wold S.PLS-regression:a basic tool of chemometrics.Chemometrics and Intelligent Laboratory Systems.2001;58(108-130).
33. Fisher RA.The Use of Multiple Measurements in Taxonomic Problems.Annals of Eugenics.1936;(2)(179-188).
34. D W.EVALUATION OF THE POTENTIAL RISK OF DRUG INDUCED MOOD DISTURBANCE AND SUICIDE:USE OF A DEDICATED PLATFORM.Patent WO2008/152146.2008.
35. Cavarec L, Vincent L, Le Borgne C, Plusquellec C, Ollivier N, Normandie-Levi P, et al.In vitro screening for drug-induced depression and/or suicidal adverse effects:a new toxicogenomic assay based on CE-SSCP analysis of HTR2C mRNA editing in SH-SY5Y cells.Neurotoxicity research.2013 Jan;23(1):49-62.PubMed PMID:22528247.

Claims (16)

  1. 患者において化合物が特定の効果を誘発する確率をインビトロで予測するためのアルゴリズムであって、前記アルゴリズム又はモデルが、
    a)−pre−mRNAがADAR酵素(RNAに作用するアデノシンデアミナーゼ(Adenosine Deaminases Acting on RNA))の基質であるRNAの、少なくとも1つの編集部位に対する前記ADARの作用により異なるアイソフォーム又は部位の生成がもたらされるA−to−I編集を示す少なくとも1つの標的を選択する工程、
    −前記少なくとも1つの標的及び少なくとも前記ADAR酵素を内因的に発現する少なくとも1つの細胞株を選択する工程、
    −前記細胞株の細胞が処理された場合、前記標的アイソフォーム又は編集部位の相対的割合を用量依存的に変化させることができる陽性対照化合物を選択する工程、
    −前記特定の効果を誘発するリスクスコアを注釈付けされた、一定の化合物比で構成された分子の集合を選択する工程、
    b)前記細胞株の細胞を、陰性対照及び前記陽性対照と併せて、前記分子の集合の個々の単一分子で処理する工程、
    c)前記標的のRNA編集レベルの割合を、その編集アイソフォーム及び/又は部位のそれぞれについて、前記集合の分子それぞれについて得るために、その集合の分子で処理された個々の試料の前記少なくとも1つの標的RNA編集プロファイルを解析する工程、
    d)−i)単変量解析統計法により、個々のアイソフォーム/又は編集部位について、分子が前記特定の効果を誘発するリスクを識別するその正診率及びその能力を評価する工程;及び/又は
    −ii)多変量解析統計法により、アイソフォーム/又は編集部位の各組み合わせについて、分子が前記特定の効果を誘発するリスクを識別するその正診率及びその能力を評価する工程、及び
    −iii)最良の識別能を示す組み合わせを選択する工程、
    e)アイソフォーム/又は編集部位の前記選択された組み合わせを用いてアルゴリズムを構築し、このようにして得られた前記アルゴリズムを、前記化合物が患者において前記特定の効果を誘発する確率を予測するために使用する工程、
    を含む方法によって得られる、アルゴリズム。
  2. 前記効果が、有害な副作用又は所望の副作用、好ましくは有害な副作用から選択される副作用である、請求項1に記載のアルゴリズム。
  3. RNAのA−to−I編集を示す前記標的が、5−HT2cR、PDE8A(ホスホジエステラーゼ8A)、GRIA2(グルタミン酸受容体2)、GRIA3、GRIA4、GRIK1、GRIK2、GRIN2C、GRM4、GRM6、FLNB(フィラミンB)、5−HT2A、GABRA3、FLNA、CYFIP2からなる群から選択される、請求項1又は2に記載のアルゴリズム。
  4. 前記特定の効果が有害な精神医学的副作用である、請求項1〜3のいずれか1項に記載のアルゴリズム。
  5. 前記細胞株が、
    −前記標的及び(1又は複数の)ADARを内因的に発現し、
    −神経芽細胞腫細胞株、好ましくはヒト細胞株、
    −陰性対照又はビヒクル対照に対して正規化された場合に、陽性対照が少なくとも4、5、又は6の誘発倍率でADAR1a発現を誘発した神経芽腫細胞株、並びに
    −ヒトSH−SY5Y細胞株、
    からなる群から選択される、請求項1〜4のいずれか1項に記載のアルゴリズム。
  6. 工程b)において、前記細胞株の細胞が、12時間〜72時間、好ましくは48時間±4時間の間、被験分子又は対照で処理される、請求項1〜5のいずれか1項に記載のアルゴリズム。
  7. 前記陽性対照がインターフェロンアルファである、請求項1〜6のいずれか1項に記載のアルゴリズム。
  8. 工程c)が、個々の分子及び個々の編集アイソフォーム/又は編集部位について、前記標的のRNA編集レベルの平均/中央値の相対的割合を得るために、ビヒクル処理対照細胞と比較した、前記細胞株中の個々のアイソフォーム/又は部位についてのRNA編集の基底レベルを決定する工程を含む、請求項1〜7のいずれか1項に記載のアルゴリズム。
  9. 前記方法が、薬物又は化合物が特定の効果を誘発する確率を、リスクなし、若しくは低リスク、又は高リスクによりインビトロで予測する方法である、請求項1〜8のいずれか1項に記載のアルゴリズム。
  10. 前記分子の集合が、それぞれが前記特定の効果の誘発に対して高リスク及び低リスクのスコアを注釈付けされた、平衡化された治療クラスの分子比で構成される、請求項1〜8のいずれか1項に記載のアルゴリズム。
  11. 工程1)d)i)が、個々のアイソフォーム又はそれらの組み合わせについて:
    −前記特定の効果について、少なくとも60の敏感度(Se%)及び少なくとも60%の特異度(Sp%)の最適閾値;
    −真の存在[真陽性/(真陽性+偽陽性)]及び真の不在[真陰性/(真陰性+偽陰性)]の割合を評価するための陽性的中率(PPV、%)及び陰性的中率(NPV、%)
    を計算する工程を含む、請求項1〜10のいずれか1項に記載のアルゴリズム。
  12. 工程c)において、RNA編集プロファイルが、
    −好ましくは2ステップPCR法を用いて標的の(編集部位を含む)目的の配列断片を選択的にシーケンシングする、NGSライブラリーの調製を含むNGS法(次世代シーケンシング);
    −得られた全てのNGSライブラリーのシーケンシング;
    並びに任意選択で
    −前記シーケンシングデータのバイオインフォマティクス解析
    を含み、
    前記バイオインフォマティクス解析が、好ましくは以下の工程:
    −配列のプリアライメント処理と品質管理
    −参照配列に対するアライメント;及び
    −編集レベルのコーリング
    を含む方法によって実行され、標的の編集プロファイルを得る、請求項1〜11のいずれか1項に記載のアルゴリズム又はモデル。
  13. 工程1)d)i)及び1)d)ii)並びに工程1)e)において、前記アルゴリズム又はモデルの取得を可能にする前記統計的方法が:
    −mROCプログラム、特に、AUC(曲線下面積)ROCを最大化し、それぞれの組み合わせについての方程式が新しい仮想マーカZとして以下:
    Z=a1(アイソフォーム1)+a2(アイソフォーム2)+...ai(アイソフォームi)+...an(アイソフォームn)
    (式中、a1は計算された係数であり、(アイソフォームi)はアイソフォームの標的の個々のRNA編集レベルの相対的割合である)
    のように提供され、使用され得る、線形結合を特定するためのプログラム;及び/又は
    −単変量及び多変量解析に適用され、異なるアイソフォームの値で分子の相対的リスクを推定するロジスティック回帰モデル;及び/又は
    −アイソフォームの組み合わせを評価するために適用されるCART(分類ツリー及び回帰ツリー(Classification And Regression Trees))アプローチ;及び/又は
    −アイソフォームの組み合わせを評価し、特に、編集アイソフォームの重要度をランク付けし、最良のアイソフォームを組み合わせて、分子の「相対的リスク」を分類するために適用されるランダムフォレスト(RF)アプローチ、及び/又は任意選択で、
    −以下の群:
    −サポートベクターマシン(SVM)アプローチ;
    −人工ニューラルネットワーク(ANN)アプローチ;
    −ベイジアンネットワークアプローチ;
    −wKNN(重み付けk近傍)アプローチ;
    −部分最小二乗法−判別分析(PLS−DA);並びに
    −線形及び二次判別分析(LDA/QDA)
    からなる群から選択する、分子の「相対的リスク」についてアイソフォームの組み合わせを評価するために適用される多変量解析、
    を含む方法によって実行される、請求項1〜12のいずれか1項に記載のアルゴリズム。
  14. −前記1つの標的が5−HT2cRであり、
    −前記特定の効果が精神医学的有害副作用であり、
    −細胞株がヒトSH−SY5Y神経芽細胞腫細胞株であり、
    −陽性対照がインターフェロンアルファであり、
    並びに、
    −試験薬物による前記精神医学的有害副作用の誘発が低リスクであるか、又は高リスクであるかを識別することができる部位の組み合わせが、以下の5−HT2cRの部位:
    A、B、C、D、及びE、
    好ましくは前記部位の少なくとも3、4、又は5つの組み合わせ、
    で構成される群から選択される少なくとも2、3、4、又は5つの単一部位の少なくとも1つの組み合わせを含むか、
    −又は試験薬物による前記精神医学的有害副作用の誘発が低リスクであるか、又は高リスクであるかを識別することができるアイソフォームの組み合わせが、以下の5−HT2cRのアイソフォーム:
    A、B、AB、ABC、AC、C、D、AD、AE、ACD、AEC、ABCD、及びNE、
    好ましくは前記アイソフォームの少なくとも5、6、又は7つの組み合わせ
    で構成される群から選択される少なくとも2、3、4、5、6、7、8、9、10、11、12、又は13の、単一アイソフォームの少なくとも1つの組み合わせを含み、並びに任意選択で、
    前記アルゴリズムまたはモデルの取得を可能にする前記統計的方法が、
    −mROCプログラム、ランダムフォレストアプローチ、及び/又はCartアルゴリズム
    を含む方法により実行される、請求項1〜13のいずれか1項に記載のアルゴリズム。
  15. 薬物、化合物、又は分子が、患者において特定の効果、好ましくは副作用、より好ましくは有害又は所望の副作用を誘発する確率又はリスクを予測するインビトロの方法であって、pre−mRNAがADAR酵素の基質であるRNAの、前記ADARの作用により異なるアイソフォームまたは部位の生成がもたらされるA−to−I編集を示す標的を使用し、
    A)前記標的のRNA編集レベルの割合を、その編集アイソフォームのそれぞれについて得るために、前記薬物又は化合物または分子で処理された試料の標的RNA編集プロファイルを解析し、
    前記標的RNA編集プロファイルが、前記特定の効果について得られた、請求項1〜14のいずれか1項に記載のアルゴリズム又はモデルにおける分子集合の分子について得られたものとして得られる工程、
    B)請求項1〜14のいずれか1項に記載の、前記標的及び前記特定の効果について得られたアルゴリズム又はモデルを用いて、前記薬物若しくは化合物について得られた最終値又は適用されたアルゴリズム若しくはモデルを計算する工程;並びに、
    C)工程B)で得られた結果を考慮して、前記薬物又は化合物が、患者において前記特定の効果を誘発するリスクがあるかどうかを、特に低リスク対高リスクで判定する工程:
    を含む、方法。
  16. 薬物が、患者において有害副作用を誘発するリスクがあるか、特に低リスク又はリスクなし対高リスクであるかどうかを判定するためのキットであって、
    1)解析により、患者において前記有害副作用を誘発するリスクが判定される最終値を前記試験薬物について得るために、請求項1〜15のいずれか1項に記載のアルゴリズム又はモデルを使用するための説明書、請求項1〜15のいずれか1項に記載の方法を適用するための説明書であり、任意選択でROC曲線又はCart決定木を含む説明書;並びに、
    2)1)の前記説明書の前記アルゴリズム又は前記モデルを判定するために使用される分子集合の個々の分子について編集RNAプロファイルを得るために必要な試薬に準じた、前記試験薬物について得られる編集RNAプロファイルを判定するための試薬;
    を含む、キット。
JP2018548117A 2016-03-11 2017-03-13 活性化合物によって誘発される特定の効果を選択するための、rna編集に基づくアルゴリズム及びインビトロの方法 Active JP7227006B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2022084981A JP2022134134A (ja) 2016-03-11 2022-05-25 活性化合物によって誘発される特定の効果を選択するための、rna編集に基づくアルゴリズム及びインビトロの方法

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP16000600 2016-03-11
EP16000600 2016-03-11
PCT/IB2017/000417 WO2017153849A1 (en) 2016-03-11 2017-03-13 Algorithm and an in vitro method based on rna editing to select particular effect induced by active compounds

Related Child Applications (1)

Application Number Title Priority Date Filing Date
JP2022084981A Division JP2022134134A (ja) 2016-03-11 2022-05-25 活性化合物によって誘発される特定の効果を選択するための、rna編集に基づくアルゴリズム及びインビトロの方法

Publications (2)

Publication Number Publication Date
JP2019515651A true JP2019515651A (ja) 2019-06-13
JP7227006B2 JP7227006B2 (ja) 2023-02-21

Family

ID=55587999

Family Applications (2)

Application Number Title Priority Date Filing Date
JP2018548117A Active JP7227006B2 (ja) 2016-03-11 2017-03-13 活性化合物によって誘発される特定の効果を選択するための、rna編集に基づくアルゴリズム及びインビトロの方法
JP2022084981A Pending JP2022134134A (ja) 2016-03-11 2022-05-25 活性化合物によって誘発される特定の効果を選択するための、rna編集に基づくアルゴリズム及びインビトロの方法

Family Applications After (1)

Application Number Title Priority Date Filing Date
JP2022084981A Pending JP2022134134A (ja) 2016-03-11 2022-05-25 活性化合物によって誘発される特定の効果を選択するための、rna編集に基づくアルゴリズム及びインビトロの方法

Country Status (9)

Country Link
US (1) US20190249251A1 (ja)
EP (1) EP3426801B1 (ja)
JP (2) JP7227006B2 (ja)
CN (1) CN108779494B (ja)
BR (1) BR112018068128A2 (ja)
CA (1) CA3016657A1 (ja)
ES (1) ES2848949T3 (ja)
IL (1) IL261579B2 (ja)
WO (1) WO2017153849A1 (ja)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
BR112018068128A2 (pt) * 2016-03-11 2019-01-15 Alcediag algoritmo e um método in vitro com base em edição de rna para selecionar o efeito particular induzido por compostos ativos
EP3819386A1 (en) * 2019-11-08 2021-05-12 Alcediag Diagnosing of mood disorders using blood rna editing biomarkers

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010070074A2 (en) * 2008-12-17 2010-06-24 Biocortech Evaluation of the potential risk of drug induced mood disturbance and suicide: use of a dedicated platform
JP2013529465A (ja) * 2010-06-24 2013-07-22 ビオコールテク PDE8AmRNA前駆体の編集プロファイリング:治療効力および/もしくは有効性または潜在的薬物副作用を診断するため、ならびに予測および評価するための、ヒト組織における特異的バイオマーカーとしてのADARs活性の使用
JP2014527036A (ja) * 2011-06-27 2014-10-09 ザ ジャクソン ラボラトリー 癌および自己免疫疾患の処置のための方法および組成物

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3272881B9 (en) 2007-06-13 2023-10-04 Alcediag Peripherical tissue sample containing cells expressing the 5htr2c and/or adars as markers of the alteration of the mechanism of the 5htr2c mrna editing and its applications
WO2010064231A1 (en) * 2008-12-02 2010-06-10 Tel Hashomer Medical Research Infrastructure And Services Ltd. Methods of analyzing a-i rna editing and nucleic acid constructs capable of same
US20150056622A1 (en) * 2013-08-06 2015-02-26 Biocortech Editing profiling of pde8a pre -mrna: use as specific biomarker of adars activities in human tissues to diagnose and to predict and assess therapeutic efficacy and/or efficiency or potential drug side effects
BR112018068128A2 (pt) * 2016-03-11 2019-01-15 Alcediag algoritmo e um método in vitro com base em edição de rna para selecionar o efeito particular induzido por compostos ativos

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010070074A2 (en) * 2008-12-17 2010-06-24 Biocortech Evaluation of the potential risk of drug induced mood disturbance and suicide: use of a dedicated platform
JP2012511923A (ja) * 2008-12-17 2012-05-31 ビオコールテク 気分障害および自殺を誘発する、薬剤の潜在的なリスクの評価:専用プラットフォームの使用
JP2013529465A (ja) * 2010-06-24 2013-07-22 ビオコールテク PDE8AmRNA前駆体の編集プロファイリング:治療効力および/もしくは有効性または潜在的薬物副作用を診断するため、ならびに予測および評価するための、ヒト組織における特異的バイオマーカーとしてのADARs活性の使用
JP2014527036A (ja) * 2011-06-27 2014-10-09 ザ ジャクソン ラボラトリー 癌および自己免疫疾患の処置のための方法および組成物

Also Published As

Publication number Publication date
CN108779494B (zh) 2024-05-14
JP7227006B2 (ja) 2023-02-21
EP3426801B1 (en) 2020-10-14
BR112018068128A2 (pt) 2019-01-15
IL261579A (en) 2018-10-31
CA3016657A1 (en) 2017-09-14
WO2017153849A1 (en) 2017-09-14
JP2022134134A (ja) 2022-09-14
EP3426801A1 (en) 2019-01-16
IL261579B2 (en) 2023-07-01
CN108779494A (zh) 2018-11-09
ES2848949T3 (es) 2021-08-13
US20190249251A1 (en) 2019-08-15
IL261579B1 (en) 2023-03-01

Similar Documents

Publication Publication Date Title
Zhao et al. Misuse of RPKM or TPM normalization when comparing across samples and sequencing protocols
Sheng et al. Multi-perspective quality control of Illumina RNA sequencing data analysis
Robins et al. Genetic control of the human brain proteome
Bousman et al. Preliminary evidence of ubiquitin proteasome system dysregulation in schizophrenia and bipolar disorder: convergent pathway analysis findings from two independent samples
Lea et al. Genetic and environmental perturbations lead to regulatory decoherence
Zhang et al. Combining micro-RNA and protein sequencing to detect robust biomarkers for Graves’ disease and orbitopathy
Lee et al. Principles and methods of in-silico prioritization of non-coding regulatory variants
Aguet et al. Molecular quantitative trait loci
Scheubert et al. Tissue-based Alzheimer gene expression markers–comparison of multiple machine learning approaches and investigation of redundancy in small biomarker sets
JP2022134134A (ja) 活性化合物によって誘発される特定の効果を選択するための、rna編集に基づくアルゴリズム及びインビトロの方法
US20230348980A1 (en) Systems and methods of detecting a risk of alzheimer's disease using a circulating-free mrna profiling assay
Hu et al. Microarray meta-analysis identifies acute lung injury biomarkers in donor lungs that predict development of primary graft failure in recipients
Mo et al. Detection of lncRNA-mRNA interaction modules by integrating eQTL with weighted gene co-expression network analysis
Sokolov et al. Methylation in MAD1L1 is associated with the severity of suicide attempt and phenotypes of depression
Deng et al. Cross-platform analysis of cancer biomarkers: a Bayesian network approach to incorporating mass spectrometry and microarray data
Torkenczy et al. Integrated single-cell analysis reveals treatment-induced epigenetic homogenization
Guo et al. Novel meta-analysis pipeline of heterogeneous high-throughput gene expression datasets reveals dysregulated interactions and pathways in asthma
Ramirez et al. Modeling a gene regulatory network of EMT hybrid states for mouse embryonic skin cells
Tsafnat et al. Gene–disease association with literature based enrichment
Wu et al. Establishment of a discriminant mathematical model for diagnosis of deficiency-cold syndrome using gene expression profiling
Sait Computational Analysis of Autism Spectrum Disorder Biomarkers
Andrayas Epigenetic biomarkers of smoking, inflammation, and social differences
Zhou et al. Xin Ren1, 2· Zhuxiao Feng3· Xiaodong Ma1, 2· Lijuan Huo1, 2· Huiying Zhou1, 2· Ayu Bai1, 2· Shujie Feng4·
Chen et al. RNA Sequencing and Genetic Disease
Seyfried Chloe Robins, Yue Liu, Wen Fan, Duc M. Duong, 2 Jacob Meigs, Nadia V. Harerimana, Ekaterina S. Gerasimov, Eric B. Dammer, 2 David J. Cutler, 3 Thomas G. Beach, 4 Eric M. Reiman, 5 Philip L. De Jager, 7 David A. Bennett, 8 James J. Lah, Aliza P. Wingo, 9, 10 Allan I. Levey

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20200213

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20210129

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20210303

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20210512

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20210902

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20220125

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20220525

C60 Trial request (containing other claim documents, opposition documents)

Free format text: JAPANESE INTERMEDIATE CODE: C60

Effective date: 20220525

C11 Written invitation by the commissioner to file amendments

Free format text: JAPANESE INTERMEDIATE CODE: C11

Effective date: 20220607

A911 Transfer to examiner for re-examination before appeal (zenchi)

Free format text: JAPANESE INTERMEDIATE CODE: A911

Effective date: 20220818

C21 Notice of transfer of a case for reconsideration by examiners before appeal proceedings

Free format text: JAPANESE INTERMEDIATE CODE: C21

Effective date: 20220823

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20221108

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20221111

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: 20230131

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20230209

R150 Certificate of patent or registration of utility model

Ref document number: 7227006

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150