IC50の4パラメータロジスティック(4PL)フィッティングの求め方

そのときの標準的なモデルが4パラメータロジスティック(4PL)です。上限・下限・Hill係数・IC50の4つで、濃度と応答の関係を作るS字（シグモイド）を表します。式は概ね「応答 ＝ 下限 ＋ (上限 － 下限) ÷ (1 ＋ (濃度 ÷ IC50) の Hill係数 乗)」の形で書かれ、この4つを同時に動かして残差平方和が最小になる曲線を探します。

問題は、あてはめ方しだいで同じデータからでも違うIC50が出てくる点です。上下限を勝手に固定してよいのか、点は何段階どこに置くのか、片側の平坦域が見えていないときにどう扱うのか。本稿では4PLの意味と求め方を確認し、制約設定・希釈系列・外挿の注意、そしてIC50とEC50・Ki・Kdの違いまでを整理します。

4PLの4つのパラメータが表すもの

まず4つのパラメータが曲線のどこを担うのかを押さえると、あてはめの良し悪しを判断しやすくなります。

ここで大切なのは、IC50が「上限と下限の中点」で定義される相対的な量だという点です。上限や下限の推定がずれれば、中点であるIC50も連動して動きます。Hill係数は遷移の急さを表し、絶対値が大きいほど狭い濃度域で応答が切り替わります。1から大きく外れる場合は、複数の結合部位や非特異的な効果、あるいはあてはめの不安定さを疑う材料になります。

4PLはHill式（シグモイド用量反応）と同じ枠組みで、細胞生存アッセイや酵素阻害アッセイでIC50を出す際の定番として広く使われています。 IC50は生データではなく、上限と下限を含む曲線全体から相対的に決まる推定量です 。

上限・下限を0/100に固定してよいか

正規化したデータでは、無処理を100%、完全阻害を0%に置き、上限100・下限0に固定して2パラメータで解く方法がよく採られます。手軽ですが、固定してよいかは条件に依存します。

固定が妥当なのは、対照ウェル（ブランクと無処理）で平坦域が実測により裏づけられている場合です。逆に、データがどちらかの平坦域まで届いていないのに上下限を勝手に固定すると、見かけ上きれいに合っても、IC50が実態からずれる原因になります。

• 上限：生データの頭打ちが対照値と一致しているかを確認します。合っていれば固定の根拠になります。
• 下限：最大阻害の平坦域が見えているかを見ます。見えていなければ、下限は対照から拘束するか、推定として慎重に扱います。
• どちらも平坦域が観測されていない場合、上下限は実データに支えられず、モデルが一意に決まりにくくなります。

無処理（濃度ゼロ）の対照は正規化の基準として役立ちますが、対数横軸には載せられないため、あてはめ自体からは外すのが無難とされます。 平坦域が実測で支えられているときに限って上下限の固定は正当化される、というのが基本の考え方です 。片側でも平坦域が欠けているなら、その側は対照から拘束するか、範囲を広げて取り直すのが安全です。

点数・希釈系列・外挿の注意

あてはめの質は、どの濃度に何点置くかで大きく変わります。用量反応の設計では、少なくとも3桁（3 log）以上の濃度域にわたって概ね8〜12点、各濃度に2〜3反復を置く構成が目安とされます。

希釈系列には落とし穴があります。段階希釈は誤差が末端に向けて積み上がりやすく、ある試算では、ピペッティングに10%の誤差があると12点系列の反対端で40倍もの濃度ずれになり得るとされます。これは希釈設計しだいでIC50が系統的に動く余地があることを意味します。

外挿はさらに注意が要ります。応答が遷移域に入っていない、つまり平坦域が取れていないデータでは、測定ノイズしだいで多数の異なるS字がほぼ同等に当てはまり、推定されるIC50が大きく揺れます。この場合に出るIC50は、実測に支えられた値ではなく曲線外挿の産物になりがちです。こうしたばらつきの評価には、アッセイ品質の指標としてZ因子などで系そのものの再現性を先に確かめておくと、フィッティングのブレとアッセイのブレを切り分けやすくなります。 両側の平坦域を実測で捉えることが、外挿に頼らないIC50を得る前提です 。

IC50とEC50・Ki・Kdの違い

近い記号が並ぶため混同しやすいのですが、指す対象が違います。

IC50は阻害系、EC50は活性化系という向きの違いがまずあります。より本質的なのは、IC50が実験条件に左右される相対値である一方、KiやKdはより条件に依存しにくい定数だという点です。IC50は基質濃度や標的量などの条件で動くため、系をまたいで直接比較しにくい面があります。

競合阻害では、IC50とKiはCheng-Prusoff式で結びつけられます。基質（またはリガンド）濃度とその親和性が分かっていれば、IC50からKiを見積もれます。この関係から、競合阻害ではIC50は常にKiより大きくなります。ただしCheng-Prusoff式は条件によって過大・過小評価を生むとされ、低濃度域などでは別の解析が推奨されることもあります。

KdはSPRやBLIといった結合速度論の測定から結合・解離の速度定数を通じて求められ、酵素阻害系ではKiがKdに近い意味を持つ場面もあります。 IC50は使い勝手のよい相対指標ですが、条件をまたいだ親和性の比較にはKiやKdへ橋渡しして考える必要があります 。数値を報告するときは、どの系・どの条件で得たIC50かを併記しておくと、後の比較で誤解を避けられます。

まとめ

IC50は用量反応をひとつの数字にまとめる便利な指標ですが、その値は4PLで引いた曲線の上に乗った相対量です。上限・下限・Hill係数を含む曲線全体を、両側の平坦域まで実測で支えられているかが精度を左右します。上下限を0/100に固定してよいのは平坦域が対照で裏づけられているときに限られ、希釈系列の誤差や平坦域を欠いた外挿は、見た目のあてはまりに反してIC50を系統的に動かします。

そして条件をまたいで親和性を比べたいなら、IC50単独ではなくKiやKdへ橋渡しして解釈します。IC50は「どの系で、どんな制約で引いた曲線から得た値か」まで含めて意味を持つ、という前提を外さないことが、ふるい分けの結果を再現可能にする近道です。

参考文献

• OECD, OECD Guidelines for the Testing of Chemicals
• ICH, ICH Guidelines
• FDA, Analytical Procedures and Methods Validation for Drugs and Biologics
• EMA, Scientific Guidelines