Quantification of Hydrogen Bond Donating Ability of Biologically Relevant Compounds
J. Org. Chem. | April 2024 | DOI: 10.1021/acs.joc.3c02939
ピラジノンUV-Vis滴定で79化合物のHBD能を実験定量化。同じNH/OHでも置換基で値が大きく動くことを示し、計算pKaで区別できない部分構造差を可視化したCADD向け基礎データベース。
① 背景と課題

水素結合は薬物-受容体相互作用・膜透過性・溶解性・経口バイオアベイラビリティを左右する中心的な分子間力でありながら、HBD 能(hydrogen bond donating ability)を実験的に定量化する手段はこれまで限られていた。NMR 滴定・IR 振動数解析・ITC は精度を出せるが装置と操作が重く、計算側の代理指標である pKa・PSA・TPSA・HBD カウントは官能基の差を粗くしか拾えない。

同じ NH 基でも インドール / ピロール / アミン で受容体への寄与は大きく違うが、TPSA や HBD カウントでは 同一値で扱われてしまう
既存の実験法(NMR シフト・ITC)は 1 化合物あたり数時間〜の手間がかかり、ヘテロ環ライブラリ全体をスキャンするのは非現実的。

→ UV-Vis 滴定だけで HBD 能を Ka として読み取れる「ピラジノンセンサー法」を、薬物関連 79 化合物に適用してデータベースを拡張する。

② 手法の概要
  • Penn グループ開発のピラジノン比色センサーを CDCl₃ / CD₂Cl₂ 中で被測定化合物(analyte)と混合
  • analyte 濃度を段階的に増やし、アミド伸縮・n→π* 遷移由来の UV-Vis 吸光度シフトを取得
  • Benesi-Hildebrand プロットから複合体形成定数 Ka を決定(log Ka が HBD 強度の指標)
  • ヘテロ環(フラン・チオフェン・オキサゾール・ピリジン・インドール・ピロール・イミダゾール・トリアゾール)/ アルコール・フェノール・アミン・スルホンアミドを横断的にカバー
UV-Vis 滴定 → Benesi-Hildebrand → log Kₐ analyte + sensor UV-Vis 滴定 B-H plot 1/ΔA vs 1/[A] log Kₐ = HBD 強度 溶媒: CDCl₃ / CD₂Cl₂ センサー: ピラジノン (n→π*, アミド伸縮) 対象: ヘテロ環・薬物部分構造 79 化合物 UV-Vis 装置のみで HBD 強度ベンチマークを高スループット取得
③ 本研究で示したこと
  • 合計 79 化合物の HBD 能(log Kₐ)を一貫した条件で実測
  • 既知標準化合物の Kₐ を再測定し、先行報告値と定量的に再現(手法信頼性の担保)
  • 同一カテゴリ内 NH / OH でも、置換基次第で大きな HBD 強度差が生じることを実証
  • 計算 pKa(B3LYP/6-31G* レベル)とは定性相関に留まり、異性体間差を取り逃がしていることを実験で可視化
  • 電子求引性置換基で HBD 強化、電子供与性置換基で HBD 弱化という置換基効果のスケールを提供
④ (a) 官能基クラス別 HBD 強度マップ
官能基カテゴリ別の log Kₐ レンジ(相対スケール) 非常に強 脂肪族 OH フェノール OH フラン / チオフェン CH ピロール NH インドール NH イミダゾール / トリアゾール NH スルホンアミド NH 脂肪族 NH (アミン)

スルホンアミド NH や アゾール NH は計算 HBD カウントでは「1 個」と単純に数えられるが、実測では脂肪族 OH と比べ明確に強い HBDとして振る舞う。

④ (b) 置換基効果(インドール NH)
置換基 σ → log Kₐ 変化(模式) 置換基 Hammett σ log Kₐ (相対) 無置換 5-OMe 5-Me indole (基準) 5-Br 5-CN 5-NO₂ EWG → HBD 強 EDG → HBD 弱

電子求引性 (EWG) と HBD 強度はおおむね単調増加。インドール骨格内の 5 位置換だけで HBD 能が一段階以上動き、計算 HBD カウントでは見えない設計余地が広い。

④ (c) 計算 pKa vs 実測 log Kₐ
計算 pKa は定性相関のみ・異性体差は捉えきれない 計算 pKa (B3LYP/6-31G*) 実測 log Kₐ 理想 1:1 同一 pKa でも log Kₐ が散る 脂肪族 OH フェノール ピロール/インドール NH スルホンアミド

計算 pKa は弱→強の順序はおよそ合うが、同一 pKa の異性体間で実測 log Kₐ が大きく散る領域があり、ここを実測値で埋めることに本論文の価値がある。

④ (d) データセット構成
カテゴリn(概算)
5員環ヘテロ環 NHピロール / インドール / イミダゾール / トリアゾール~22
5員環 ヘテロ環 (CH/O/S)フラン / チオフェン / オキサゾール~12
6員環 N ヘテロ環ピリジン誘導体・関連~10
スルホンアミド NHaryl-SO₂NH-R~8
フェノール OH置換フェノール (EDG/EWG)~10
脂肪族 OH / NHアルコール・アミン~17
合計79 化合物
79 CDCl₃ / CD₂Cl₂ 中で取得した log Kₐ エントリ

同一手法・同一条件で取得されているため、官能基間の横串比較が可能。スルホンアミド NH は計算 HBD カウント=1 と扱われがちだが、実測では強 HBD 群に属する。

⑤ テイクホームメッセージ
HBD は「数」ではなく「強さ」で見る
同じ NH / OH でも置換基により log Kₐ が大きく動く。HBD カウントや TPSA に頼った最適化は粒度が粗く、置換基チューニングの設計余地を取りこぼす。
UV-Vis 1 台でベンチマーク化可能
NMR・ITC を持ち出さずとも、ピラジノンセンサー + UV-Vis 滴定 + Benesi-Hildebrand で log Kₐ が出る。シリーズ展開時の HBD パラメータ取得が一気に現実的になる。
計算 pKa の死角を埋めるデータ
計算 pKa は強弱の順序は当てられるが、同一 pKa でも実測 HBD が散るゾーン(特にアゾール NH と スルホンアミド NH)が存在し、その差を実測で初めて切り分けられた。
スルホンアミド NH は過小評価されがち
計算側で弱く見積もられやすいスルホンアミド NH が、実測では強 HBD 帯に位置する。リガンド設計で HBA とのペアリング候補として再評価する価値が高い。
応用: 計算化学への取り込み
  • lib/docking: ProLIFCalculator の HBD 接触に log Kₐ を重み付けし、単純カウントでなく強度加重スコアへ拡張
  • lib/molgen: MolgenYaml の制約に「HBD 強度ターゲット」を追加し、同一 SMARTS の異なる置換基を log Kₐ 推定で差別化
  • 新スコアラー: 79 化合物データを fragment-level lookup として持ち、生成分子の HBD 強度プロファイルを即時推定するフィルター
インパクト
  • HBD カウント / TPSA を超える、置換基まで効く HBD 強度パラメータが手に入る
  • スルホンアミド NH の再評価により、CADD のスコア関数を局所的に再校正できる
  • UV-Vis のみで増設可能な実験プロトコル → 社内シリーズで HBD ベンチマーク取得を継続できる
限界: 溶媒は CDCl₃ / CD₂Cl₂ で水性条件は未評価/79 化合物は ML モデル学習にはまだ小さい/受容体ポケット幾何依存性は反映されない。