モスネツズマブの公開PopPKモデルを再実装する：段階的増量の曝露比較とRO–CRS外挿の限界
On this page
何を調べたか
モスネツズマブ（mosunetuzumab）の段階的増量について、公開された母集団薬物動態（PopPK）モデルをPythonで再実装し、初回投与量・2回目投与量・投与間隔を変えたときの曝露を比較した。さらに、CD20受容体占有率（RO）を介したサイトカイン放出症候群（CRS）のモデルが、どこまで解釈できるかを調べた。
結果として、投与量や間隔によるPKの変化は公開モデルから計算できた。一方、CRSの発現率を変更レジメンの臨床リスクとして予測するための情報は不足していた。 公表値に近い計算結果が得られることと、患者集団への予測が検証されたことは区別する必要がある。
本記事は公開情報を使った研究・学習用の解析であり、査読を受けていない。患者の治療、臨床試験の用量決定、個人のCRSリスク判定には使用しない。以下の変更レジメンはモデルの性質を調べるための計算条件であり、投与の提案ではない。
背景：段階的増量をPKだけで説明できるか
モスネツズマブはCD20とCD3を標的とする二重特異性抗体である。米国の静注製剤LUNSUMIOの投与法では、Cycle 1のDay 1、8、15にそれぞれ1、2、60 mg、Cycle 2のDay 1に60 mgを投与し、Cycle 3以降は21日ごとに30 mgを投与する。Cycle 1は少なくとも4時間の点滴であり、その後は忍容性に応じて2時間とされる。本解析はこの静注製剤を対象とする。米国の製品情報
初期投与量が小さいと初回の血中濃度ピークが下がる。しかし、「血中濃度が下がる」ことから「CRSがどれだけ減る」かを直ちに計算できるわけではない。本記事では、投与入力からPK、PKからRO、ROからCRSという各段階で、公開情報が支える範囲を切り分ける。
| 段階 | 問いたいこと | この解析で扱うもの |
|---|---|---|
| 投与 → PK | 初回量・増量幅・間隔が曝露をどう変えるか | モデル内の定量比較 |
| PK → RO | 残存anti-CD20が結合競合にどう影響するか | 仮定を明示した感度解析 |
| RO → CRS | 変更レジメンのCRS発現率を予測できるか | 公表回帰式の探索的な出力。臨床予測としては未検証 |
QOIとCOU
Question of interest（QOI）は、初回投与量、2回目投与量、step-up間隔を変更したとき、初回ピーク、初回60 mg投与直前濃度、42日間の累積曝露がどれだけ変わるか、そしてその違いをRO・CRSに結び付けられるか、である。
Context of use（COU）は、公開モデルの再実装と研究・学習目的のレジメン比較である。臨床的に最適なレジメンを選ぶ解析ではない。
公開情報とモデルの再実装
使用した資料
構造モデルとNONMEMコード、ROの式はBenderらのPopPK論文と補足資料から取得した。最終推定値、個体間変動（IIV）、CRS回帰式はLiらのExposure–Response論文の補足資料でも確認した。Bender et al., 2024、Li et al., 2025
| 資料 | 主な用途 | 確認箇所 |
|---|---|---|
| PopPK論文の補足資料（cts13825） | 構造式、NONMEM実装、競合結合、既存患者に基づく曝露要約 | Equation 1、Table S2、最終モデルコード |
| E–R論文の補足資料（cpt3445） | 最終固定効果、Ω、RO–CRS回帰 | Tables S3、S12 |
| LUNSUMIO静注製剤の製品情報 | 投与日、点滴時間、Cycle別の曝露値 | Dosage and Administration、Clinical Pharmacology |
補足資料はEurope PMCの公開アーカイブから取得した。掲載するのは独自に作成した図、計算コードと結果であり、論文の図表の転載ではない。取得元とファイルのSHA-256を資料マニフェストに記録した。
time-varying clearanceを持つ2-compartment model
中心・末梢コンパートメント内の薬物量をそれぞれ \(A_1,A_2\)（mg）、中心濃度を \(C\)（mg/L = µg/mL）とする。
\[ C(t)=\frac{A_1(t)}{V_1} \]
\[ CL(t)=CL_{ss}+(CL_{base}-CL_{ss})\exp\left(-\frac{\ln 2}{HL_{trans}}t\right) \]
\[ \frac{dA_1}{dt}=R_{in}(t)-\frac{CL(t)+Q}{V_1}A_1+\frac{Q}{V_2}A_2 \]
\[ \frac{dA_2}{dt}=\frac{Q}{V_1}A_1-\frac{Q}{V_2}A_2 \]
\(R_{in}\)は点滴速度（mg/day）である。時刻は初回投与開始を0日とし、CLの時間変化はその時刻から連続して計算する。投与ごとにCLの時計をリセットしない。
| パラメータ | 典型値 | 公表95% CI | 単位 |
|---|---|---|---|
| CLbase | 1.08 | 0.962–1.20 | L/day |
| CLss | 0.584 | 0.561–0.607 | L/day |
| HLtrans | 16.3 | 14.026–18.6 | day |
| V1 | 5.49 | 5.221–5.76 | L |
| V2 | 6.17 | 5.729–6.61 | L |
| Q | 1.46 | 1.354–1.57 | L/day |
固定効果はE–R補足資料Table S3の最終推定値を使用した。NONMEMコードの初期値を最終推定値として読み込まないようにした。
このCLの減少は、元のPopPK解析の時間依存性を再現する経験的な表現である。今回の再実装では病態、腫瘍量や治療反応を時間とともに解く機構モデルにはしていない。したがって、変更レジメンでも同じ時間依存CLが成立するという仮定が含まれる。
共変量と仮想患者
基本計算は、男性、体重78 kg、アルブミン39 g/L、腫瘍SPD 2,970 mm²（平方根54.5 mm）の参照条件に固定した。患者の背景分布を推定して作った集団ではない。
IIVは、\(\eta\sim N(0,\Omega)\)、\(\theta_i=\theta_{typ}\exp(\eta_i)\)として実装した。パラメータの順序をCLbase、V1、CLss、HLtrans、V2とすると、使用した行列は次の通りである。QにはIIVを付けていない。
\[ \Omega=\begin{pmatrix} 0.426&0.1804&0&0&0\\ 0.1804&0.0981&0&0&0\\ 0&0&0.0343&-0.08925&0\\ 0&0&-0.08925&0.739&0\\ 0&0&0&0&0.0621 \end{pmatrix} \]
Table S3のCLbase–V1共分散には記号欄に不整合があるため、行名とNONMEMのランダム効果ブロックの順序も照合した。Ωの正定値性を確認し、乱数seedを20261002として2,000例を生成した。Ω自体の推定不確実性と測定誤差は加えていない。
解析条件と数値検証
参照レジメンは解析時刻0、7、14、21日に1、2、60、60 mg、その後42、63日に30 mgとした。通常の診療上のDay 1は解析時刻0日、Day 8は7日に対応する。Cycle 1のAUCは0–21日、Cycle 4のAUCは63–84日である。投与間隔末の濃度は次回投与直前の値を取り、42日AUCには42日開始の30 mg点滴を含めない。
| 比較条件 | 初回量 | 2回目量 | 最初の2間隔 | 初回60 mgの時刻 |
|---|---|---|---|---|
| 参照 | 1 mg | 2 mg | 7日・7日 | 14日 |
| 初回量変更 | 0.5 / 2 mg | 2 mg | 7日・7日 | 14日 |
| 2回目量変更 | 1 mg | 1 / 4 mg | 7日・7日 | 14日 |
| 間隔短縮 | 1 mg | 2 mg | 3.5日・3.5日 | 7日 |
| 間隔延長 | 1 mg | 2 mg | 14日・14日 | 28日 |
間隔変更では、最初の2間隔を両方変更し、その後の投与日も初回60 mgを基準に動かす。2回目の60 mgは初回60 mgの7日後、最初の30 mgは28日後である。この定義では、短縮条件は35日までに30 mgも投与されるが、参照条件ではその投与は42日になる。固定した0–42日のAUC比較には、このスケジュール差が入る。
ODEは点滴開始・終了時刻に整合した15分刻みの4次Runge–Kutta法で解き、AUCは台形則で積分した。刻みを7.5分にした計算と比較し、初回Cmax、60 mg直前濃度、AUCの差が0.01%未満であることを確認した。また、イベントごとに区間を分割した独立の適応刻みソルバー（SciPy solve_ivp）でも42日AUCを計算した。125.082680と125.082679 µg·day/mLが得られ、数値積分の差は結論に影響しない大きさだった。
この検証は実装と数値計算の検証である。新たな臨床データを使ったモデル検証ではない。
結果1：公表曝露要約との照合
図1．参照条件の濃度時間曲線。Cycle 1の点滴は4時間、以降は2時間として計算した。独自計算。
| 指標・単位 | 公表の幾何平均 | 典型値の再計算 | 差 (%) |
|---|---|---|---|
| Cycle 1 AUC：µg·day/mL | 35.2 | 34.80 | -1.14 |
| Cycle 1 Cmax：µg/mL | 11.1 | 10.66 | -4.00 |
| Cycle 1投与間隔末の濃度：µg/mL | 2.6 | 2.59 | -0.42 |
| AUC 0–42：µg·day/mL | 126 | 125.08 | -0.73 |
| Cycle 4 AUC：µg·day/mL | 52.9 | 52.70 | -0.38 |
参照値は製品情報とPopPK補足資料Table S2に基づく。投与間隔末の参照値には、2025年FDAラベルTable 10にも記載されている2.6 µg/mLを使用した。
AUCと投与間隔末濃度は近かったが、Cycle 1 Cmaxには約−4.0%の差があった。したがって「全指標が3%以内」とはしない。Cmaxは点滴時間にも依存するため、点滴を瞬時投与として扱うことは避けた。
さらに、公表値は患者集団のモデル予測曝露の幾何平均であり、ここで計算した参照共変量・η=0の典型曲線とは統計量が異なる。公表の患者別パラメータ、共変量分布と同じものを用いた再現ではない。この近接性は実装の整合性を確認する参考材料であり、独立した外部検証や患者集団の再現を意味しない。
結果2：同一仮想患者でのレジメン比較
各仮想患者のパラメータを固定し、参照レジメンと変更レジメンを両方与えた。患者ごとの変化率を計算し、その中央値を表に示した。これにより、比較群ごとに異なる患者を抽出したことによる差を除ける。
| 変更条件 | 初回Cmax | 初回60 mg直前濃度 | AUC 0–42 |
|---|---|---|---|
| 初回 1 → 0.5 mg | -50.0% | -11.0% | -0.4% |
| 初回 1 → 2 mg | +100.0% | +22.1% | +0.8% |
| 2回目 2 → 1 mg | +0.0% | -39.0% | -0.9% |
| 2回目 2 → 4 mg | +0.0% | +77.9% | +1.7% |
| 間隔 7 → 3.5日 | +0.0% | +54.1% | +18.6% |
| 間隔 7 → 14日 | +0.0% | -36.3% | -30.8% |
図2．点はpaired変化率の中央値、横線は仮想患者間の2.5–97.5パーセンタイル。推定値の95%信頼区間ではない。すべて同じ2,000例を使った。独自計算。
初回量と2回目量は、主として直近の曝露を動かす
このPKモデルは、同じ時間依存CLの下では投与入力に対して線形である。したがって、ほかの投与前に測る初回Cmaxは初回量に比例し、0.5 mgでは−50%、2 mgでは＋100%となる。この厳密な比例はモデルの性質であり、臨床的なCRSの比例関係ではない。
2回目量を1 mgに減らすと初回60 mg直前濃度の中央値は約−39%、4 mgに増やすと約＋78%だった。一方、初回量・2回目量だけを変えた場合の42日AUC変化は約−0.9%から＋1.7%である。累積曝露では、その後の60 mg投与の寄与が大きい。
間隔変更では累積曝露も変わる
間隔を3.5日に短縮すると60 mg直前濃度は約＋54%、14日に延長すると約−36%だった。ただし、同じ60 mg直前であっても、比較する時刻は7、14、28日と異なる。残存量、分布・消失、時間依存CLが合わさった結果である。
42日AUCは短縮で約＋19%、延長で約−31%となった。初回・2回目量の変更と異なり、小さな差ではない。大用量が入る時刻と、42日までに含まれる維持投与の有無も変わるため、「投与間隔がCLに及ぼす純粋な効果」とは解釈しない。
結果3：何の不確実性を調べたのか
固定効果の推定不確実性
6つの固定効果を一つずつ、公表95% CIの下限・上限に置き換えた。これはone-at-a-time感度解析であり、固定効果を同時に変える確率的解析ではない。固定効果推定値の共分散行列を持たないため、結果を曝露の95%信頼区間とは呼ばない。
図3．他のパラメータを固定した局所的な感度解析。各点の条件と数値はCSVで公開する。独自計算。
最大の絶対変化は、初回Cmaxで約4.9%、60 mg直前濃度で約7.9%、42日AUCで約3.9%だった。初回量変更によるピークや、2回目量・間隔変更による直前濃度の変化は、今回の局所的な固定効果感度より大きい。
ただし、これは全指標に共通する結論ではない。初回量・2回目量だけを変えた42日AUCの差は、このパラメータ感度より小さい場合がある。また、これらのCIでは構造モデルの誤指定や患者背景の不足を評価できない。
個体間変動
参照レジメンを仮想患者に与えた絶対曝露は次の通りだった。
| 指標 | 中央値 | 2.5–97.5パーセンタイル |
|---|---|---|
| 初回Cmax (µg/mL) | 0.175 | 0.096–0.333 |
| 初回60 mg直前濃度 (µg/mL) | 0.101 | 0.032–0.191 |
| AUC 0–42 (µg·day/mL) | 123.830 | 55.519–193.473 |
絶対曝露にはばらつきがある一方、paired解析では2回目量や間隔変更の効果の方向は維持された。しかし、この集団は参照共変量を固定してΩだけを反映したものであり、実患者の体重、腫瘍量、アルブミン、残存anti-CD20とIIVのjoint distributionを再現していない。患者集団の曝露分布として検証されたものではない。
結果4：PKからRO、ROからCRSへ
残存anti-CD20との結合競合
公開ROモデルを、すべてµg/mLの単位で次のように実装した。
\[ RO(t)=100\frac{C_M(t)}{C_M(t)+K_{D,M}+\frac{K_{D,M}}{K_{D,R}}C_R(t)+\frac{K_{D,M}}{K_{D,G}}C_G(t)} \]
ここでMはモスネツズマブ、Rはrituximab（RTX）、Gはobinutuzumab（OBI）である。\(K_{D,M}=10.2\)、\(K_{D,R}=0.675\)、\(K_{D,G}=0.600\) µg/mLとした。RTXとOBIは開始時濃度から、それぞれ半減期24日、28日で減衰させた。開始時濃度を42日間一定に保つ計算ではない。PopPK補足資料のRO式・実装
RO式は競合結合に基づくCD20占有の近似である。血中濃度を入力して計算したROを、腫瘍組織内の占有率やT細胞活性化の直接測定と同一視しない。
公表CRS式に入るのは42日間の最大RO
Study GO29781の再発・難治性非ホジキンリンパ腫に対する段階的増量群（Group B）について、ASTCT基準によるGrade ≥2 CRSの公表ロジスティック式を用いた。これは承認用量の濾胞性リンパ腫患者だけから構成されるモデルではない。
\[ \operatorname{logit}(p)=-2.64+0.0196\,RO_{max,0-42} \]
ROの入力は0–100の百分率であり、0–1の割合ではない。また、初回投与の直後の最大ROではなく、最初の42日間全体の最大ROである。E–R補足資料Table S12
Table S12の見出しはGrade 2と記載されているが、論文本文のGrade ≥2という定義に従った。回帰係数のSEは切片0.279、傾き0.00804と公表されているが、両者の共分散がないため、ここではCRS式出力の信頼区間を作っていない。
この違いは重要である。後半の60 mg投与が最大ROを決める場合、初回量を変更してもこの指標は大きく動かない可能性がある。しかし、それは初回投与のprimingや投与履歴がCRSに無関係という意味ではない。この回帰式には、変更した投与履歴を機構的に表現する状態変数がない。
anti-CD20の仮定を二通りに分ける
まずPKを固定して、競合結合だけを変えた。次に、公開PKモデルに含まれる開始時anti-CD20のCLbase共変量効果も反映した。公表実装の共変量項は、開始時anti-CD20をng/mLとして、\([\log(ACD20)/\log(500)]^{-0.573}\)である。
残存薬なしの結合条件では、共変量計算上は参照値500 ng/mLを使った。これはログに0を入れないための本解析上の参照条件であり、元のNONMEMの欠測補完を患者データなしで完全再現したものではない。コードの下限処理を、定量限界未満の患者一般への検証済みルールとは解釈しない。
| 開始時anti-CD20 | 最大RO：結合のみ | CRS式の出力：結合のみ | 最大RO：CLも変更 | CRS式の出力：CLも変更 |
|---|---|---|---|---|
| 残存薬なし | 56.63% | 17.80% | 56.63% | 17.80% |
| RTX 10 µg/mL | 12.60% | 8.37% | 12.81% | 8.40% |
| OBI 10 µg/mL | 10.71% | 8.09% | 10.89% | 8.12% |
| RTX 305 µg/mL | 0.53% | 6.73% | 0.54% | 6.73% |
図4．参照レジメン・典型患者での感度解析。点を結ぶ線は比較を見やすくするためであり、薬剤間の連続的な用量反応曲線ではない。CRS式の出力は個人のリスク予測ではない。独自計算。
この条件では、CLbase共変量の反映による差より、RO式の競合結合による差が大きかった。ただし、4つの濃度条件を同じ割合で集団に混ぜたり、記事の表をそのまま実患者の発現率としたりする根拠はない。
特に、ROが0に近づいても回帰式の出力は約6.7%となる。これは切片の性質であり、薬理学的な残余リスクを独立に実証した数字ではない。単純な式が出す確率には、元の解析集団と解析条件の影響が含まれる。
何が言えて、何が言えないか
| 判断対象 | この解析から支持できること | 支持できないこと |
|---|---|---|
| PK実装 | 構造式・投与イベントを独立に実装し、数値計算の整合性を確認した | 新しい患者データによる外部検証 |
| 投与量比較 | 同じモデル・同じ患者パラメータ内での曝露変化 | 変更後もモデルが臨床で同じ精度を持つ保証 |
| 不確実性 | 局所的な固定効果感度と、Ωに基づく個体差 | 固定効果・共変量・モデル構造を統合した完全な予測分布 |
| RO–CRS | 残存薬・回帰式の仮定に結果がどれだけ依存するか | 変更レジメンの集団CRS発現率、個人リスク、最適レジメン |
患者別のPK観測、共変量と残存anti-CD20濃度の同時分布、CRSの時点と投与履歴を結び付けたデータは今回の公開資料から得られなかった。また、元のRO–CRS関係が、仮想的な3.5日・14日間隔に輸送可能かを確認するデータもない。
公表値との近接性だけを根拠に「PKモデルをqualifyできた」と結論することは、この証拠だけではできない。本解析の到達点は、公開PKモデルの計算実装を検証し、そのモデル内でのレジメン曝露比較を行ったことである。RO–CRSへの翻訳は感度解析として探索できるが、患者集団の発現率予測としては検証できない。
再実行用コードと結果
Python 3.12で実行した。仮想患者の乱数、使用ライブラリ、数値積分条件を実行記録に保存した。
- 解析コード：run_analysis.py
- 依存ライブラリ：requirements.txt
- 実行方法と前提：README
- 公表値との照合CSV
- pairedレジメン比較CSV
- 仮想集団の参照曝露CSV
- 固定効果感度CSV
- anti-CD20感度CSV
python -m pip install -r requirements.txt
python run_analysis.py
図表は公表データを入力した本記事独自のシミュレーション結果であり、実測した患者データではない。
参考文献・資料
- Bender B, et al. Population pharmacokinetics and CD20 binding dynamics for mosunetuzumab in relapsed/refractory B-cell non-Hodgkin lymphoma. Clinical and Translational Science. 2024;17:e13825. DOI: 10.1111/cts.13825. 補足資料.
- Li CC, et al. A Novel Step-Up Dosage Regimen for Enhancing the Benefit-to-Risk Ratio of Mosunetuzumab in Relapsed or Refractory Follicular Lymphoma. Clinical Pharmacology & Therapeutics. 2025;117:465–474（オンライン公開2024年）. DOI: 10.1002/cpt.3445. 補足資料.
- LUNSUMIO (mosunetuzumab-axgb), intravenous prescribing information. DailyMed; FDA 2025 label.
資料確認日：2026年10月2日。解析・執筆：Daiki Kumakura。所属先の非公開データは使用していない。