モスネツズマブの公開PopPKモデルを再実装する：段階的増量の曝露比較とRO–CRS外挿の限界
On this page
要旨
モスネツズマブ（mosunetuzumab）の段階的増量について、公開された母集団薬物動態（population pharmacokinetics、PopPK）モデルをPythonで再実装し、初回投与量・2回目投与量・投与間隔を変えたときの曝露を比較した。さらに、CD20受容体占有率（receptor occupancy、RO）を介したサイトカイン放出症候群（cytokine release syndrome、CRS）のモデルが、どこまで解釈できるかを調べた。
結果として、投与量や間隔によるPKの変化は公開モデルから計算できた。一方、CRSの発現率を変更レジメンの臨床リスクとして予測するための情報は不足していた。 公表値に近い計算結果が得られることと、患者集団への予測が検証されたことは区別する必要がある。
本記事は公開情報を使った研究・学習用の解析である。患者の治療、臨床試験の用量決定、個人のCRSリスク判定には使用しない。以下の変更レジメンはモデルの性質を調べるための計算条件であり、投与の提案ではない。
背景と目的
モスネツズマブは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時間とされる。本解析はこの静注製剤を対象とする (1)。
初期投与量が小さいと初回の血中濃度ピークが下がる。しかし、「血中濃度が下がる」ことから「CRSがどれだけ減る」かを直ちに計算できるわけではない。本記事では、投与入力からPK、PKからRO、ROからCRSという各段階で、公開情報が支える範囲を切り分ける。
表1に以下の項目を示す。
表1．解析の段階と対象
| Step | Question | What this analysis handles |
|---|---|---|
| Dose → PK | How do the first dose, step-up size and interval change exposure? | Quantitative comparison within the model |
| PK → RO | How does residual anti-CD20 affect binding competition? | Sensitivity analysis with explicit assumptions |
| RO → CRS | Can CRS incidence for a modified regimen be predicted? | Exploratory output of a published regression; not validated as clinical prediction |
解析の問いと用途
Question of interest（QOI）は、初回投与量、2回目投与量、step-up間隔を変更したとき、初回ピーク、初回60 mg投与直前濃度、42日間の累積曝露がどれだけ変わるか、そしてその違いをRO・CRSに結び付けられるか、である。
Context of use（COU）は、公開モデルの再実装と研究・学習目的のレジメン比較である。臨床的に最適なレジメンを選ぶ解析ではない。
方法
使用した資料
構造モデルとNONMEMコード、ROの式はBenderらのPopPK論文と補足資料から取得した。最終推定値、個体間変動（interindividual variability、IIV）、CRS回帰式はLiらの曝露反応（exposure–response、E–R）論文の補足資料でも確認した (2, 3)。
表2に以下の項目を示す。
表2．使用資料と確認箇所
| Source | Main use | Location |
|---|---|---|
| Supplement of the PopPK paper (cts13825) | Structural equations, NONMEM implementation, competitive binding, exposure summaries based on actual patients | Equation 1, Table S2, final model code |
| Supplement of the E–R paper (cpt3445) | Final fixed effects, Ω, RO–CRS regression | Tables S3, S12 |
| LUNSUMIO intravenous product information | Dosing days, infusion durations, exposure values by cycle | Dosage and Administration, Clinical Pharmacology |
補足資料はEurope PMCの公開アーカイブから取得した。掲載するのは独自に作成した図、計算コードと結果であり、論文の図表の転載ではない。取得元とファイルのSHA-256を資料マニフェストに記録した。
時間依存クリアランスを持つ2コンパートメントモデル
中心・末梢コンパートメント内の薬物量をそれぞれ \(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の時計をリセットしない。
表3に以下の項目を示す。
表3．公開モデルのパラメータ
| Parameter | Typical value | Published 95% CI | Unit |
|---|---|---|---|
| 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 |
固定効果は原著補足資料のTable S3の最終推定値を使用した (3)。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点滴を含めない。
表4に以下の項目を示す。
表4．比較する投与条件
| Comparison | First dose | Second dose | First two intervals | Time of first 60 mg dose |
|---|---|---|---|---|
| Reference | 1 mg | 2 mg | 7 days, 7 days | Day 14 |
| Changed first dose | 0.5 / 2 mg | 2 mg | 7 days, 7 days | Day 14 |
| Changed second dose | 1 mg | 1 / 4 mg | 7 days, 7 days | Day 14 |
| Shorter interval | 1 mg | 2 mg | 3.5 days, 3.5 days | Day 7 |
| Longer interval | 1 mg | 2 mg | 14 days, 14 days | Day 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が得られ、数値積分の差は結論に影響しない大きさだった。
この検証は実装と数値計算の検証である。新たな臨床データを使ったモデル検証ではない。
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日で減衰させた。原著補足資料のRO式と実装に従い、開始時濃度を42日間一定に保つ計算ではない (2)。
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である（原著補足資料のTable S12） (3)。
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の欠測補完を患者データなしで完全再現したものではない。コードの下限処理を、定量限界未満の患者一般への検証済みルールとは解釈しない。
結果
公表曝露要約との照合
図1．参照条件の濃度時間曲線。Cycle 1の点滴は4時間、以降は2時間として計算した。独自計算。
表5に以下の項目を示す。
表5．公表曝露要約と典型値による計算
| Metric and unit | Published geometric mean | Recalculated typical value | Difference (%) |
|---|---|---|---|
| Cycle 1 AUC: µg·day/mL | 35.2 | 34.80 | -1.14 |
| Cycle 1 Cmax: µg/mL | 11.1 | 10.66 | -4.00 |
| Concentration at end of Cycle 1 interval: µ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を使用した (1, 2, 4)。
AUCと投与間隔末濃度は近かったが、Cycle 1 Cmaxには約−4.0%の差があった。Cmaxは点滴時間にも依存するため、点滴を瞬時投与として扱うことは避けた。
さらに、公表値は患者集団のモデル予測曝露の幾何平均であり、ここで計算した参照共変量・η=0の典型曲線とは統計量が異なる。公表の患者別パラメータ、共変量分布と同じものを用いた再現ではない。この近接性は実装の整合性を確認する参考材料であり、独立した外部検証や患者集団の再現を意味しない。
同一仮想患者でのレジメン比較
各仮想患者のパラメータを固定し、参照レジメンと変更レジメンを両方与えた。患者ごとの変化率を計算し、その中央値を表に示した。これにより、比較群ごとに異なる患者を抽出したことによる差を除ける。
表6に以下の項目を示す。
表6．同一仮想患者における曝露の変化率
| Change | First Cmax | Concentration before first 60 mg | AUC 0–42 |
|---|---|---|---|
| First dose 1 → 0.5 mg | -50.0% | -11.0% | -0.4% |
| First dose 1 → 2 mg | +100.0% | +22.1% | +0.8% |
| Second dose 2 → 1 mg | +0.0% | -39.0% | -0.9% |
| Second dose 2 → 4 mg | +0.0% | +77.9% | +1.7% |
| Interval 7 → 3.5 days | +0.0% | +54.1% | +18.6% |
| Interval 7 → 14 days | +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に及ぼす純粋な効果」とは解釈しない。
固定効果の感度と個体間変動
固定効果の推定不確実性
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では構造モデルの誤指定や患者背景の不足を評価できない。
個体間変動
参照レジメンを仮想患者に与えた絶対曝露は次の通りだった。
表7に以下の項目を示す。
表7．仮想患者集団の絶対曝露
| Metric | Median | 2.5th–97.5th percentile |
|---|---|---|
| First Cmax (µg/mL) | 0.175 | 0.096–0.333 |
| Concentration before first 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を再現していない。患者集団の曝露分布として検証されたものではない。
PK・RO・CRSの接続
ROの競合結合のみを変える条件と、CLbase共変量も変える条件を比較した。結果を表8と図4に示す。
表8．残存anti-CD20の条件別感度
| Baseline anti-CD20 | Max RO: binding only | CRS equation output: binding only | Max RO: CL also changed | CRS equation output: CL also changed |
|---|---|---|---|---|
| No residual drug | 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%となる。これは切片の性質であり、薬理学的な残余リスクを独立に実証した数字ではない。単純な式が出す確率には、元の解析集団と解析条件の影響が含まれる。
考察
表9．解析が支持する判断と限界
| Subject of judgment | What this analysis supports | What it does not support |
|---|---|---|
| PK implementation | Independently implemented the structural equations and dosing events and confirmed numerical consistency | External validation with new patient data |
| Dose comparison | Exposure changes within the same model and the same patient parameters | A guarantee that the model has the same clinical accuracy after the change |
| Uncertainty | Local fixed-effect sensitivity and individual differences based on Ω | A full predictive distribution integrating fixed effects, covariates and model structure |
| RO–CRS | How much results depend on assumptions about residual drug and the regression | Population CRS incidence for a modified regimen, individual risk, the optimal regimen |
患者別の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
図表は公表データを入力した本記事独自のシミュレーション結果であり、実測した患者データではない。
参照
- Genentech, Inc. LUNSUMIO (mosunetuzumab-axgb) injection, for intravenous use: prescribing information [Internet]. DailyMed [cited 2026 Oct 2]. Available from: DailyMed.
- Bender B, Li CC, Marchand M, Turner DC, Li F, Vadhavkar S, et al. Population pharmacokinetics and CD20 binding dynamics for mosunetuzumab in relapsed/refractory B-cell non-Hodgkin lymphoma. Clin Transl Sci. 2024;17(6):e13825. doi:10.1111/cts.13825. Supplementary materials.
- Li CC, Bender B, Wilkins J, Li F, Turner DC, Wang B, et al. A novel step-up dosage regimen for enhancing the benefit-to-risk ratio of mosunetuzumab in relapsed or refractory follicular lymphoma. Clin Pharmacol Ther. 2025;117(2):465–474. Published online 2024 Sep 27. doi:10.1002/cpt.3445. Supplementary materials.
- Genentech, Inc. LUNSUMIO (mosunetuzumab-axgb) injection, for intravenous use: prescribing information [Internet]. US Food and Drug Administration; 2025 [cited 2026 Oct 2]. Available from: FDA label.