＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞

 ＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞＞

ABSTRACT 

アロステリック制御の物理的起源は長らく未解明であり、特にポケット形態の変化がどのように遠隔部位の機能変化へと結びつくのかは明確な説明が存在しなかった。本研究では、タンパク質ポケット内に存在する“拘束水”に着目し、酵素アロステリーにおいてポケットの形態変化とエネルギー伝搬の関係にたいして酵素中の水の挙動が連動し随伴してふるまうのではないかという考えのもとにポケット水の運動方程式を解析した。拘束水を過減衰系として扱い、減衰方程式に基づく解析を行った結果、酵素蛋白の構造変化に伴って拘束水の運動エネルギ（量子的揺らぎを含む）が電場エネルギへと変換されることを数式的に示すに至った。さらに、ポケット体積・質量分布・距離に基づく自由度指標 λ が導かれ、構造変化に伴うλ の変化比（λ′/λ）がアロステリック効果の強度を定量的に記述しうることを発見した。これに基ずき水分子の振動運動エネルギーが変換され電場形成が発揮され、このことがアロステリの本体であることを既存データと照合し実証した。この物理モデルにもとづき、複数の生化学現象を統一的に説明することができた。

INTRODUCTION

タンパク質のアロステリック制御は、生体分子に普遍的に見られる現象であるにもかかわらず、その駆動源となる物理的エネルギーが何であるのかは未解明のままであった。結晶構造解析や分子動力学（MD）計算により、局所的な揺らぎがタンパク質全体へ伝播しうることは示されてきたが、これらの研究は主として“運動の相関”を記述するにとどまり、ポケット形態の変化がどの程度アロステリック強度に寄与するのかという根本問題には答えていない。

特に見落とされてきたのは、タンパク質ポケット内部に存在する“拘束水”の物理的役割である。ポケット内の水はバルク水とは異なり、配向自由度が低下し、双極子が整列しやすく、強い局所電場を形成しうる。しかし従来の生化学・MD などの研究では、水は単なる溶媒として扱われ、その電場形成能はアロステリック機構の要素として考慮されてこなかった。そのため、既報で観察されてきた Arg や Lys などの荷電側鎖の位置変化は、因果的説明を欠いた“偶発的揺らぎ”として解釈されてきた。

タンパク質ポケットの体積が減少すると、内部の水分子は表面効果により その配向次元が　3D から 2D、さらに 1D へと次元低下し、配向秩序を獲得して電気双極子として振る舞う。この現象はアクアポリン内部の 1D 水ワイヤーで観察される双極子整列と同一の物理機構であり、挟小空間における“運動エネルギー → 電気エネルギー”変換の一般形である。

本研究では、酵素の会合や薬剤結合に伴うポケット形態の変化が拘束水の双極子配向と電場形成を誘導し、その電場がアロステリック制御の駆動源となるという新しい物理モデルを提案し実証する。拘束水を過減衰系として扱い、古典的減衰方程式に基づく解析を行うことで、拘束水の運動エネルギ（量子的揺らぎを含む）が電場エネルギへと変換されることを数式的に示した。

さらに、ポケット体積・質量分布・特徴距離に基づく自由度指標 λ を導入し、構造変化に伴う λ の変化比（λ′/λ）がアロステリック効果の強度を定量的に記述しうることを示す。これにより、従来の MD 研究が扱えなかった「どの程度ポケットが変化すると、どの程度アロステリック効果が生じるのか」という“程度の問題”を、物理的に一貫した枠組みで解決する。

本研究は、拘束水の電場形成を介したエネルギ変換という視点から、ポケット形態とアロステリック制御を統合的に説明する物理モデルを提示し、既存データとの整合を実証した。これにより、アロステリー研究における長年のギャップを埋める定量的基盤を提供する。

酵素はバルク水中でブラウン運動を行い、並進運動だけでなく固有の振動モードを常に励起している。 本研究では、この酵素振動がポケット内部の水分子に与える 微小トルク に着目し、 拘束水の減衰挙動を物理モデルとして定式化し、 その作用形態がアロステリック制御とどのように結びつくかを検証する。

