AI創薬：薬効や副作用を予測するデータ駆動型アプローチ

山西芳裕 (yamanishi@bioreg.kyushu-u.ac.jp)

九州大学生体防御医学研究所 科学技術振興機構さきがけ

ケモインフォマティクスや創薬における機械学習の世界初の本

Lodhi, H. and Yamanishi, Y., Chemoinformatics and Advanced Machine Learning Perspectives, IGI Global, 2010.

創薬の現状

費用も時間もかかる。•  平均1000億円•  10年以上

ゼロからの創薬は、ほとんど失敗する。•  成功率は３万分の１•  体内動態が不十分•  毒性などの問題が途中で判明

ドラッグリポジショニング（DR） 薬 再配置する（違う病気に）

n 既承認薬や過去に医薬品開発に失敗した化合物の新しい効能を発見し、別の疾患の治療薬として開発

n ヒトでの安全性や体内動態が確認され、製造法も確立

n 高速・低コスト・低リスク

医薬品開発の流れ

1.  創薬標的の同定 2.  大規模化合物スクリーニング 3.  化学構造 適化 4.  薬物動態試験 5.  非臨床試験 6.  第１相臨床試験（安全性の検証） 7.  第２、３相臨床試験（効果の検証） 8.  申請

従来 （１０〜１７年）

ドラッグリポジショニング （３〜１２年）

1.  創薬標的の同定 2.  既存薬スクリーニング 3.  化合構造 適化 4.  薬物動態試験 5.  非臨床試験 6.  第１相臨床試験（安全性の検証） 7.  第２、３相臨床試験（効果の検証） 8.  申請

スキップ

(Nature Reviews Drug Discovery, 3, 673-683, 2004)

早く、安く、安全な薬を患者さんに届けることができる

実例n シルデナフィル（バイアグラ）

¨ 狭心症治療薬　→　男性機能障害の薬 　　　　　　　　　　　→　肺高血圧症の薬

n ミノキシジル（リアップ、ロゲイン） ¨ 高血圧治療薬　→　発毛薬

より多くの患者さんを救済できる!

問題：これまでは偶然の発見に大きく依存していた

本研究では医薬ビッグデータに基づく ドラッグリポジショニングを提案

オブジェクト データの例

薬物 薬理資料、臨床情報、化学構造、副作用報告、治療標的タンパク質、オフターゲット、薬物応答遺伝子発現情報、既知の効能など

低分子化合物

化学構造、標的タンパク質、生理活性情報など

遺伝子 タンパク質

アミノ酸配列、３次元立体構造、機能モチーフ、パスウェイ、タンパク質間相互作用、分子機能、生理学的役割、病理学的役割など

疾患 臨床情報、レセプト、電子カルテ、病因遺伝子、環境因子、バイオマーカー、合併症情報、患者の遺伝子発現情報、異常パスウェイなど

薬と病気の関連を予測する機械学習法（人工知能の基盤技術）を開発

薬１ 病気1

病気2

病気3

すでに知られている効果

コンピュータ上で予測する新しい効果

薬２

薬３

薬４

セレンディピティから脱却したい！

f (x, y) =wTφ(x, y)

DRのためのAI創薬

多様な疾患の分子レベルでの理解が進んできた

n  病因遺伝子 n  発現異常遺伝子 n  疾患関連パスウェイ n  診断マーカー n  環境因子

共通する 分子的特徴

疾患A 疾患B

遺伝子１

遺伝子2

遺伝子3

提案手法の目的

様々な医薬ビッグデータから薬物・タンパク質・疾患ネットワークを大規模予測する

x1

x2

x3

y1

y2

y3

薬（８千個） 疾患（２千個）

既知の関連エッジ

未知の関連エッジ (本研究で予測する)

タンパク質（３万個）

z1

z2

z3

提案手法の目的

様々な医薬ビッグデータから薬物・タンパク質・疾患ネットワークを大規模予測する

x1

x2

x3

y1

y2

y3

薬（８千個） 疾患（２千個）

既知の関連エッジ

未知の関連エッジ (本研究で予測する)

タンパク質（３万個）

z1

z2

z3

既承認薬でも、薬効に関与する標的タンパク質が不明なものは多い（約６割）

タンパク質 etc.

薬物・化合物の化学構造

etc. ゲノム空間 ケミカル 　空間

ケモゲノミクス法

フェノミクス法

人体へのフェノタイプ

薬理 空間 頭痛、吐き気、気分高揚、血圧の変化、

疾患マーカーの変動、 etc.

Targeted Proteins Drugs Side-effects

Molecular scale Phenotypic scale

