本文へスキップ

コンピュータのなかの実験室 Web 版 第 3 章

薬の候補を「探す」

インフルエンザの薬は、ウイルスのタンパク質のくぼみにはまって働きを止めます。くぼみに合う分子を、何億もの候補からコンピュータで探します。

←目次を見る手法:ドッキング

インフルエンザの薬は、どうやって生まれるのでしょうか。インフルエンザの薬にかぎらず、薬は、何十万個もの化合物を作って試した中から、ようやく見つかります。 新しい薬が1つできるまでには、9〜17年かかります。 作って試す前に、当たりそうな分子をコンピュータの中で選べないでしょうか。 この章では、コンピュータが何億個もの分子を試して、薬の候補を絞りこむしくみを紹介します。

MISSION

インフルエンザの「飲み薬」を、コンピュータで設計せよ

第 3〜5 章は、1 つのミッションとして通して読めます。応用化学概論の授業でも、このストーリーをたどります。

合格の条件は2つ

  1. ① ウイルスの「はさみ」のくぼみに、ぴったりはまること(STEP 1・2)
  2. ② 飲んだあと、ウイルスのいる場所まで体の中を届くこと(STEP 3)
  1. STEP 1 敵を知る →
  2. STEP 2 くぼみにはまる分子を探す →
  3. STEP 3 体に届く分子に直す(第 4 章)→
  4. STEP 4 データから規則を学ぶ(第 5 章)→
テオ
まずは敵を知ることだ。ウイルスのどこが弱点か、見てみよう。

🖐 さわってみよう

カギをにぎる、表面の2つの「手」

ウイルスの表面には、2種類のとげが生えています。ヘマグルチニン(HA)は、細胞をつかむ「くっつく手」。ノイラミニダーゼ(NA)は、増えたウイルスを細胞から切りはなす「はさみ」です。

インフルエンザウイルスが増えるしくみ

まず、何を止めればウイルスが広がらなくなるのかを見ておきます。ウイルスの表面には、2種類のタンパク質が突き出ています。1つはヘマグルチニンです。 細胞の表面には糖の鎖がたくさん生えていて、その端にはシアル酸という分子がついています。 ヘマグルチニンは、このシアル酸に結合して、ウイルスを細胞につなぎとめます。

🖐 さわってみよう本の図 3-1

インフルエンザウイルスが増えるサイクル

細胞にくっついて中に入り、コピーを作らせて、外へ出ていくまでの場面がくり返し流れます。

①でヘマグルチニン、⑥でノイラミニダーゼがはたらきます。場面の区切り方は、本の図 3-1 と少しちがいます。

細胞の中で増えたウイルスは、細胞の膜をかぶって外へ出てきます。ところが、新しいウイルスのヘマグルチニンも、細胞の表面のシアル酸に結合してしまいます。 このままでは、ウイルスは細胞から離れられません。そこで、もう1つのタンパク質であるノイラミニダーゼが、糖の鎖からシアル酸を切りはなします。 こうしてウイルスは細胞から離れ、まわりの細胞へ広がっていきます。

ノイラミニダーゼのように、ほかの分子の反応を進めるタンパク質を酵素と呼びます。 酵素の表面には、相手の分子(基質)を受けとめて反応を進めるくぼみがあります。この場所を活性部位と呼びます。 ノイラミニダーゼの活性部位に別の分子が結合してシアル酸が入れなくなれば、新しいウイルスは細胞に結合したままになり、まわりへ広がれません。 このように酵素のはたらきを止めることを阻害といい、止める分子を阻害剤と呼びます。インフルエンザの薬を探すには、ノイラミニダーゼの阻害剤を探せばよいのです。

🖐 さわってみよう本の図 3-3

はさみ(ノイラミニダーゼ)を止める

「薬なし」と「薬あり」を切りかえてみてください。阻害剤がはさみのくぼみにはまると、増えたウイルスが細胞から離れられなくなります。

薬の候補を実験で試す

ある分子が阻害剤になるかどうかは、実験で確かめられます。ノイラミニダーゼと、切られると光る目印をつないだ基質を、試験管で混ぜます。 候補の分子を入れたとき、光が強くなる速さが落ちれば、その分子はノイラミニダーゼを阻害しています。

ただし、この実験をするには、候補の分子を1つずつ作らなければなりません。 日本の製薬会社が2015〜2019年の5年間に、あらゆる病気の薬の候補として合成した化合物は、545,967個です。 このうち、薬として承認されたのは何個だと思いますか。