ポケット内の水分子の時間依存位置を X(t) とすると、 水分子は酵素振動により駆動されつつ、ポケット形態に由来する抵抗（インピーダンス）を受ける分子量や空間的サイズが大きく異なる多様な酵素スーパーファミリー間で普遍的に適用可能な枠組みを構築するため、我々はタンパク質マトリックスの固有振動タイムスケールに対して時間軸を再スケーリングした基準化モデル（Normalized Model）を導入する。この標準化されたレジーム下では、水アンサンブルの有効質量、および構造復元の基準力定数はともに 正規化される。。 この状況を振動減衰方程式として

と表す。

• 減衰項 X′ の係数はポケット体積と相関するため V を用い、 その他の減衰要因を λ にまとめた。

• 位置項 X の係数は、酵素の分子量 M と、酵素重心からポケット中心までの距離 d に基づくモーメント dM を用いた。

• さらに振動の大きさと相関する距離 d を掛け合わせ、ポケット形態依存性を明示した。

この方程式は、酵素が振動を続けていても、 ポケット内の水分子は相対的には振動を停止しうる ことを意味する。

Xを X=exp(−αt) と置換すると、

となり、特性方程式の判別式は

である。

この判別式は、拘束水の物理状態が D<0 / D=0 / D>0 の 3 つの領域に分かれることを示す。

• 解は複素数

• 振動しながら減衰

• エネルギーは熱として散逸

• 水分子は 3D 的に自由度を保持

これは バルク水の通常の減衰 に対応する。

• 振動モードが消失する直前

• バルク水的減衰と拘束水的秩序化が混在

• 酵素によってはアロステリーが成立し、別の酵素では成立しない

つまり アロステリーの有無が分岐する境界状態。

3．D＞０：過減衰（拘束水の秩序化＞電気双極子形成）

ここが本研究の核心である。

D>0 では：

• 解は実数のみ

• 振動モードが完全に消滅

• 自由度が「揺れ」から 「偏極」へと切り替わる

• 水分子は 3D → 2D → 1D へと次元低下

• 準結晶的な整列（ordering）が生じる

• 酵素振動エネルギーが 熱ではなく秩序化へ変換 される

• その秩序化が 双極子整列 → 局所電場形成 を生む

つまり D>0 は、　“運動エネルギー → 電気エネルギー” への変換が起こる領域である。CNT（カーボンナノチューブ）内の水の固体化現象と類似するが、 タンパク質ポケットではアミノ酸側鎖による表面科学的影響が強く、 正四面体的な氷構造ではなく 2D/1D の準結晶的整列 が形成される。この整列が 双極子ワイヤ（water wire） を生み、 アロステリック電場の源となる。

実効的にタンパク質へのアロステリック効果を検討する場合は D = 0 のときの λ を計算する。 つまり、

を用いる。

ここで重要なのは、λ は減衰の強さを表す量ではないという点である。

λ は 振動的にも非振動的にも振る舞える境界点（critical structural freedom） を表す。 実際に過減衰となるかどうかは λ ではなく、

• 温度

• 圧力

• 溶媒粘性

• 双極子応答

といった 環境パラメータ によって決定される。

したがって λ は 環境依存の減衰係数ではなく、構造そのものが持つ自由度の最大点 を示す指標である。

PDB 構造は結晶化条件により揺らぎが極限まで抑制されており、事実上この臨界点近傍に押し込まれている。このため λ は、生体環境の影響を受けない 純粋な構造的指標として抽出可能 となる。

現状で得られる λ の絶対値は、環境因子を含まない「構造自由度の臨界値」であるため、 実用上は相対比較が最も重要な意味を持つ。

• 酵素ごとの λ の分布

• 2 状態間の λ の大小

• Δλ（構造変化に伴う自由度の変化）

これらは、どの状態がより電場形成をなしうるか を示す強力な指標となる。

λ 計算に用いたポケット体積 V は、主として既報の MD 計算済みデータおよび CASTp による既報値を使用した。自ら MD 計算を行った例は少数である。既報で活性化状態が 1 種類しか得られない場合には、ポケット形成アミノ酸の特定の Cα–Cα 距離を長径・短径として複数の幾何学状態を設定し、それぞれにリガンドを付加した構造を AMBER により最適化した。得られた複数状態の比率から、既報に存在しない V を推定した。