Component 1 Component 3Component 2

Targeted proteins ^ŝĚĞͲĞīĞĐƚƐ

Dru

gs

Dru

gs

A

C

B

Targeted Proteins Drugs Side-effects

Molecular scale Phenotypic scale

Component 1 Component 3Component 2

Targeted proteins ^ŝĚĞͲĞīĞĐƚƐ

Dru

gs

Dru

gs

A

C

B

Targeted Proteins Drugs Side-effects

Molecular scale Phenotypic scale

Component 1 Component 3Component 2

Targeted proteins ^ŝĚĞͲĞīĞĐƚƐ

Dru

gs

Dru

gs

A

C

B

Targeted Proteins Drugs Side-effects

Molecular scale Phenotypic scale

Component 1 Component 3Component 2

Targeted proteins ^ŝĚĞͲĞīĞĐƚƐ

Dru

gs

Dru

gs

A

C

B

転写 空間

薬物応答遺伝子発現 トランスクリプトミクス法

ゲノムワイドな薬物・タンパク質相互作用予測

ケモゲノミクスの手法 方針：化学構造が似ている薬は同様のタンパク質に

相互作用すると予測 (Yamanishi et al, Bioinformatics, 2008; Faulon et al., Bioinformatics, 2008; Jacob et al, Bioinformatics, 2008,

Keiser et al, Nature, 2009; Yabuuchi et al, Mol Sys Bio, 2011)

タンパク質のアミノ酸配列・ドメイン・モチーフの類似度

possible chemical substructures 薬の化学構造の類似度

475,692 KCF-S substructures (Kotera et al, BMC Syst. Biol., 2013)

kx (xi,x j ) for i, j =1, 2,...,nx

kz (zi,z j ) for i, j =1, 2,...,nz

2012年度「理論分子生物学」講義予定表 !   ゲノム解析、ポスト・ゲノム解析とバイオインフォマティクス !   配列アライメント、ダイナミックプログラミング法 !   ホモロジー検索、FASTA、BLASTアルゴリズム !   マルチプルアライメント、系統樹解析 !   配列モチーフ !   二次構造予測、膜貫通部位予測、立体構造予測 !   遺伝子の機能アノテーション、比較ゲノム解析 !   ネットワーク解析 !   分子生物学データベース !   演習

!   http://goto.kuicr.kyoto-u.ac.jp/lecture/bioinfo.html

1

配列アライメント

!  配列アライメント（sequence alignment） !  ２つのタンパク質または遺伝子の配列を並べて、進化的な関連があるかどうかを調べること

!  ２つの遺伝子が進化的に関連があるか？ !  異なる生物種間で同じ機能を持つ遺伝子 !  一つの生物種内で類似した機能を持つ遺伝子

真正細菌

古細菌

真核生物

原生生物

植物

菌類

動物

分子レベル（配列レベル）の情報：16S rRNA

生物種の系統関係

a

a1 a2 遺伝子重複

a1 a2 a1 a2

種分岐

ホモログ（Homolog） 進化的な起源を同じくする遺伝子

オーソログ（Ortholog） 種分岐の際に同じ遺伝子だったもの 通常同じ機能を持つ

パラログ（Paralog） 遺伝子重複によってできた類似遺伝子 通常異なる機能を持つ

ゼノログ（Xenolog） 水平移動によって得られた類似遺伝子

a

a1 a1’ 種分岐

種１ 種２

オーソログとパラログ

配列アライメント !  配列アライメント

!  ２つのタンパク質または遺伝子の配列を並べて、ホモログ（相同）かどうかを調べること

!  実際には類似性を調べる !  文字の一致（マッチ）、不一致（ミスマッチ）、挿入、欠失を考慮する

!  アライメントのキーポイントは !  アライメントの種類 !  アライメントの方法・アルゴリズム !  アライメントを評価するためのスコア !  スコアの重要性を評価するための統計的基準

!  グローバルアライメント ! 配列全体を並べる

!  ローカルアライメント ! 局所的によく似た部分を探す

LGPSSKQTGKGW-SRIWDN!| +| ||| |+ |!LN-ITKSAGKGAIMRLGDA!

-------TGKG--------! ||| !-------AGKG--------!

マッチ（+, |） ミスマッチ ギャップ・挿入（-）

配列アライメントの種類

Local sequence alighnment kernel (Saigo et al, Bioinformatics, 2004)

Generalized Jaccard index

Drug space Feature space

Protein space

薬物・タンパク質相互作用の予測

A pairwise model for any drug-protein pair ( ʹx , ʹz ) :

f ( ʹx , ʹz ) = aijk((xi,z j ), ( ʹx , ʹz ))j=1