合成した化合物545,967個のうち、臨床試験に進んだのは約53個、薬として承認されたのは24個。合わせて約2.3万個に1個。
薬ができるまでのじょうご(2015〜2019年、日本の製薬会社。あらゆる病気の薬の合計)。帯の幅は形だけで、数の比ではありません。

答えは24個です。およそ2.3万個を作って、ようやく1個が薬になる計算です。

テオ
形が合うか、手をつなげるか。計算なら、作る前に確かめられる。

活性部位に結合する2つの条件

酵素と基質の関係は、よく鍵と鍵穴にたとえられます。鍵穴に入って錠を回せるのは、決まった形の鍵だけです。 同じように、活性部位の形に合う分子だけが、活性部位に入れます。

ただし、形が合うだけでは足りません。活性部位に入った分子は、まわりの原子と引き合わないと、すぐに外へ出ていってしまいます。 分子が活性部位に結合するには、形が合うことと、引き合うことの2つが必要です。

左は鍵と鍵穴で、形が合えば入ることを表す。右は活性部位に入った分子で、OHが水素結合で、マイナスの部分がプラスの相手と塩橋で引き合っている。
活性部位に結合する2つの条件。鍵と鍵穴(左)は、形が合うことだけを表します。右は活性部位に分子が結合したようすの模式図で、緑の点線は水素結合、青の太い点線は塩橋です。うすく重ねた形は、熱運動で少し動いたときの位置です。

鍵と鍵穴のたとえには、合わない所もあります。金属の鍵穴とちがって、本当の活性部位も分子も、熱運動で少しずつ形を変えながら動いています(第 1 章)。

どの原子がどこにあるかが分からなければ、形が合うかどうかも、どこで引き合うかも分かりません。タンパク質の中で原子が3次元に並んだようすを立体構造と呼びます。 立体構造は、タンパク質の結晶にX線を当てて調べる、X線結晶構造解析という方法で決めます。ノイラミニダーゼの立体構造は1983年に決まり、活性部位の形が分かりました。

活性部位で分子が引き合う主なしくみは、水素結合と静電気力(プラスとマイナスの電荷が引き合う力)です。 ノイラミニダーゼの活性部位には、プラスの電荷をもつアルギニンというアミノ酸が3つ並んでいます。 シアル酸のカルボキシ基(COOH)は、体の中の水ではH+を放して、マイナスの電荷をもつCOO−になっていて、3つのアルギニンと引き合います。 このような、プラスとマイナスのイオンどうしの結びつきを塩橋と呼びます。

活性部位に結合したシアル酸の模式図。シアル酸のCOO−が3つのアルギニン(プラス)と塩橋をつくり、OHやNHがまわりのアミノ酸と水素結合をつくる。
ノイラミニダーゼの活性部位に結合したシアル酸(模式図)。緑の地がタンパク質、淡い水色の所が活性部位です。シアル酸は、線の折れ目と端に C がある書き方(骨格構造式。C と、C についた H を省きます)で描いています。青の太い点線が塩橋、緑の点線が水素結合、白い丸がまわりのアミノ酸です。
テオ
第 2 章では、電子の広がり方で電荷が決まるのを見た。その電荷どうしの引き合いが、ここでは結合の決め手になる。

コンピュータの中で1つ試す:ドッキング

形が合うかどうかと、どこで引き合うかは、原子の位置が分かれば計算できます。 候補の分子を活性部位に置いてみて、結合の強さを点数にする計算をドッキングと呼びます。 分子は、単結合のまわりで原子のまとまりが回るので、形を少しずつ変えられます。コンピュータは、向きや形を変えながら点数を付けていきます。

活性部位の立体構造と候補の分子をコンピュータに渡すと、置く・回す・点数にするを何千〜100万通りもくり返し、いちばん点数の高い置き方と、その点数(4点)を返す。
ドッキングで、コンピュータに渡すもの(左)と、返ってくるもの(右)。活性部位と分子は、考え方を見せるための簡単な絵です。

点数の付け方を、手で1回やってみます。ここでは、考え方が分かるように簡単な決め方を使います。 水素結合が1本できるごとに1点、塩橋が1本できるごとに2点を加えます。原子どうしがぶつかって重なる所が1か所あるごとに、3点を引きます。

活性部位に置いた分子が、2本の水素結合と1本の塩橋をつくり、ぶつかりはない。点数は1+1+2で合計4点。
手で1回、点数を付ける。緑の点線が水素結合、青の太い点線が塩橋です。

この置き方では、水素結合が2本、塩橋が1本でき、ぶつかりはありません。点数は1+1+2=4点です。 同じ分子でも、向きを変えると、引き合う相手が入れかわったり、壁にぶつかったりして、点数が上下します。