既報値が全く得られない場合には、Vina の GRID を用いてポケット領域が立方体 GRID 内に収まるように設定し、GRID 体積からアミノ酸占拠体積を差し引くことでポケット水が占有可能な体積を算出した。アミノ酸占拠体積は、GRID 内のアミノ酸数を取得し、10 残基あたり 1300 ų として近似した。ポケット中心は GRID 中心とした。

注目すべき点は、この GRID 内アミノ酸差分体積近似 を用いても λ 計算値が既報の酵素活性状態を良好に説明したことである。これは λ が M,d,V の組み合わせで決定されるにもかかわらず、活性状態の変化が主として V の比率に依存する ことを示している。したがって、アロステリック効果はポケット水の状態変化によって規定されるという本研究の中心命題が、構造データにより強く支持された。

λ は M, d, V の組み合わせで決定されるが、 実際の活性状態の変化は主として V の比率に依存していた。 これは、アロステリック効果が構造変化そのものではなく、 ポケット水の状態変化（占有体積の変化）によって規定される という本研究の中心命題を強く支持する。

また、電場形成ポケットには 恒常的に電場を形成するポケットと 状況依存的に電場を形成するポケットが存在する。 λ 計算の対象となるのは後者であり、 刺激の有無によるアロステリック効果の差異は、 状況依存型ポケットの電場形成の相転移 によって説明できる。

判別式＝0 の λ 計算で活性化の有無を比較した際、 活性化が起こる場合には、 既知のポケット以外の電場形成ポケットの寄与 または ポケット形態ではなく性状（疎水性・水和自由エネルギー）の変化 により、臨界条件が過減衰から臨界減衰へと移行したことを示唆する。

拘束水が一次元双極子カラムを形成しうる「ポケット」を操作的に定義するため、ポケット表面の 曲率 と 連続性 の 2 条件を課した。

1. 曲率条件 ポケット表面を主軸に垂直な平面へ投影し、得られた輪郭線の局所曲率半径 Ri を三点円法により推定した。 Ri≤10 Å の領域を「水を拘束しうる十分な凹面」と定義した。

2. 連続性条件 cavity をポケットと分類するためには以下の 2 条件を満たす必要がある：

   • 全輪郭長の 50%以上 が Ri≤10 Å を満たす

   • Ri≤10 Å の連続領域が 12–15 Å 以上 存在する（単列水分子 4–5 個に相当）

これらの条件を満たさない cavity は、全体の凹面率が 50% を超えていても、双極子カラムを安定に保持できない断片的空間 として除外した。

従来の MD 計算では、酵素近傍の水分子の取り扱いが統一されておらず、計算負荷の制約から水の物理的寄与は十分に評価されてこなかった。本研究で導出した 過減衰条件 λ の解析により、ポケット構造変化が拘束水の電場形成を規定する主要因であることが明確になった。

さらに、この λ 条件は既報の多数の立体構造において満たされており、水の電場応答がアロステリック制御に普遍的に組み込まれている ことが示された。進化的に保存されたアミノ酸配列の“乱雑性”の背後に、電場形成の様式が保存されている可能性 も示唆される。

Dimensional analysis of λ

振動減衰方程式

の次元解析により、λ の物理的意味が明確化された。

減衰項 λVX′ の単位を整えると、式全体は

の次元を持つ。

基準体積のポケット水に対して、ポケット中心–重心距離の時間変化率を 1 としたとき、 単位時間あたりに λ g の水状態が変化する量 として読み替えることができる。

この量は 拘束水とバルク水の比率 に対応し、 すなわち ポケット内で新たに形成される電場の強さを直接表す指標 となる。

RESULTS

ここでは各酵素のD=０とした条件で既報の活性状態もしくは不活状態のλの値を比較することによって、アロステリックが機能している場合のほうがλ値が高いことを確認し実証する。λは構造変化情報以外の不測状態を含めたパラメーターでありアロステリック効果と関係することを示す。酵素ごとにカットオフ値はあるもののλは酵素活性の指標としての意義があることを確認する。