nz

∑i=1

nx

∑ = aijkx (xi, ʹx )kz (z j, ʹz )j=1

nz

∑i=1

nx

∑

Step 1: Pairwise learning

Interacting pair

Non-interacting pair

Learning a model

Drug similarity

Protein similarity

(Yamanishi et al, Bioinformatics, 2008; ; Yamanishi, Adv Neural Inf Process Syst, 2009; Yamanishi et al, Nucleic. Acid Res., 2014)

薬物・タンパク質相互作用の予測

A pairwise model for any drug-protein pair ( ʹx , ʹz ) :

f ( ʹx , ʹz ) = aijk((xi,z j ), ( ʹx , ʹz ))j=1

nz

∑i=1

nx

∑ = aijkx (xi, ʹx )kz (z j, ʹz )j=1

nz

∑i=1

nx

∑

Step 2: Predicting new interactions

Drug similarity

Protein similarity

Feature space

Interacting pair

Non-interacting pair

Prediction

New pairs

(Yamanishi et al, Bioinformatics, 2008; ; Yamanishi, Adv Neural Inf Process Syst, 2009; Yamanishi et al, Nucleic. Acid Res., 2014)

タンパク質 etc.

薬物・化合物の化学構造

etc. ゲノム空間 ケミカル 　空間

ケモゲノミクス法

フェノミクス法 薬理 空間 頭痛、吐き気、気分高揚、血圧の変化、

疾患マーカーの変動、 etc.

Targeted Proteins Drugs Side-effects

Molecular scale Phenotypic scale

Component 1 Component 3Component 2

Targeted proteins ^ŝĚĞͲĞīĞĐƚƐ

Dru

gs

Dru

gs

A

C

B

Targeted Proteins Drugs Side-effects

Molecular scale Phenotypic scale

Component 1 Component 3Component 2

Targeted proteins ^ŝĚĞͲĞīĞĐƚƐ

Dru

gs

Dru

gs

A

C

B

Targeted Proteins Drugs Side-effects

Molecular scale Phenotypic scale

Component 1 Component 3Component 2

Targeted proteins ^ŝĚĞͲĞīĞĐƚƐ

Dru

gs

Dru

gs

A

C

B

Targeted Proteins Drugs Side-effects

Molecular scale Phenotypic scale

Component 1 Component 3Component 2

Targeted proteins ^ŝĚĞͲĞīĞĐƚƐ

Dru

gs

Dru

gs

A

C

B

転写 空間

薬物応答遺伝子発現 トランスクリプトミクス法

人体へのフェノタイプ

ゲノムワイドな薬物・タンパク質相互作用予測

トランスクリプトミクス法 方針：薬物応答遺伝子発現パターンが似ている薬（化

合物）は同様のタンパク質に相互作用すると予測 (Isker et al, Mol Syst Biol, 2013; Wong et al, Plos Comp Bio, 2013; Hizukuri et al, BMC Med Genomics, 2015; Iwata et al, Sci rep, 2017)

x = (x1, x2,, x22276 )T

各薬物に対して、薬物応答の遺伝子発現比 （コントロールに対する）プロファイル LINCS（16268薬物、77細胞、22276遺伝子）

化合物 の種類

細胞 の種類

drug

gene expression profile

human cell line

高い閾値：類似構造の薬物が多いデータ低い閾値：多様な構造の薬物だけのデータ

ベンチマークデータ（4870個の相互作用）での性能評価：化学構造多様性を考慮

トランスクリプトミクス法の性能は化学構造に依存しない

ベンチマークデータの 化学構造多様性

薬物の適応可能疾患の新規予測

n  8270個の薬物（日本・欧米での既承認薬）に対して、標的タンパク質（酵素、GPCR、イオンチャネルなど）を推定

n  1325個の疾患（がん、神経変性疾患、免疫系疾患、精神疾患など）に対する効能を予測

薬と疾患の196048個の新しい関連を予測 （6531薬と1132疾患）

ケモゲノミクスによる予測の例 n  ピオグリタゾン（２型糖尿病の薬）

¨ 新しい適応疾患: パーキンソン病 ¨ 推定されたタンパク質: MAOB (monoamine oxidase B)

Pioglitazone hydrochloride

Similar compound in the learning set

(Quinn LP, et al. Br J Pharmacol., 2008)

トランスクリプトミクスによる予測例

n フェノチアジン（抗精神病薬） ¨ 新しい適応疾患: 前立腺がん ¨ 推定タンパク質： AR (androgen receptor)

Phenothiazine Similar compound in the learning set

Enzalutamide

