PRISMのEC50をそのまま次の実験濃度に使ってよいか

Daiki Kumakura · 2026-10-28 · Technical note · 1.0

要旨

PRISM 19Q4から事前選択した3薬・肺由来6細胞株を、制約付きHill曲線、条件付き残差profile、technical redoで再解析した。予定30組のうち1組は全欠測で、残る29曲線の12曲線はEC50 profileが探索境界まで閉じなかった。最適化の収束は識別性を保証せず、閉じた範囲も独立した生物学的再現性を意味しない。この対象からPRISM全体の信頼性は推定できない。

背景と目的

公開スクリーニングのEC50は、薬剤と細胞株を絞る入口になる。しかし、数値が表にあることと、その値の近くを追試すればよいことは別である。濃度反応の中点が測定範囲にあるか、反復が残っているか、別screenでも曲線が保たれるかを確認する必要がある。

本記事ではPRISMの公式19Q4データから、事前に選んだ3薬・肺由来6細胞株を再解析した。30組のうち1組は全欠測、残る29曲線では12曲線のEC50 profileが探索範囲の端まで閉じなかった。この12曲線では、点推定を精密な感受性の順位や追加濃度の根拠として扱えない。 一方、erlotinib–HCC827のようにtechnical redoでも遷移が保たれる例もあった。目的はPRISM全体の良し悪しを判定することではなく、次に何を測るべきかを曲線ごとに分けることである。

方法

問いと使用範囲

QOIは「公開EC50を次の濃度・再実験の選択に使えるか、その判断を変えるのは何か」とする。COUは非臨床の追試設計に限る。患者のPKを推定した解析ではなく、細胞株の結果から臨床用量や治療選択は決めない。誤ったEC50の採用は、不適切な濃度に追試を集中させる危険がある。

データと前処理

PRISMは細胞株を混合し、barcodeを通じて薬剤処理後の相対的なviabilityを測るスクリーニングである (1)。 公式19Q4 version 4を利用した。Figshare上のライセンスはCC BY 4.0であり、出典を明示する (2)。

読込前に、肺由来・STR確認済みの細胞株をDepMap ID昇順で6株選び、erlotinib、palbociclib、doxorubicinを固定した。薬剤反応や遺伝子変異で選んでいない。ただし、この小さな便宜的対象はPRISM全体を代表しない。

表1　選択対象と観測値

入力はreplicate未集約のlog2 fold changeであるが、すでにQC・対照正規化・ComBat補正を受けている (2, 3)。 生の蛍光測定値ではない。検出plateによって残る測定が異なり、最大3plateが残る曲線と1plateしかない曲線を同列の独立反復数として数えない。

erlotinibとpalbociclibにはHTS002とMTS010の両方がある。READMEはMTS010をtechnical redoと記載し、利用可能なら優先するよう勧めている (2)。 したがって、両screenの比較を独立した生物学的追試や外部検証とは呼ばない。

曲線推定と条件付きprofile

相対viabilityを \(v=2^{\mathrm{log2FC}}\) とし、クリッピングせずに次を推定した。

\[ v(c)=L+\frac{U-L}{1+(c/EC_{50})^h}+b_{plate} \]

\(L,U\) は高濃度・低濃度側の漸近値、\(h\) はHill slopeである。\(b_{plate}\) はscreen内の平均が0になるplate offsetとした。少数plateから集団分散を推定する混合効果モデルではない。\(0\le L\le U\le2.5\)、\(L\le1.5\)、\(0.1\le h\le5\)とし、EC50の探索を測定範囲の両側3桁まで延ばした。複数の開始点で残差平方和を最小化した。

EC50は上下の漸近値の中点であり、viabilityが0.5になるIC50とは異なる。 \(L>0.5\)なら、EC50が有限でも絶対50%阻害濃度はこのモデル上存在しない。

次にEC50を固定し、残りのパラメータを再推定するprofileを計算した。全域101点に推定値近傍81点を追加し、各点を4開始点で最適化した。残差増加をF分布に基づく95%相当の基準と比べた。この基準は独立・等分散で近似的に正規な残差を想定し、制約付き非線形推定では近似となる。これは補正済み技術反復、独立残差、指定曲線を条件とした近似的な許容範囲であり、独立した実験の95%信頼区間ではない。境界制約や補正後の相関もあるため、coverageを保証しない。