CASE①：GLP-１Rについて、PDB-ID:5VEXはLIGANDが結合していないholo状態で活性のない状態で、総分子量（M)106.32ｋD、pocket体積（V)：153459立方Å、pocket中心/分子重心距離（ｄ）：33.993Åである。λ＝0.004568である。これにGsとGLP-1が結合した状態はID:５VAIでＭ：161.76ｋD,　Ｖ：3950.218立方Å、ｄ：23.641Åで、λは0.1522であった。比でいうと約33倍活性化していた。LIGNDの結合により活性化したレセプターの様子をλ値の上昇でよく表している

CASE②：Rhodopsinについて、ID：3CAPは不活性型で、M:82.68kD, V:2469.283立方Å、d:24.575Åでλは0.2559であった。ID:1U19は11-cis-retinolが結合しているが不活性型で、M:82.68kD, V:2820.250立方Å、d:25.942でλは0.2404である。この蛋白ではこの程度のλでは活計化できない。ID：3POXは活性化型でMeta-２が結合している。M:41.33kD, V:2870.000立方Å、dは重心とポケット中心がほぼ重なるためもともとの振動方程式のモデルの物理的意味から１とした。λは12.85で活性化をよく表している。

CASE③：D2ドーパミンRについて、7LJCはly3154207が結合した　D2　dopamine　レセプターでGs蛋白が結合しているが不活性型である。λは低値で0.1698である。6CM4は λ = 0.2435 で　リスペリドンが結合したD2ドーパミンRであり、ややλ値が増加しているが活性化はしていない。6VMSは λ = 0.3094で、　D2ドーパミンRで活性があるタイプであり、λと酵素活性の相関を保証している

CASE④：β２アドレナリンRについて：β₂アドレナリン受容体（β₂AR）について、逆作動薬結合状態（2RH1）と G 蛋白複合体（3SN6）を比較したところ、自由度指標 λ は 0.1695から 0.7168へと約 4.2倍に増大した（λ′/λ = 5.27）。これは、β₂AR の活性化が TM6 の大規模な外向き変位と細胞内側ポケットの開口を伴うことと一致しており、GPCR におけるアロステリック活性化が「ポケット形態の変化を介した自由度の劇的な増大」として定量的に記述できることを示している。

CASE⑤：PDZ3 ドメイン（1BE9/1BFE）では、リガンド結合により λ は 0.01420 から 0.005887 へと減少し、λ′（1BFE)/λ (1BE9)= 0.0.41 であった。これは、結合ポケット周辺の自由度がリガンド結合によってむしろ制限される「拘束型アロステリー」であることを示唆し、PDZ における allosteric preorganization や entropic stabilization の報告と整合的である。

CASE⑥：AT1R における状態依存的 λ の変化、アンジオテンシン II 受容体 AT1R の活性型構造（6DO1）と不活性型構造（4YAY）について、 内部ポケットの幾何学量から算出した λ を比較した。 活性型 6DO1 の λ は 0.4788、不活性型 4YAY の λ は 0.1325 であり、 活性状態では約 3 倍の λ 増大が認められた。この差は、活性化に伴うポケットの伸長と体積減少により、 拘束水の配向性が高まり、局所電場が強化されることを示唆する。 AT1R は、構造状態の違いが λ を通じて拘束水電場に反映される典型例である。

CASE⑦：ヘモグロビンにおける T 状態と R 状態の λ の比較、ヘモグロビンの脱酸素型（T 状態, PDB 2DN2）と 酸素結合型（R 状態, PDB 1HHO）について、 内部ポケットの幾何学量から算出した λ を比較した。T 状態の λ は 0.36370、 R 状態の λ は 0.08342 であり、 R 状態より約 4.3 倍の λ 増大が認められた。T→R 変換に伴うポケットの伸長と体積増加により、 拘束水の配向性が弱まり、局所電場が低下されることを示す。これにより酸素を放出しやすくなる。ヘモグロビンは、構造状態の違いが λ を通じて 拘束水電場に反映される典型的なアロステリ例である。