(Iwata et al, Sci Rep, 2017)

予測されたARの阻害効果はウェット実験で確認できた

(Iwata et al, Sci Rep, 2017)

遺伝子発現変動から薬物が作用するパスウェイが分かる

Regulated genes GenesIn a pathway i kNot in a pathway r - i l - kTotal r l

n  hypergeometric test

drug

gene expression profile

cell line

down-regulated genes

up-regulated genes

163 biological pathways in KEGG

inactivated pathway

activated pathway

P-value

疾患 A

疾患 B

　薬 (候補化合物)

既知の 治療効果

新しい 治療効果

分子機序の 類似性

疾患の類似性に基づくドラッグ リポジショニング

化合物の毒性・副作用予測

　候補化合物

etc.

インシリコ予測 で代替できる？

人体への毒性（数百種類）

心臓毒性、肝臓毒性、腎臓毒性、発癌性 etc.

Targeted Proteins Drugs Side-effects

Molecular scale Phenotypic scale

Component 1 Component 3Component 2

Targeted proteins ^ŝĚĞͲĞīĞĐƚƐ

Dru

gs

Dru

gsA

C

B

Targeted Proteins Drugs Side-effects

Molecular scale Phenotypic scale

Component 1 Component 3Component 2

Targeted proteins ^ŝĚĞͲĞīĞĐƚƐ

Dru

gs

Dru

gs

A

C

B

Targeted Proteins Drugs Side-effects

Molecular scale Phenotypic scale

Component 1 Component 3Component 2

Targeted proteins ^ŝĚĞͲĞīĞĐƚƐ

Dru

gs

Dru

gs

A

C

B

Targeted Proteins Drugs Side-effects

Molecular scale Phenotypic scale

Component 1 Component 3Component 2

Targeted proteins ^ŝĚĞͲĞīĞĐƚƐ

Dru

gs

Dru

gs

A

C

B

•  医薬品・食品・化粧品開発では、化合物の効能だけでなく、毒性の適切な評価が重要

•  哺乳動物を用いた毒性試験は極めて高価（１件につき数千万円）で、世界的に哺乳動物を用いた毒性試験が禁止される流れ

1339種類の副作用を同時に予測したときの精度を評価Relating drug-protein interaction network with drug side-effects

Table 2. Most frequently appearing enriched pathways.

ID KEGG pathway maps

map04080 Neuroactive ligand-receptor interactionmap04020 Calcium signaling pathwaymap04728 Dopaminergic synapsemap04010 MAPK signaling pathwaymap04260 Cardiac muscle contractionmap04727 GABAergic synapsemap04970 Salivary secretionmap04725 Cholinergic synapsemap00590 Arachidonic acid metabolismmap04270 Vascular smooth muscle contraction

All component-based enriched pathways and associatedenrichment scores are shown in Supplementary Table S2.

outcomes of biological pathway perturbations by drugs targeting theproteins that appear in the component.

4.3 Performance evaluationIt is difficult to evaluate the performance of the feature extractionmethod, because there is little knowledge about true associationbetween targeted proteins and side-effects. However, if theextracted components are biologically meaningful, they shouldcontain some general properties which could be exploited for side-effect prediction. We evaluate the performance of the methodby recovering known drug side-effect profiles from drug proteinbinding profiles, using the extracted canonical components.

In previous literature, chemical structure fingerprints were usedfor predicting side-effect profiles in the framework of OCCA (Atiasand Sharan, 2011) and SCCA (Pauwels et al., 2011). Therefore, wemade a comparison between chemical structure-based approach andtargeted protein-based approach in the framework of both OCCAand SCCA by performing the following 5-fold cross-validation.First, drugs in the gold standard set were split into five subsets ofroughly equal sizes, and each subset was used in turn as a test set.Second, the CCA model was trained on the remaining four sets.Third, the prediction score was computed from the test set, based onthe components extracted from the training set. Finally, the modelwas evaluated for prediction accuracy over the five folds.

We evaluated the performance of the methods by the ROC curve(receiver operating characteristic curve) and the Precision-Recallcurve. The ROC curve is the plot of true positives as a function offalse positives based on various prediction score thresholds, wheretrue positives are correctly predicted side-effects and false positivesare incorrectly predicted side-effects. The Precision-Recall curve isthe plot of ”precision” (positive predictive value) as a function of”recall” (sensitivity) based on various thresholds.

