分子モデリングの分野では、ドッキング(英: Docking)は、安定なタンパク質複合体を形成するために互いに結合したときに、ある分子の第2の分子に対する好ましい配向を予測する方法である。好ましい配向の知識を使用すれば、例えばスコアリング関数を使用して、2つの分子間の会合の強さや結合親和性を予測することができる。

タンパク質、ペプチド、核酸、炭水化物、脂質などの生物学的に関連する分子間の関連付けは、シグナル伝達において中心的な役割を果たしている。さらに、相互作用する2つのパートナーの相対的な配向は、生成されるシグナルの種類(例えば、アゴニスト対アンタゴニスト)に影響を与える可能性がある。したがって、ドッキングは、生成されるシグナルの強度と種類の両方を予測するのに有用である。

分子ドッキングは、低分子リガンドの適切なターゲット結合部位への結合コンホメーションを予測できるため、構造に基づいた医薬品設計 (structure-based drug design; SBDD) において最も頻繁に使用される手法の一つである。結合挙動の特性評価は、基本的な生化学的プロセスを解明するだけでなく、薬剤の合理的な設計においても重要な役割を果たしている。

## 問題の定義

分子ドッキングは、「錠前 (lock)」を開ける「鍵 (key)」の正しい相対的な向き(錠前の表面のどこに鍵穴があるか、鍵を挿入した後に鍵をどの方向に回すかなど)を見つけたいという、「鍵と鍵穴説」(lock-and-key) の問題と考えることができる。ここでは、タンパク質を「錠前」、リガンドを「鍵」と考えることができる。分子ドッキングは最適化問題として定義されることがあり、興味のある特定のタンパク質に結合するリガンドの最も当て嵌りのよい(ベストフィットな)配向を記述することになる。しかし、リガンドとタンパク質の両方が柔軟であるため、「錠前と鍵」よりも「手袋の中の手」(hand-in-glove) の例えがより適切である。ドッキングプロセスの間、リガンドとタンパク質は立体配座を調整して全体的な「ベストフィット」を達成する。この種の配座調整により、全体的な結合が生じることを誘導適合と呼ぶ。

分子ドッキング研究では、分子認識プロセスを計算機的にシミュレーションすることに焦点を当てている。それは、タンパク質とリガンドの両方の最適な配座、およびタンパク質とリガンドの相対的な配向を達成し、系全体の自由エネルギーを最小化することを目的としている。

## ドッキングアプローチ

分子ドッキングのコミュニティでは、2つの手法が特に人気がある。1つは、タンパク質とリガンドを相補的な表面として記述するマッチング技術を使用している。2つ目のアプローチは、実際のドッキングプロセスをシミュレーションし、リガンドとタンパク質のペア毎の相互作用エネルギーを計算する。どちらのアプローチにも大きな利点があり、いくつかの制限もある。これらを以下に概説する。

### 形状の相補性

幾何学的マッチング (geometric matching)/形状相補性法 (shape complementarity methods) は、タンパク質とリガンドをドッキング可能にする特徴の集合として記述する。これらの特徴には、分子表面/相補的表面記述子 (complementary surface descriptors) を含んでもよい。この場合、受容体の分子表面は、その溶媒に接触可能な表面積の観点から説明され、リガンドの分子表面は、その相補的な(一致する)表面記述の観点から記述される。2つの表面間の相補性は、標的分子とリガンド分子をドッキングする相補的な姿勢を見つけるのに役立つ可能性のある形状マッチングの説明になる。別のアプローチは、主鎖原子のターンを使用してタンパク質の疎水性の特徴を記述することである。さらに別のアプローチは、フーリエ形状記述子の技術を使用することである。形状相補性に基づくアプローチは一般的に高速で堅牢であり、最近の開発によりこれらの方法でリガンドの柔軟性を調べることができるようになったものの、通常はリガンド/タンパク質の配座の動きや動的変化を正確にモデル化することはできない。形状相補性法は数千個のリガンドを数秒でスキャンし、実際にそれらがタンパク質の活性部位に結合できるかどうかを把握することができ、通常はタンパク質とタンパク質の相互作用にまで拡張可能である。それらはまた、最適な結合を見つけるためにリガンドの幾何学的記述を使用するため、ファーマコフォアに基づくアプローチにも適している。

### シミュレーション

ドッキングプロセスのシミュレーションははるかに複雑である。このアプローチでは、タンパク質とリガンドは物理的な距離だけ隔てられており、リガンドはそのコンフォメーション空間内で一定数の「移動」を行った後、タンパク質の活性部位にその位置を見つけだす。この移動には、並進や回転などの剛体変換だけでなく、ねじれ角回転などのリガンドの構造への内部変化も含まれている。リガンドの配座空間におけるこれらの移動のそれぞれは、系の総エネルギーコストを誘発する。したがって、系の総エネルギーは、すべての移動の後に計算される。