CASE⑧：PKA C-subunit における活性型と不活性型の λPKA C-subunit の活性型構造（PDB 1ATP）と、阻害ペプチド PKI 結合型の不活性構造（PDB 1CMK）について、内部ポケットの幾何学量から算出した λ を比較した。 1ATP の λ は 0.3092、1CMK の λ は 0.02026 であり、活性状態では約 15倍の λ 増大が認められた。この増大は、活性化に伴うポケットの伸長と有効体積の減少により、拘束水の配向性が高まり、局所電場が強化されることを示す。 PKA C-subunit は、キナーゼにおいても構造状態の違いが λ を通じて拘束水電場に反映されることを示す代表例である。

CASE⑨：SRC kinase における λ の状態依存性、SRC kinase の自己阻害型構造（PDB 2SRC）と、活性化ループが開いた不活性型構造（PDB 1YI6）について、内部ポケットの幾何学量から算出した λ を比較した。 2SRC の λ は 0.09429、1YI6 の λ は 0.1156 であり、双方とも不活性であった。既報とも整合的であった。

CASE⑩：CCR5 における λ の状態依存性、CCR5 のアンタゴニスト結合型不活性構造（PDB 4MBS）と、Gi タンパク質結合型活性構造（PDB 6AKX）について、内部ポケットの幾何学量から算出した λ を比較した。 4MBS の λ は 0.03156、6AKX の λ は 0.03569であり、活性状態では約 1.13 倍の λ 増大が認められた。この増大は、活性化に伴うポケットの伸長と有効体積の減少により、拘束水の配向性が高まり、局所電場が強化されることを示す。 CCR5 は、GPCR においてもアロステリが λ を通じて拘束水電場に反映されることを示す一例である。

CASE⑪：CDK2 における λ の状態依存性、CDK2 の活性型構造（PDB 1FIN）と、不活性型構造（PDB 1HCK）について、内部ポケットの幾何学量から算出した λ を比較した。 1FIN の λ は 0.09743、1HCK の λ は 0.04321 であり、活性状態では約 2.3 倍の λ 増大が認められた。この増大は、活性化に伴うポケットの伸長と有効体積の減少により、拘束水の配向性が高まり、局所電場が強化されることを示す。 CDK2 は、細胞周期キナーゼにおいてもアロステリが λ を通じて拘束水電場に反映されることを示す一例である。

CASE⑫：KIT における λ の状態依存性、KIT 受容体の活性型構造（PDB 6DO1）と、不活性型構造（PDB 4YAY）について、内部ポケットの幾何学量から算出した λ を比較した。 6DO1 の λ は 0.4267、4YAY の λ は 0.1325 であり、活性状態では約 3.2 倍の λ 増大が認められた。この増大は、SCF 結合に伴うポケットの伸長と体積減少により、拘束水の配向性が高まり、局所電場が強化されることを示す。 KIT は、受容体型チロシンキナーゼにおいてもアロステリが λ を通じて拘束水電場に反映されることを示す代表例である。

CASE⑬：ZAP70 における λ の状態依存性、ZAP70 の不活性型構造（PDB 1U59）と、活性型構造（PDB 2OZO）について、内部ポケットの幾何学量から算出した λ を比較した。 1U59 の λ は 0.08765、2OZO の λ は 0.2144 であり、活性状態では約 2.5 倍の λ 増大が認められた。この増大は、活性化に伴う activation loop の開放とポケット体積の減少により、拘束水の配向性が高まり、局所電場が強化されることを示す。 ZAP70 は、免疫シグナル伝達におけるアロステリが λ を通じて拘束水電場に反映される典型例である。

CASE⑭：FYN における λ の状態依存性、Src-family kinase である FYN のリガンド非結合構造（PDB 2DQ7）と リガンド結合構造（同 2DQ7）について、内部ポケットの幾何学量から算出した λ を比較した。 リガンド非結合状態の λ は 0.3672、結合状態の λ は 0.3356 であり、 結合により λ は約 9% 減少した。この小さな減少は、FYN が HCK のような“ロック型スイッチ”ではなく、 揺らぎを保ったままシグナルを調整する“チューニング型キナーゼ”であることを示す。 FYN は、同じ Src-family 内でもアロステリの設計思想が λ により明確に区別される例である。

