Unveiling the Activation Mechanism of GLP-1 Receptor by an Ago-Allosteric Modulator via MD Simulations
J. Chem. Inf. Model. | 2026 | DOI: 10.1021/acs.jcim.6c00224
GLP-1R ago-PAM compound 2 がペプチド系/非ペプチド系アゴニストに対し異なるアロステリック経路で活性化を増強する機構を、μs MD + OPES + 動的ネットワーク解析で原子レベル解明
① 背景と課題:GLP-1R 活性化の動的描像

GLP-1R はクラス B GPCR の代表格で、2 型糖尿病・肥満治療の主要標的。Ago-PAM「compound 2」は単独で受容体を活性化しつつ、内因性ペプチド GLP-1 と非ペプチド低分子 LY3502970 の両者に対するアフィニティ・エフィカシーを増強する二重機能を示す稀有なリガンド。

クライオEM(PDB 7DUR/7DUQ/7E14)で複合体の静止画は得られているが、以下が未解明だった:

直交性アゴニストの種類によりアロステリック経路が変わるのか(経路依存性)、静的構造からは推定不可
不活性 (6LN2) → 完全活性 (7DUR) 遷移上の中間体の存在・自由エネルギーは未同定

→ μs スケール MD・OPES 拡張サンプリング・動的ネットワーク解析を統合して、二重機能の動的機構を初めて定量化することを目的とする。

② 手法の概要:4 系統合解析
  • 4 系の膜埋込 MD: 7DUR / 7DUQ / 7E14 / 6LN2
  • POPC 160 + コレステロール 40, 0.15 M NaCl, CHARMM36m, 310 K NPT
  • 主系で 5 μs プロダクション MD
  • OPES で不活性↔活性遷移の自由エネルギー面構築
  • 動的相関行列 → Floyd-Warshall 最短経路
  • Girvan-Newman でコミュニティ分解
解析パイプライン 4 cryoEM 系 7DUR/7DUQ/7E14/6LN2 CHARMM-GUI POPC/Chol 膜 5 μs MD CHARMM36m OPES サンプリング 不活性↔活性 FES 動的相関行列 Floyd-Warshall 経路 + コミュニティ Girvan-Newman → ago-PAM 二重機能の機構を抽出
③ 本研究で示したこと
  • Compound 2 単独で活性コンフォメーションを安定化(TM6 細胞内端の外向き変位を維持)
  • OPES で不活性→活性遷移上に 2 つの中間状態 を同定
  • Compound 2 結合部位 (C347 周辺) → 直交性ポケットへのアロステリック経路が GLP-1 系と LY3502970 系で異なる
  • 動的コミュニティ構造の差が、ペプチド/非ペプチド両アゴニストへの選択的増強の動的根拠
  • 静的 cryoEM 像では捉えられない経路依存性を定量化
④ (a) シミュレーション系構成
PDB結合リガンド
① compound 2 単独活性7DURcompound 2
② Compound 2 + ペプチド7DUQ+ GLP-1 (R36)
③ Compound 2 + 低分子7E14+ LY3502970
④ 不活性アポ6LN2
5 μs
プロダクション MD(系①)
POPC 160 / Chol 40
膜組成 + 0.15 M NaCl, CHARMM36m, 310 K NPT
④ (b) 不活性→活性 FES の中間状態
OPES 自由エネルギー断面 活性化座標 (TM6 外向き変位) ΔG 不活性 6LN2 中間 I1 中間 I2 活性 7DUR 2 つの中間状態を同定

Compound 2 結合系では活性側ベイスンへの集団偏りが顕著で、TM6 外向き変位・ECD 活性様配向・Gタンパク質界面の開放が安定維持された。

④ (c) アゴニスト別アロステリック経路
C347 → 直交性ポケット 最短経路 C347 compound 2 部位 R36 (GLP-1) LY3502970 部位 経路 P-pep(TM6 上方→ECD) 経路 P-sm(TM5/TM6 横断) ペプチド系/非ペプチド系で異経路

同一の compound 2 結合サイトから出発しても、相手アゴニストによって受容体内の動的相関ネットワークは異なる残基群を伝搬する。

④ (d) 活性コンフォメーション特徴の維持
活性指標スコア(compound 2 系) TM6 外向き 維持 G界面開放 維持 ECD 活性配向 維持 ECL 配向変化 確認 活性 (7DUR) 不活性 (6LN2)

5 μs MD で 4 つの活性コンフォメーション特徴がいずれも活性側に維持され、compound 2 単独でも活性化を引き起こす ago 機能を裏付けた。

⑤ テイクホームメッセージ
経路依存的アロステリー
ペプチド GLP-1 系と非ペプチド LY3502970 系で C347 から直交性ポケットに至る最短経路が異なる。同じ ago-PAM でも相手リガンドにより伝達ネットワークが切り替わる。
2 中間状態の自由エネルギー描像
OPES が不活性 (6LN2) → 活性 (7DUR) 遷移上に 2 つの中間状態を可視化。compound 2 結合は活性ベイスン側へ集団を偏らせる。
Compound 2 単独活性化の動的根拠
5 μs MD で TM6 外向き変位・G 界面開放・ECD 活性配向・ECL 配向変化が安定維持され、ago 機能を構造ダイナミクスで説明。
クラス B GPCR 創薬の指針
動的ネットワーク解析を ago-PAM スクリーニングや経路選択的モジュレーター設計の基盤として活用できる枠組みを提示。
応用補足:lib/md への展開
  • OPES ランナーを lib/md に追加(PLUMED 連携、CV 設計テンプレ)
  • 動的ネットワーク解析モジュール: 相関行列 → Floyd-Warshall → コミュニティ抽出のパイプライン化
  • クラス B GPCR テンプレートのベンチマーク系として 7DUR/7DUQ/7E14/6LN2 を採用
  • HBondAnalyzer / RMSDAnalyzer と連結して活性指標スコアを定量出力
インパクト
  • ago-PAM の二重機能を動的観点で初めて定量化
  • 2 型糖尿病/肥満治療のリード最適化に活かせる経路依存的設計指針
  • 動的ネットワーク解析を他クラス B GPCR へ系統的応用する道筋