実際のドッキングでも、考え方は同じです。ただし、形の重なり、水素結合、静電気力などの項を、原子どうしの距離から計算して足し合わせる、もっとていねいな式を使います。 この点数の式をスコア関数と呼びます。

では、大きすぎる分子、小さすぎる分子、形は合うが引き合わない分子、形も引き合いも合う分子では、どれの点数がいちばん高いでしょうか。 答えは、形も引き合いも合う分子です。大きすぎる分子は活性部位の壁にぶつかり、小さすぎる分子は、水素結合の相手のいる壁まで原子が届きません。 上の簡単な決め方では、形がぴったりでも、引き合う原子をもたない分子の点数は上がりません。

下の動く図は、上から見た2次元に簡単にした絵です。形も手も合う分子を入れたら「回す」で向きを変えてみてください。手がつながらなくなります。

🖐 さわってみよう本の図 3-6・3-13

くぼみに結合する分子はどれ?

4つの候補(大きすぎる・形は合うが手がない・形も手も合う・小さすぎる)のボタンを押すと、その分子がくぼみ(活性部位)に入ります。順に試してみてください。形が合うことと、引き合う相手(水素結合の手)がいることの2つがそろうと、はまり度が高くなります。

はまり度は説明用の目安で、本の図 3-11 の点数の付け方とは別のものです。

テラ
手で1回やってみたから、しくみはもう分かったね! あとはこれを、コンピュータが100万回以上くり返すだけ。

何億個をふるいにかける:バーチャルスクリーニング

ドッキングでは、向きと形を変えながら、たくさんの置き方を試します。 2019年に報告された研究では、細菌の酵素へのドッキングで、分子1個あたり平均で約4,000の向きと、それぞれの向きで約280の形を試しました。 合わせると、分子1個につき100万通りを超える置き方を点数にしています。手で数えた1回を、100万回以上くり返すのです。

この計算を、たくさんの分子でくり返します。候補になる分子を集めておき、1つずつドッキングして、点数の高いものだけを残します。 このコンピュータの中でのふるい分けを、バーチャルスクリーニングと呼びます。

🖐 さわってみよう本の図 3-14

コンピュータで、一気にふるいにかける

「スタート」を押すと、候補の分子をコンピュータが次々に調べて、点数の高いものだけを皿に残します。

数の増え方はイメージです。本物のバーチャルスクリーニングでは、何億個もの分子を調べます。

2019年の研究では、注文すれば作れる1億個を超える分子を集め、細菌の酵素では9,900万個を試しました。 その多くは、まだ作られたことのない分子です。合成するのは、点数が上位の分子だけです。

コンピュータの中で、9900万個の候補の分子の集まりを、1分子に約4,000の向き×約280の形でドッキングし、点数の順に並べる。上位から44個を作って阻害の実験をすると、当たりは5個(11%)だった。
バーチャルスクリーニングのふるい(2019年の例)。

2020年には、10億個を超える分子をふるいにかけた報告もあります。分子1個のドッキングに平均15秒かかるとすると、コア(計算をする部品)1個では475年かかります。 1万個のコアを同時に使えば、約2週間で終わる計算です。

では、点数が上位の分子は、本当に阻害剤になるのでしょうか。細菌の酵素では、上位から選んで作った分子の中から阻害剤が見つかりました。 ただし、上位の分子でも、8〜9割は外れです。計算では当たりに見えたのに、実験では結合しない分子を、偽の当たり(偽陽性)と呼びます。 ドッキングでは、分子の置き方は比較的よく当たりますが、結合がどれだけ強いかや、分子の順位はあまり当たりません。そのため、候補を絞ったら、最後は実験で確かめます。

偽の当たりが出る理由は、モデルが省いたものにあります。ドッキングでは、活性部位をほとんど動かない形とみなし、まわりの水は省くか簡単に表すだけで、結合の強さは点数の式で見積もっています。 どのくらい当たるのかは、第 6 章「計算を信じてよいか」で、ほかの計算とあわせて答え合わせをします。

コンピュータは、活性部位の立体構造をもとに、候補の分子が結合する強さを点数にします。何億個もの分子をふるいにかけ、作って試す価値のある候補を絞りこみます。ただし点数には外れがあるので、最後は実験で確かめます。

計算から生まれた薬

今も使われているインフルエンザの薬のうち2つは、活性部位の立体構造をもとに、コンピュータを使って設計されました。