ドッキングシミュレーションの明らかな利点は、リガンドの柔軟性を簡単に取り入れることができることであるのに対し、形状相補性技術では、リガンドの柔軟性を取り入れるには独創的な方法を用いなければならない。また、形状補完技術がより抽象的であるのに対し、シミュレーションはより正確に現実をモデル化することができる。

明らかに、シミュレーションは計算コストが高く、大規模なエネルギー地形 (energy landscape) を探索しなければならない。グリッドベースの技術、最適化手法、コンピュータの高速化により、ドッキングシミュレーションがより現実的なものになった。

## ドッキングの機構

ドッキング・スクリーン (docking screen) を実行するために、最初の要件は、関心のあるタンパク質の構造である。通常、構造は、X線結晶構造解析、NMR分光法、低温電子顕微鏡法(クライオEM)などの生物物理学的手法を用いて決定されているが、ホモロジーモデル構築からも導き出すことができる。このタンパク質構造と潜在的なリガンドのデータベースは、ドッキングプログラムの入力として機能する。ドッキングプログラムが成功するかどうかは、検索アルゴリズムとスコアリング関数の2つの要素に依存する。

### 検索アルゴリズム

理論的には、探索空間は、リガンドと対になるタンパク質のすべての可能な配向 (orientations) とコンホメーション (conformations) から構成されている。しかし、実際には現在の計算資源では、探索空間を網羅的に探索することは不可能である。これには、各分子のすべての可能な歪み(分子は動的であり、コンホメーション状態のアンサンブルの中に存在している)と、与えられた粒度レベルでのタンパク質に対するリガンドのすべての可能な回転および並進の方向を列挙する必要がある。使用されているほとんどのドッキングプログラムは、リガンド(柔軟性のあるリガンド)のコンホメーション空間全体を考慮しており、いくつかは柔軟性のあるタンパク質受容体をモデル化しようとしている。ペアの各「スナップショット」は、ポーズ (pose) と呼ばれている。

さまざまな立体配座の探索戦略が、リガンドと受容体に適用されてきた。これらには以下が含まれる:

- 回転可能な結合に関する系統的または確率的なねじれ角探索
- 分子動力学シミュレーション
- 遺伝的アルゴリズムは、新しい低エネルギーコンホメーションを「進化」させ、各ポーズのスコアが次の反復のための個体を選択するために使用されるフィットネス関数として機能する。

#### リガンドの柔軟性

リガンドのコンホメーションは、受容体の不在時に生成され、その後ドッキングされてもよいし、受容体結合空洞 (receptor binding cavity) の存在時にオン・ザ・フライ (on-the-fly) で生成されてもよいし、フラグメントベースのドッキングを使用して、すべての二面体角の完全な回転の柔軟性を持つコンホメーションが生成されてもよい。力場エネルギー評価は、エネルギー的に合理的なコンホメーションを選択するために最も頻繁に使用されるが、知識ベースの方法も使用されている。

ペプチドは非常に柔軟性が高く、比較的大きな分子であるため、その柔軟性をモデル化することは困難な課題となっている。タンパク質-ペプチドのドッキング時にペプチドの柔軟性を効率的にモデル化するために、多くの手法が開発されてきた。

#### 受容体の柔軟性

計算能力は過去10年間 で飛躍的に向上し、コンピュータ支援薬物設計におけるより洗練された計算集約的な手法を使用できるようになった。しかし、ドッキング方法論における受容体の柔軟性を扱うことはまだ厄介な問題である。この困難さの背後にある主な理由は、この種の計算で考慮しなければならない多数の自由度である。しかし、それを無視すると、いくつかの場合では、結合ポーズの予測面でドッキング結果が乏しくなる可能性がある。

異なるコンホメーションの同じタンパク質について実験的に決定された複数の静的構造は、しばしば、受容体の柔軟性を模倣するために使用される。あるいは、結合空洞を取り囲むアミノ酸側鎖の回転異性体ライブラリを検索して、代替的ではあるがエネルギー的に合理的なタンパク質のコンホメーションを生成してもよい。

### スコアリング関数

ドッキングプログラムは多数の潜在的なリガンドのポーズを生成するが、その中にはタンパク質との衝突のために即座に拒否されるものもある。残りはスコアリング関数を使って評価される。この関数は入力としてポーズを取り、そのポーズが好ましい結合相互作用を表す可能性(likelihood; 尤度; もっともらしさ)を示す数値を返し、あるリガンドを別のリガンドに対して相対的なランク付けを行う。

ほとんどのスコアリング関数は、物理学に基づいた分子力学的な力場であり、結合部位内のポーズのエネルギーを推定する。結合への様々な寄与は、加法方程式として書くことができる:

$\bigtriangleup G_{bind}=\bigtriangleup G_{solvent}+\bigtriangleup G_{conf}+\bigtriangleup G_{int}+\bigtriangleup G_{rot}+\bigtriangleup G_{t/t}+\bigtriangleup G_{vib}$