CASE⑮：PDGFRβ における λ の状態依存性、PDGFRβ のリガンド PDGF-B 結合型構造（PDB 3MJG）について、 Ig ドメイン 1–3 の溝構造から算出した λ を評価した。 3MJG の λ は 0.009285 であり、 同じ計算手法で得られた他の RTK（KIT、EGFR など）と比較して 低い λ を示した。これは、PDGFRβ の Ig 溝が比較的広く、 拘束水の配向性が弱く、局所電場が小さいことを反映する。 PDGFRβ は、RTK の中でも “低 λ 型” に分類され、 KIT（3.2倍増）、EGFR（高 λ 型）とは対照的な アロステリ設計の多様性を示す例である。

CASE⑯：Frizzled における λ の状態依存性、Frizzled の不活性型構造（PDB 4F0A）と、Wnt 結合型活性構造（PDB 4F0B）について、内部ポケットの幾何学量から算出した λ を比較した。 4F0A の λ は 0.01452、4F0B の λ は 0.03944 であり、Wnt 結合により λ は約 2.7 倍に増大した。この増大は、Wnt 結合に伴う CRD の大規模な再配置とポケット体積の減少により、拘束水の配向性が高まり、局所電場が強化されることを示す。 Frizzled は、RTK と GPCR の中間に位置する受容体であっても、アロステリが λ を通じて拘束水電場に反映されることを示す例である。

CASE⑰：CASE：Caspase‑9 における巨大複合体依存的 λ の劇的上昇、Caspase‑9 について、不活性状態（PDB 1NW9）と Apaf‑1 apoptosome に結合した活性状態（PDB 5WVE）の λ を比較した。 1NW9 における λ は 0.004034 と極めて低値であり、単体の Caspase‑9 がほぼアロステリック活性を示さないことを反映していた。一方、5WVE における λ は 16.516 であり、不活性状態の約 4,000 倍（λ′/λ ≈ 4096）という、これまでの解析例の中で最大級の増大を示した。この劇的な λ の上昇は、Caspase‑9 が単体では広いポケットと小さな有効質量を持ち、拘束水が液体的揺らぎ（D < 0）を維持しているのに対し、Apaf‑1 apoptosome に組み込まれることでポケットが強制的に圧縮され、拘束水が過減衰状態（D > 0）へと遷移することを示唆する。これにより双極子整列が進み、局所電場が急激に強化される結果、Caspase‑9 の活性化がスイッチ的に誘導される。Caspase‑9 は単体では活性化しないという生化学的特徴を持つが、本結果はその分子機構が 「巨大複合体依存的な λ の急峻な増大」として定量的に説明できることを示す。すなわち、Caspase‑9 は本研究で解析した中でも最も極端な “巨大複合体依存型アロステリー” を示す酵素であり、拘束水の過減衰化がアロステリック活性化の必須条件であることを強く支持する。

CASE⑱：SERCA Ca²⁺-ATPase における λ の状態依存性SERCA Ca²⁺-ATPase の E1（Ca²⁺結合・開状態）と E2（Ca²⁺非結合・閉状態）について、内部ポケットの幾何学量から算出した λ を比較した。E1（開・活性型, PDB 1SU4）：λ = 0.08015　E2（閉・不活性型, PDB 1IWO）：λ = 0.01334活性状態では 約 6 倍の λ 増大（λ′/λ ≈ 6.0） が認められた。 この増大は、E1 状態でポケットが伸長し有効体積が減少することで拘束水の配向性が高まり、局所電場が強化されることを示す。 SERCA は、酵素・受容体とは異なる 輸送体（transporter） においても、 アロステリック状態の違いが λ を通じて拘束水電場に反映されることを示す代表例である。