1つ目は、ザナミビル(商品名リレンザ)です。研究者たちは、活性部位のどこに、どんな原子のまとまりを置くと強く引き合うかを計算しました。 その結果、環にある O–H のうち1つ(下の図の点線の丸)を、プラスの電荷をもつまとまりに置きかえると、結合が強くなると予測されました。 この予測にしたがって、N と H がたくさんついたグアニジノ基をつけたのが、ザナミビルです。

シアル酸に似た分子の環のO–Hを、計算の予測にしたがってグアニジノ基に置きかえると、ザナミビルになる。グアニジノ基は活性部位のマイナスのアミノ酸(グルタミン酸など)と引き合う。
環の O–H を、グアニジノ基に置きかえる(環のまわりを切り抜いた模式図)。線と丸の意味は、シアル酸の図と同じです。

2つ目は、オセルタミビル(商品名タミフル)です。分子の環から伸びた枝を、側鎖と呼びます。ザナミビルではO–Hが3つ並んだ側鎖がある場所に、オセルタミビルは、炭素と水素の多い側鎖をもっています。 この側鎖が、活性部位の中にある疎水性(水になじみにくく、油になじみやすい性質)の小さなすき間に入りこむので、オセルタミビルは活性部位に強く結合します。

ザナミビルとオセルタミビルが活性部位に結合した模式図。どちらも3つのアルギニンと塩橋をつくる。ザナミビルはグアニジノ基がマイナスのアミノ酸と引き合い、オセルタミビルは炭素と水素の多い側鎖が疎水性のすき間に入る。
ザナミビル(左)とオセルタミビル(右)が、ノイラミニダーゼの活性部位に結合したようす(模式図)。線と丸の意味は、シアル酸の図と同じです。黄色の帯が疎水性のすき間です。オセルタミビルは、飲んだ後に体の中で一部が外れた形で描いています。

🖐 さわってみよう本の図 3-18

ザナミビルとオセルタミビルを回す

2つの薬の分子を、同じ向きに並べました。ドラッグすると、2つが同じように回ります。真ん中の環と、まわりの飾りの並び方を比べてみてください。

オセルタミビルは飲む前の形(端がエステル)です。体の中でエステルが外れて、効く形になります(第 4 章)。

どちらの薬も、同じ活性部位に結合して、ノイラミニダーゼを阻害します。ところが、リレンザは粉末を吸いこむ吸入薬で、タミフルはカプセルを飲みこむ経口薬です。 活性部位に強く結合しても、飲んだ薬が体の中で作用する場所まで届かなければ効きません。なぜ薬を体に入れる道筋がちがうのかは、第 4 章で考えます。

ウイルスが増えるとき、タンパク質のアミノ酸が入れかわる(変異する)ことがあります。変異によって薬が効きにくくなることを薬剤耐性と呼びます。 タミフルが効きにくくなったウイルスでは、ノイラミニダーゼのアミノ酸が1つ入れかわり、タミフルの側鎖が入る疎水性のすき間がうまくできなくなっていました。リレンザは変わらず強く阻害しました。

山本研究室の生物物理チームは、北海道大学 人獣共通感染症国際共同研究所と共同で、この薬剤耐性が分子のレベルでどう起きるのかを調べています。 薬とタンパク質の結合の強さを、分子動力学を使ったていねいな計算(自由エネルギー計算)で見積もり、ひとつの変異が結合をどれだけ弱めるのかを確かめます。チームのページでは、変異でくぼみの「手」が失われ、結合が弱くなるしくみのイメージを動く図で試せます。 創薬ミッションの続きは、次の章です。見つけた分子を、体に届く形に直します。

3つの問いで読む

ドッキングのまとめ

何をモデルにしたか(省いたもの)
活性部位はほぼ動かない形に。水は省く。結合の強さは点数の式で見積もる
何をくり返したか
置く・回す・点数にするを、1分子で100万通り以上
どこまで信じられるか
上位でも当たりは1〜2割。偽の当たりが混じる
Web 版の目次
  1. はじめに 計算化学への招待
  2. 序章 コンピュータのなかの実験室
  3. 第 1 章 分子の動きを「見る」(分子動力学)
  4. 第 2 章 電子を「見る」(量子化学)
  5. 第 3 章 薬の候補を「探す」(ドッキング)
  6. 第 4 章 体に届くかを「予測する」(ケモインフォマティクス)
  7. 第 5 章 データから「学ぶ」(機械学習)
  8. 第 6 章 計算を信じてよいか(実験との答え合わせ)
  9. 終章 計算化学が切り拓く未来(材料の設計)