A DFT study of the axial coordination of metomidate to iron-porphyrin

技術要約：鉄ポルフィリンへのメトミデートの軸方向配位に関するDFT研究

問題提起
麻酔薬エトミデートの構造類似体であるメトミデートは、乱用の可能性とそれに伴う公衆衛生上のリスクから懸念される物質として浮上しており、中国において規制物質に指定されている。現在の検出法はクロマトグラフィー質量分析に依存しているが、これは高コストであり、専門的な実験環境や訓練を受けた人員を必要とするため、迅速なオンサイト・スクリーニングには限界がある。費用対効果の高い代替検出戦略を開発するためには、メトミデートと潜在的なセンシングプラットフォームとの間の特定の分子相互作用に関する基礎的な理解が必要である。鉄ポルフィリンは、選択的な軸方向リガンド結合能と、イミダゾールのような窒素含有塩基に対する明確な電子応答性を有することで知られており、このような用途のための有望な枠組みを提示している。しかし、メトミデートの鉄ポルフィリンへの配位における具体的な結合メカニズム、エネルギー論、および電子的な帰結については、理論的に未解明のままである。

手法
本研究では、密度汎関数理論（DFT）計算を用いて、5種類の鉄ポルフィリンモデル（FeTPPCl、FeTPyPCl、FeTBnPCl、FeT(p-NO2Ph)PCl、およびFeT(p-OMePh)PCl）に対するメトミデートの軸方向配位を系統的に調査した。これらのモデルは、単純なアリール環、ヘテロアリール環、メチレン分離されたアリール環、および電子吸引性または電子供与性基を含む、幅広いメソ置換体を代表するように選択された。

すべての量子化学計算はORCA 6.1.1を用いて行われた。構造最適化にはPBE0-D4関数とdef2-SVP基底系を用い、一点エネルギー計算にはwB97M-V関数とdef2-TZVP基底系を用いた。基底関数重なり誤差を排除するために、GrimmeのGCP補正を適用した。励起状態計算は、ジクロロメタン条件下をシミュレートするために、CPCM溶媒モデル、B3LYP-D4関数、およびdef2-TZVP基底系を用いたTDA法を用いて実施された。

その後の波動関数解析では、Multiwfn 26.2.2を用い、以下の手法によってシステムを特性評価した：

• 静電ポテンシャル（ESP）およびFukui関数： メトミデートの領域選択性と反応性を予測するため。
• 原子内解析（AIM）： 結合臨界点（BCP）におけるトポロジカルな特性を解析し、結合の性質を確認するため。
• 相互作用領域インジケーター（IRI）： 微弱な非共有結合的相互作用を可視化および特定するため。
• 電荷移動スペクトル（CTS）およびホール・エレクトロン解析： 励起状態の特性および電子移動メカニズムを調査するため。

主な結果

1. 反応性と領域選択性： ESP解析および凝縮Fukui関数計算により、メトミデート中のイミダゾール窒素原子（N4）が最も親核的な部位であることが特定された。この部位は、イミダゾール自体の窒素と同等、あるいはそれ以上に反応性が高い。この部位が配位の主要な箇所となる。
2. 結合エネルギーとスピンクロスオーバー： 鉄ポルフィリンへのメトミデートの配位は自発的かつ発熱的であり、結合エネルギーは-15.34から-19.52 kcal/molの範囲であった。明確なスピンクロスオーバー現象が観察された。すなわち、5配位の前駆体錯体は高スピンのセクステット状態を好むのに対し、6配位のメトミデート錯体は低スピンのダブレット基底状態をとる。これは、メトミデートの配位が鉄中心におけるリガンド場を強化することを示している。
3. 結合の性質（AIM解析）： AIM解析により、鉄とメトミデート窒素の間の配位結合の形成が確認された。電子密度（$\rho$）が0.076から0.079 a.u.の間であり、負のエネルギー密度（$H$）が約-0.012 a.u.である(3,-1)結合臨界点（BCP）の存在は、安定化相互作用を裏付けている。$|V|/G$比（1.116から1.123の範囲）は、Fe–N結合が有意な部分共有結合性を有していることを示している。ラプラシアン結合次数（LBO）解析は、Fe–Nの共有結合強度と全体的な結合エネルギーとの間に強い相関があることを示した。
4. 微弱な相互作用（IRI解析）： 主要なFe–N結合を超えて、IRI解析はポルフィリンの電子特性によって調整される微弱な相互作用のネットワークを明らかにした：
   • イミダゾール基とポルフィリンマクロサイクル間のCH-$\pi$相互作用がすべての錯体に存在した。
   • $\pi$-$\pi$スタッキングのパターンは多様であった：電子不足のポルフィリン（例：FeT(p-NO2Ph)PCl）では、メトミデートのフェニル基とフェイス・トゥ・フェイス（面対面）のスタッキングが起こったが、電子豊富なポルフィリン（例：FeT(p-OMePh)PCl）ではT字型スタッキングが観察された。
   • これらの二次的な相互作用、特に電子不足系におけるフェイス・トゥ・フェイスの$\pi$-$\pi$スタッキングは、結合エネルギーの変動に寄与している。
5. 励起状態の特性： 励起状態計算により、メトミデートの配位がポルフィリン由来の遷移（ソレー帯およびQ帯）の振動子強度を増強することが示された。電荷移動スペクトル解析は、リガンドとポルフィリンマクロサイクル間の双方向の電子移動を明らかにし、ポルフィリンからメトミデートへの電子供与が支配的であった。この効果は、電子豊富なFeT(p-OMePh)PCl錯体において最も顕著であった。

意義および主張
著者らは、本研究が原子および電子レベルにおいて、メトミデート–鉄ポルフィリン・ホストゲスト系の結合メカニズムと電子相互作用を解明したと主張している。本研究は、相互作用が主要な部分共有結合的Fe–N配位結合によって安定化され、複数の補助的な微弱な相互作用（CH-$\pi$および$\pi$-$\pi$スタッキング）および明確なスピン状態遷移を伴うものであることを確立している。

本論文は、これらの知見が、鉄ポルフィリンによるメトミデートの分子レベルでの認識基盤を理解するための体系的な理論的基礎を提供するものであると断じている。したがって、これらの結果は、当該システムの特異的な配位化学およびその結果生じる明確な光学的・電子的応答（振動子強度の増強や電荷移動特性など）を活用した、潜在的なメトミデート・センサー設計のためのガイドとして提示されている。本研究は、具体的な実験的センサーデバイスを提案するものではなく、むしろそれらの将来の開発に必要な理論的根拠を提供するものである。