MEPプラグインで分子の静電ポテンシャルマップを読む

なぜ構造だけでなく電荷を見るのか？

骨格構造式では電子がどこにあるかはわかりません。しかし電荷分布は多くの化学を説明します。水が極性を持つ理由、水素結合ができる場所、求核剤が攻撃する炭素などです。

教科書では、カラフルな静電ポテンシャルマップでこれを示しています。これらはふつう量子化学ソフトで作られ、準備にも計算にも時間がかかります。MEPプラグインを使えば、描いたどんな分子についても、同じ図の近似版をすぐに得られます。

MEPマップが示すもの

ある点での分子静電ポテンシャル（MEP）とは、そこに置いた小さな正の試験電荷が、分子のすべての原子核と電子から受けるエネルギーです。マップはこの値を分子の周りの表面上に色で表します。

• 赤は負のポテンシャルを表します。電子豊富な領域で、正電荷、プロトン、求電子剤を引きつけます。酸素や窒素の非共有電子対が典型的な赤い部分です。
• 青は正のポテンシャルを表します。電子不足の領域で、負電荷や非共有電子対を引きつけます。O–HやN–Hの水素が典型的な青い部分です。
• その中間の淡い色はほぼ中性で、たとえばC–Hの大部分の領域がこれにあたります。

MolDrawでの近似方法

プラグインはGasteiger（Gasteiger–Marsili）部分電荷を計算します。この方法では、結合した原子間で電気陰性度に応じて電荷を少しずつ分配し、収束するまで繰り返します。原子と結合だけを使うので高速で、ブラウザ内で動作します。得られた電荷で3Dモデルの周りの溶媒接触可能表面を色分けします。

これは定性的な図です。どの原子がより負か、より正かの順位はわかりますが、引用できるポテンシャルエネルギーの値は得られません。

プラグインで表示される情報

MEPプラグインのインストール方法

1. MolDrawエディターを開いてサインインします。無料アカウントで十分です。
2. Settings（歯車）をクリックし、Pluginsをクリックします。
3. Electrostatic Potential (MEP)の横のInstallをクリックします。解析パネルにMEPタブが表示されます。

MEPマップの作り方と読み方

1. 分子を描くか、PubChem検索ボックスで名前から読み込みます。
2. MEPタブを開きます。電荷がすぐに計算されます。
3. 電荷の表の上のほうを読みます。これらが最も電子豊富な原子です。
4. Show in 3Dを押します。色分けされた表面付きで3D表示が開きます。
5. 表面の下の原子が見えるまでOpacityスライダーを動かします。モデルを回転させ、赤や青の部分がどの原子に対応するかを確認します。
6. 描画を変更したらRefresh chargesを押します。終わったらTurn surface offを押します。

授業や研究室での活用法

• 極性：水、メタノール、炭化水素を順に表示します。どの分子に強く色付いた領域があるか、学生はすぐにわかります。
• 水素結合：赤い酸素原子や窒素原子はアクセプターになりやすく、青いO–HやN–Hの水素はドナーになりやすい部位です。原子をアクセプターにしている非共有電子対を示すには、1つの原子または分子全体に非共有電子対を追加しましょう。
• カルボニルの反応性：アセトンやエステルを描きます。酸素は負、カルボニル炭素は正で、そこが求核剤の攻撃する場所です。
• 置換基効果：フェノールとアニリン、あるいはニトロ基の有無で分子を比べ、電荷がどう移動するかを見ます。
• スライドと配布資料：色分けされた表面は、講義で抽象的な考え方をわかりやすく見せてくれます。後ろの席からもヘテロ原子が目立つ平面的なスライドには、2D描画を色付きの球棒で表示しましょう。

ヒントと率直な限界

• 1つの分子内の原子どうし、またはよく似た分子どうしで比較してください。数値をエネルギーとして読まないでください。
• Gasteiger電荷は描いた結合に依存し、3D形状には依存しません。結合を回転させても変わりません。また、描いたプロトン化状態にも依存するので、まずpKaプラグインで目的のpHでの主な荷電形を確認しましょう。
• 原子中心の電荷では、芳香環の電子豊富な面や重いハロゲンの正のキャップ（σホール）など、実際のポテンシャルの特徴をすべて表すことはできません。
• 論文品質のマップには量子化学プログラムを使ってください。このプラグインは直感を養い、すばやく比較するために使いましょう。

その他のMolDrawプラグイン