We summarized the performance by the AUC (area under theROC curve) score and the AUPR (area under the Precision-Recallcurve) score. To obtain robust results we repeated the cross-validation experiment five times, and computed the mean and thestandard deviation (S.D.) of the AUC scores over the five repetitions.Sparsity parameters c1, c2 ranged from 0 to 1 by 0.1 increments, andm ranged from 10 to 200 by 10 increments. The best results wereobtained with c1 = 0.1, c2 = 0.1, and with m = 80 componentsin the case of SCCA. The same cross-validation experiments wererepeated for OCCA (with no sparsity constraint), and the best resultswere obtained for m = 20 components.

Table 3. Performance evaluation based on 5-fold cross-validation.

Method AUC ± S.D. AUPR ± S.D.Chemical structure based approach:Random 0.5000 ± 0.0000 0.0556 ± 0.0000OCCA 0.8355 ± 0.0010 0.3753 ± 0.0016SCCA 0.8708 ± 0.0007 0.3766 ± 0.0030Target protein based approach:Random 0.5000 ± 0.0000 0.0556 ± 0.0000OCCA 0.8850 ± 0.0007 0.4067 ± 0.0006SCCA 0.8895 ± 0.0002 0.4103 ± 0.0018

Scores of the proposed method are highlighted in bold.

Table 3 shows the resulting AUC and AUPR scores for the fourdifferent approaches, where the prediction scores for all side-effectswere merged, and a global ROC curve and a global PR curve wereevaluated for each method. This indicates that both SCCA andOCCA produce fairly good results and SCCA is slightly better thanOCCA. It also seems that the targeted protein-based approach worksbetter than the chemical structure-based approach. Results suggestthat the targeted protein information is indeed useful for side-effectprediction.

4.4 Prediction of side-effects for uncharacterized drugsIn the DrugBank database, there are still 730 drugs whose targetprotein information is available, but side-effects are not stored inthe SIDER database. Based on their protein-binding profiles, wepredicted the potential side-effects for these uncharacterized drugsusing the SCCA model, all of the 658 reference drugs being used asa training set. All prediction results can be found in SupplementaryTable S3(A). Complete analysis of all predictions is of course out ofreach, so we focused on the side-effect predictions of highest scores.

Some of the top ranked predicted side-effects involve Cinnarizine(DB00568), an anti-histaminic drug used against motion sickness.This drug binds to the Histaminic H1 receptor, which is believedto explain its indication against vomiting in motion sickness. Itspredicted ”tremor” (cyclical movement of a body part) side-effectwas confirmed by literature (Gimenez-Roldan and Mateo, 1991).Interestingly, Cinnarizine also binds to the voltage-dependentcalcium channel involved in muscle contraction, which mightexplain the ”tremor” side effect. ”Constipation” is also predicted forCinnarizine, as we could find in the adverse effect report 6127929-0 of the Food and Drug Administration (FDA). The predicted”dry mouth” side-effect for Cinnazarine was also confirmed fromliterature (Gordon et al., 2001).

The second ranked predicted side-effect is ”diplopia” (doublevision) for benzocaine (DB01086), a surface anesthetic that acts bypreventing transmission of impulse along nerve fibers. Consistentwith this activity, benzocaine is an inhibitor of voltage-dependentsodium channel. The predicted side-effect was confirmed from theliterature (Horowitz et al., 2005). ”Syncope” was another side-effectpredicted for benzocaine with a high score (Walker et al., 2003).The fourth ranked predicted side-effect is ”tremor” for Bepridil(DB01244), an antihypertensive drug. Tremor is indeed one of themost common side-effects for this drug, reported for 5 percent of allpatients (Williams et al., 2002). ”Tachycardia” was also confirmedfor Promazine (DB00420), an antipsychotic agent (Aronson, 2007).The side-effect ”diplopia” for Nisoldipine (DB00401), a calciumchannel blocker used for the management of hypertension, was alsoconfirmed (O’Keefe and Creamer, 1987).

5

Mizutani et al, Bioinformatics, 28, i522-i528, 2012.

Prediction score:

ynew= β jρ jα jT

j=1

m

∑ xnew, where ρ j is canonical correlation

化学構造情報とタンパク質情報の 融合による副作用予測

Yamanishi et al, J. Chem. Inf. Model., 52(12),3284-3292, 2012.

化学構造情報と標的分子情報の融合で予測精度は向上する

まとめ n 機械学習によるデータ駆動型アプローチで、

医薬品候補化合物の潜在的な標的タンパク質、新規効能、副作用の予測が可能。

n 臓器別の疾患分類でなく、メカニズムに基づく疾患分類。

n 承認薬以外の化合物や、漢方薬・生薬への応用も可能。

n 様々な疾患に対する治療薬候補の発見が可能。

n 現在、多くの製薬企業と共同研究で、社会実装。