CASE⑲：Calmodulin における Ca²⁺依存的 λ の変化  Calmodulin（CaM）の Ca²⁺非結合型（Apo, PDB 1CFD）と Ca²⁺結合型（PDB 1CLL）について、内部ポケットの幾何学量から算出した λ を比較した。Apo 状態 1CFD の λ は 0.0007661、Ca²⁺結合状態 1CLL の λ は 0.001938 であり、Ca²⁺結合により λ は約 2.5 倍（λ′/λ ≈ 2.5）に増大した。この増大は、EF-hand の開口と疎水性ポケット形成に伴い、拘束水の配向性が高まり局所電場が強化されることを示す。Calmodulin は、酵素や受容体ではない Ca²⁺センサーにおいても、構造状態の違いが λ を通じて拘束水電場に反映されることを示す一例である。

CASE⑳：Hsp90 の代表的三状態に対して、構造的拘束度を示す正規化パラメータ λ を算出した。 対象としたのは、ATP 結合型の閉鎖構造（2XG9）、アポ型の開放構造（2IOP）、およびクライアント結合構造（5FWK）である。　得られた λ はそれぞれ 0.01137（2XG9）, 0.001011（2IOP）, 0.06959（5FWK） であり、三者の間に一桁以上の差が認められた。 ATP 結合型（2XG9）は、全体として均質な拘束を受けた閉鎖状態を示し、アポ型（2IOP）は最も低い λ を示すことで、広範な柔軟性を特徴とする開放状態であることが確認された。 一方、クライアント結合構造（5FWK）は最も高い λ を示し、ATP 結合とは異なる 局所的かつ非対称な拘束パターン を呈した。これらの結果は、Hsp90 が単純な「開 → 閉」の二状態モデルでは記述できず、クライアント結合により第三の力学的レジームが形成されることを示唆する。 この三相性は、Hsp90 が多様なクライアントを処理できる構造的基盤を与えると考えられる。

以上のように、LIGANDや会合酵素によりλの値が変化し、同じ酵素ではλがある一定上昇すると酵素活性が出現することがわかる。このことはすなわち、ポケット内水の挙動がアロステリック効果の遠隔的な影響に大きく関与していることを表している。

既報の MD 研究が相関的記述にとどまり、アロステリック強度の“程度の問題”を説明できなかったのに対し、本研究は λ′/λ という単一の指標がアロステリック効果の大きさを定量的に予測できることを実証した。これは、アロステリーの物理的起源を拘束水の過減衰ダイナミクスと電場エネルギ変換に求める本モデルが、既存のMD計算による説明を補完するもので、将来的にはポケット水の物理的な挙動も含めてシュミレーション計算すべきであることを示唆している。

以上より、λ はアロステリック効果の普遍的かつ定量的な指標として機能し、ポケット形態から機能変化を直接予測する新しい生物物理学的指標となっていることを実証した。本研究は、アロステリー研究の長年のギャップを埋めるとともに、創薬におけるアロステリック部位の予測とリード化合物の最適化に新たな道を開くものである。

この値は、 ポケット内の拘束水が臨界的に電場を形成する条件を満たすことを意味する。

さらに、 それぞれのポケットで生じる 電場の方向性が互いに整合するため、 2つの p52 分子は 静電的に引き寄せられ、ホモダイマーを形成すると考えられる。

一般に、タンパク質–タンパク質相互作用は MD 計算で得られる自由エネルギー低下によって説明されることが多い。 しかし本系では、 自由エネルギーという巨視的量ではなく、ポケット拘束水が生み出す“局所電場”が主要因である という点が特徴的である。

2．Glucokinase（GCK）— 活性型と不活性型における λ の比較　グルコキナーゼ（GCK）は、モノマー酵素でありながら顕著なアロステリック性を示す特異な酵素である。GCK のアロステリーは、従来より「スーパーループの構造変化」や「開閉運動」によって説明されてきたが、その物理的起源は未解明のままであった。本研究では、GCK の活性型構造（PDB: 1V4S）と不活性型構造（PDB: 1V4T）について、ポケット体積 V、分子量 M、重心–ポケット中心距離 d を算出し、拘束水の自由度指標 λ を比較した。

DISCCUSION

λ が解決した理論的ギャップ