その構成要素は、溶媒効果、タンパク質とリガンドの立体構造的変化、タンパク質とリガンドの相互作用による自由エネルギー、内部回転、単一複合体を形成するためのリガンドと受容体の会合エネルギー、振動モードの変化による自由エネルギーで構成されている。低い(負の)エネルギーの場合は、系が安定していることを示し、結合が相互作用している可能性がある。

別のアプローチとしては、蛋白質構造データバンク (Protein Data Bank; PDB) のようなタンパク質-リガンド複合体の大規模なデータベースから相互作用の知識ベースの統計的ポテンシャルを導出し、この推定されたポテンシャルに従ってポーズの適合性を評価する方法がある。

タンパク質と高親和性リガンドとの複合体については、X線結晶構造解析による構造が多数存在するが、低親和性リガンドについては比較的少ないため、安定性が低く結晶化が難しい傾向にある。このデータを用いて訓練されたスコアリング関数は、高親和性リガンドを正しくドッキングさせることができるが、結合していないリガンドに対しても、もっともらしいドッキングコンホメーションを与える。これにより、多くの偽陽性ヒット、すなわち、タンパク質に結合すると予測されたリガンドが、in vitroで実際には結合しないと予測される。

偽陽性の数を減らす一つの方法は、一般化ボルン法やポアソン-ボルツマン法のような、(潜在的に)より正確ではあるが計算量の多い手法を用いて、上位スコアのポーズのエネルギーを再計算することである。

## ドッキング評価

サンプリングとスコアリング機能の間の相互依存性は、新規化合物のための説得力のあるポーズや結合親和性を予測する際のドッキング能力に影響する。したがって、その予測能力を決定するためには、ドッキングプロトコルの評価が一般的に必要とされている(実験データが利用可能な場合)。ドッキング評価は、次のようなさまざまな戦略を使用して実行できる:

- ドッキング精度計算。
- ドッキングスコアと実験的応答または濃縮係数 (EF) の決定との相関関係。
- イオン結合部位と活性部位内のイオンとの間の距離。
- 誘導適合モデルの存在

### ドッキング精度

ドッキング精度は、リガンドの正しいポーズを実験的に観測されたものと比較して予測する能力を合理化することによって、ドッキングプログラムの適合性を定量化する一つの指標となる。

### 濃縮係数

ドッキング・スクリーンはまた、結合していないと推定される「デコイ」分子の大規模なデータベースの中から、既知のバインダー(結合剤)のアノテーションされたリガンドを濃縮することによって評価することができる。このようにして、ドッキング・スクリーンの成功は、データベース内の非常に多数のデコイ分子の中から、スクリーンの最上位にある少数の既知の活性化合物を濃縮する能力によって評価される。その性能を評価するために、受信者操作特性 (ROC) 曲線の下の面積が広く使われている。

### 見込み

ドッキング・スクリーン (docking screens) から得られた結果のヒットは、薬理学的検証(例: IC50、親和性または力価の測定)の対象となる。前向き研究のみが、特定の標的に対する技術の適合性に対し、決定的な証拠を構成する。

### ベンチマーク

X線結晶構造解析によって決定された結合様式を再現するドッキングプログラムの可能性は、さまざまなドッキング・ベンチマーク・セットによって評価できる。

低分子については、ドッキングとバーチャル・スクリーニングのためのいくつかのベンチマーク・データセットが存在する。例えば、高品質のタンパク質-リガンドのX線結晶構造で構成されるAstex Diverse Setやバーチャル・スクリーニングの性能を評価するためのDirectory of Useful Decoys (DUD)などがある。

ペプチド結合モードを再現する可能性のあるドッキングプログラムの評価は、「ドッキングとスコアリングの効率評価」(Lessons for Efficiency Assessment of Docking and Scoring (LEADS-PEP)) によって評価することができる。

## アプリケーション

低分子リガンドと酵素タンパク質との間の結合相互作用は、酵素の活性化または阻害をもたらし得る。タンパク質が受容体である場合、リガンドの結合は、アゴニズムまたはアンタゴニズムをもたらす可能性がある。ドッキングは医薬品設計の分野で最も一般的に使用されている。ほとんどの薬物は小さな有機分子であり、ドッキングは次のような場合に適用される。

- ヒット同定 - スコアリング関数と組み合わせたドッキングを使用して、潜在的な薬物の大規模なデータベースをin silicoですばやくスクリーニングし、目的のタンパク質ターゲットに結合する可能性のある分子を同定できる(バーチャル・スクリーニングを参照)。
- リード最適化 (lead optimization) - ドッキングは、リガンドがタンパク質に結合する場所と相対的な配向 (結合モードやポーズとも呼ばれる) を予測するために使用できる。この情報を使用して、より強力で選択的なアナログを設計できる。
- バイオレメディエーション - タンパク質-リガンドドッキングは、酵素によって分解される可能性のある汚染物質を予測するためにも使用できる。