結果

遷移を観測している曲線と数値だけが出る曲線

図1。 左のerlotinib–HCC827は両screenで濃度依存的な低下が見える。中央のNCIH1581では反応がばらつき、曲線の遷移を決めにくい。右のpalbociclib–HCC827では、HTS002の遷移に対してMTS010の点が大きく散らばる。凡例のnは有限な技術測定値の数であり、独立実験数ではない。

表2　代表的なEC50推定値と条件付きprofile範囲

29曲線で最適化は収束したが、12曲線のprofileが探索境界まで開放し、7曲線の点推定は測定範囲外だった。17曲線では漸近値またはslopeが設定した境界に接し、7曲線は1plateしか残っていない。収束したことだけでは、EC50がデータから識別されるとは言えない。 また、閉じたprofileでも1plateだけから推定されたものがあり、狭い許容範囲を独立再現性の証拠とはしない。

点推定の変化と推定できないことを分ける

図2。 HCC827のerlotinibでは許容されるEC50が限定される。一方、中央の平坦なprofileでは、非常に異なるEC50でも同程度に観測値を説明できる。グラフの縦軸は表示を0〜2に限定しているが、判定は保存した全profile値で行った。

同じ薬・株を比較できた11組のうち10組では、screen間のEC50点推定が4倍を超えて動いた。ただし、未識別のEC50も含むため「感受性が4倍変わった」とは解釈しない。閉じたprofileが両screenにあるのは3組に限られ、その範囲でも条件付き推定の違いと技術的再現性の問題は残る。

公式処理済みパラメータと結合できた27曲線も比較した。今回の制約、plate offset、最適化によって点推定は変わるため、公式EC50をそのまま再現した解析ではない。同じ観測源からの再処理比較であり、独立検証ではない。両表を比較できるよう、点推定・profile・反復構造を一緒に保存する。

考察

局所感度と追加濃度

図3。 上下漸近値とHill slopeが既知なら、EC50付近の濃度はその変化に敏感になる。これらも不確かな曲線では、局所感度の最大点だけで最適な追加濃度は決まらない。

erlotinib–HCC827の約0.01〜0.03 µMには、両screenで遷移が見える。一方、profileが開く曲線や1plateしか残らない曲線では、点推定の近くへ濃度を追加する根拠が弱い。図の感度は、対照や漸近値、別実験での再現性が不明という問題を解消しない。

限界

対象は3薬・6株であり、PRISM全体の未識別率を推定する標本ではない。補正前のMFIからQC・ComBatを再実装しておらず、plate offsetだけでは処理の不確実性を復元できない。残差の独立性・等分散性、Hill曲線の妥当性も十分ではなく、profileを独立実験CIとして読むことはできない。培地中の遊離濃度、曝露時間差、患者PKへの変換も扱っていない。

結論

公開EC50を追試設計へ使うには、その曲線でどの濃度範囲が観測され、どの反復が残り、EC50がどこまで識別されるかを確認する必要がある。安定した遷移の周辺へ濃度を増やす場合と、元の範囲から再測定する場合を分けられたことが、今回の実測データ解析の結論である。

データとコードの公開状況

数値はanalysis/02_analyze.pyのcurve_fits.csv、profiles.csv、audit.json、screen_comparison.csvから取得した。図はanalysis/03_figures.pyで作成し、再実行手順とNotebookを解析作業フォルダに保存した。

配布関数

解析関数と単体検査のZIPを配布する。原データや解析全体の再実行資料は含まない。検査に使う構成例を、臨床データの検証とは扱わない。

参照

1. Corsello SM, Nagari RT, Spangler RD, Rossen J, Kocak M, Bryan JG, et al. Discovering the anticancer potential of non-oncology drugs by systematic viability profiling. Nature Cancer. 2020;1(2):235–248. doi:10.1038/s43018-019-0018-6. PRISM project.
2. Broad DepMap. PRISM Repurposing 19Q4 Dataset [Dataset]. Version 4. figshare. README, annotations, uncollapsed log2 fold changes and processed curve parameters. CC BY 4.0. Versioned dataset.
3. Broad Institute. PRISM secondary screen processing script [Software]. Official processing code.

All articles