本研究で導入した指標 λ は、既存のアロステリー研究が抱えてきた複数の理論的限界を統合的に解決する枠組みを提供する。従来の MD ベースの解析手法—DCCM、相互情報量、ENM、NMA、residue network 解析—は、タンパク質内部の運動の「相関」や「経路」を抽出する点で有用であった。しかし、これらの手法はいずれも、アロステリック効果を駆動する物理的エネルギ源を特定できない　“どの程度”アロステリーが生じるのかを定量化できないという根本的限界を共有していた。

本研究の初期段階では、アロステリーの本体はアミノ酸の微細な位置変化と側鎖配向の遠隔転移にあると考え、ポケット水を代表とする酵素内水の挙動がタンパク質の挙動を反映するとの考えの上にポケット水の位置エネルギーと運動エネルギーを記述するために減衰方程式を導入した。しかし解析の過程で、狭小空間における拘束水の 相転移の必然性 が明らかとなり、方程式の再解釈がアロステリーの本体の発見につながった。

すなわち、ポケット形態が拘束水の自由度（物理的状態）を規定し、双極子配向を通じて局所電場を形成し、その電場が側鎖配置と蛋白構造を変化させる という一連の過程である。

従来のアロステリーモデル—MWC、KNF、conformational selection、population shift—は、 アロステリーの概念的枠組みを築き、状態間の平衡や協同性を理解する上で重要な役割を果たしてきた。 しかし、これらのモデルは共通して、その背後にある具体的な物理量（自由度・拘束・エネルギー差）を明示的に扱うことができず、理論的にはそこで限界に達していた。

この限界は、アロステリーの駆動源を“ポケット水の物理量”として表現する枠組みが存在しなかったことに起因する。 すなわち、従来モデルは状態間の“関係”を記述することはできても、 その関係を生み出す物理的基底（V, d, M によって規定される拘束水の状態） を扱うことができなかった。この点で λ は、従来モデルを置き換えるのではなく、 それらを物理量に接続する補完的な軸 を提供する。 ポケット体積 V、特徴距離 d、質量分布 M を通じて、 状態間の自由度（物理的状態量）の差を直接数値化できるようになった。

本研究においてポケットは単なる幾何学的なくぼみではなく、拘束水が双極子を形成し電場を生成するための位相的単位である。この理由から、ポケットの定義には以下の二段階基準を採用した    ① 曲率半径 ≤10 Å の凹部が輪郭の 50% 以上            ②  12–15 Å 以上の連続凹部（single‑file water 4–5 個分）　この基準は、浅い凹みや断片的 cavity を排除し、電場形成が可能な“機能的ポケット”のみを抽出するためのものである。これは CNT やナノポア研究で示された「サブナノメートル拘束下でのみ水が強く秩序化する」という知見と整合する。　λ＝2 d √M / V は幾何学的に定義されるが、その機能的転移を決めるのは V そのものではなく、ポケットが双極子整列した水柱を保持できるかどうかである。また、内径 < 8 Å の cavity：水が並ばず電場が形成されない → λ は上昇しない　ということや　長径 > 14 Å の cavity：4–5 個の水が整列し電場が形成される → λ が非線形に跳ね上がる　という非線形性は、幾何学から直接生じる物理的シグモイドであり、アロステリック制御の on/off と根本的に関係する。このことは反対にポケットの変化の既存データをもとに、パターン認識を行うことにより新規のアロステリック誘導部位やリード分子を高確率に発見できることを意味している。創薬開発の現場に相当のインパクトを与えることにもなる。リード分子探査においても、先に酵素表面の分子に負荷をかけてポケット変化が起きることを確認した後に分子設計してゆくという手法にシフトがあるであろう。

本研究の枠組みは酵素スケールでは普遍的に成立したが、免疫受容体（TCR–pMHC など）への適用では例外的解釈を要し、統一的説明には至らなかった。これは、分子サイズ、界面形状、拘束水の次元性（3D→2D→1D）が階層をまたぐ際に異なるためである。細胞膜を貫通するような酵素においてリガンド結合部位、ポケット部位が膜を挟んでいる場合、もともとの微分方程式から変える必要もある。、今後、これらを体系的に解析することで、酵素から受容体に至るまでの 階層横断的エネルギー論 が構築されることが期待される。