(57)【要約】【課題】アキシチニブの薬物動態を簡便に判定し、アキシニチブによる治療効果を予測する。【解決手段】被験者に関する特定の遺伝子多型と背景因子とを用いてアキシチニブの予測薬物動態パラメータを算出する工程とを含む、アキシチニブの薬物動態判定方法。【選択図】図１
2018014938.tif 000002
【特許請求の範囲】
【請求項１】
被験者に関する下記[1]〜[7]の情報を取得する工程と、
[1]  UGT1A1遺伝子の多型であって、変異型アレルが機能欠損型変異である多型
[2]  UGT1A7遺伝子の多型であって、変異型アレルが機能欠損型変異である多型
[3]  UGT1A9遺伝子の多型であって、変異型アレルが機能獲得型変異である多型
[4]  BCRP遺伝子の多型であって、変異型アレルが機能欠損型変異である多型
[5]  MDR1遺伝子の多型であって、変異型アレルが機能欠損型変異である多型
[6]  アキシチニブ投与量
[7]  前治療有無
上記[1]〜[7]の情報に基づいてアキシチニブの予測薬物動態パラメータを算出する工程とを含む、アキシチニブの薬物動態判定方法。
【請求項２】
上記[1]はUGT1A1＊6及び/又はUGT1A1＊28又はこれらに対して連鎖不平衡となる多型であることを特徴とする請求項１記載のアキシチニブの薬物動態判定方法。
【請求項３】
上記[2]はUGT1A7＊12又はこれに対して連鎖不平衡となる多型であることを特徴とする請求項１記載のアキシチニブの薬物動態判定方法。
【請求項４】
上記[3]はUGT1A9＊1b又はこれに対して連鎖不平衡となる多型であることを特徴とする請求項１記載のアキシチニブの薬物動態判定方法。
【請求項５】
上記[4]はBCRP遺伝子におけるC421A多型又はこれに対して連鎖不平衡となる多型であることを特徴とする請求項１記載のアキシチニブの薬物動態判定方法。
【請求項６】
上記[5]はMDR1遺伝子におけるT1236C多型又はこれに対して連鎖不平衡となる多型であることを特徴とする請求項１記載のアキシチニブの薬物動態判定方法。
【請求項７】
上記予測薬物動態パラメータは、全身クリアランスの予測値（予測CLtot値）であることを特徴とする請求項１記載のアキシチニブの薬物動態判定方法。
【請求項８】
上記全身クリアランスの予測値（予測CLtot値）を算出する際に、上記[1]の多型を有する場合に予測CLtot値が高くなり、上記[2]の多型を有する場合に予測CLtot値が低くなり、上記[3]の多型を有する場合に予測CLtot値が高くなり、上記[4]の多型を有する場合に予測CLtot値が低くなり、上記[5]の多型を有する場合に予測CLtot値が低くなる予測式を使用することを特徴とする請求項７記載のアキシチニブの薬物動態判定方法。
【請求項９】
上記予測式は、上記[1]の多型を有する場合に算出される予測CLtot値への影響が最も高く、次に上記[2]の多型を有する場合に算出される予測CLtot値への影響が高く、次に上記[7]の前治療が有る場合に算出される予測CLtot値への影響が高くなる式であることを特徴とする請求項８記載のアキシチニブの薬物動態判定方法。
【請求項１０】
上記請求項１〜９いずれか一項記載のアキシチニブの薬物動態判定方法により算出したアキシチニブの予測薬物動態パラメータに基づいて、アキシチニブによる抗腫瘍活性及び/又は副作用を判定する工程を含む、アキシニチブによる治療効果予測方法。
【請求項１１】
アキシチニブによる抗腫瘍活性は、腫瘍縮小効果であることを特徴とする請求項１０記載のアキシニチブによる治療効果予測方法。
【請求項１２】
アキシチニブによる副作用は、タンパク尿であることを特徴とする請求項１０記載のアキシニチブによる治療効果予測方法。
【請求項１３】
被験者に関する下記[1]〜[7]の情報を入力する入力部と、
[1]  UGT1A1遺伝子の多型であって、変異型アレルが機能欠損型変異である多型
[2]  UGT1A7遺伝子の多型であって、変異型アレルが機能欠損型変異である多型
[3]  UGT1A9遺伝子の多型であって、変異型アレルが機能獲得型変異である多型
[4]  BCRP遺伝子の多型であって、変異型アレルが機能欠損型変異である多型
[5]  MDR1遺伝子の多型であって、変異型アレルが機能欠損型変異である多型
[6]  アキシチニブ投与量
[7]  前治療有無
上記入力部で入力した[1]〜[7]の情報に基づいてアキシチニブの予測薬物動態パラメータを算出する演算部とを有する、アキシチニブの薬物動態判定装置。
【請求項１４】
上記[1]はUGT1A1＊6及び/又はUGT1A1＊28又はこれらに対して連鎖不平衡となる多型であることを特徴とする請求項１３記載のアキシチニブの薬物動態判定装置。
【請求項１５】
上記[2]はUGT1A7＊12又はこれに対して連鎖不平衡となる多型であることを特徴とする請求項１３記載のアキシチニブの薬物動態判定装置。
【請求項１６】
上記[3]はUGT1A9＊1b又はこれに対して連鎖不平衡となる多型であることを特徴とする請求項１３記載のアキシチニブの薬物動態判定装置。
【請求項１７】
上記[4]はBCRP遺伝子におけるC421A多型又はこれに対して連鎖不平衡となる多型であることを特徴とする請求項１３記載のアキシチニブの薬物動態判定装置。
【請求項１８】
上記[5]はMDR1遺伝子におけるT1236C多型又はこれに対して連鎖不平衡となる多型であることを特徴とする請求項１３記載のアキシチニブの薬物動態判定装置。
【請求項１９】
上記演算部は、上記予測薬物動態パラメータとして全身クリアランスの予測値（予測CLtot値）を算出することを特徴とする請求項１３記載のアキシチニブの薬物動態判定装置。
【請求項２０】
上記演算部は、上記[1]の多型を有する場合に予測CLtot値が高くなり、上記[2]の多型を有する場合に予測CLtot値が低くなり、上記[3]の多型を有する場合に予測CLtot値が高くなり、上記[4]の多型を有する場合に予測CLtot値が低くなり、上記[5]の多型を有する場合に予測CLtot値が低くなる予測式を使用して上記全身クリアランスの予測値（予測CLtot値）を算出することを特徴とする請求項１９記載のアキシチニブの薬物動態判定装置。
【請求項２１】
上記予測式は、上記[1]の多型を有する場合に算出される予測CLtot値への影響が最も高く、次に上記[2]の多型を有する場合に算出される予測CLtot値への影響が高く、次に上記[7]の前治療が有る場合に算出される予測CLtot値への影響が高くなる式であることを特徴とする請求項２０記載のアキシチニブの薬物動態判定装置。
【請求項２２】
被験者に関する下記[1]〜[7]の情報を入力する入力部と、
[1]  UGT1A1遺伝子の多型であって、変異型アレルが機能欠損型変異である多型
[2]  UGT1A7遺伝子の多型であって、変異型アレルが機能欠損型変異である多型
[3]  UGT1A9遺伝子の多型であって、変異型アレルが機能獲得型変異である多型
[4]  BCRP遺伝子の多型であって、変異型アレルが機能欠損型変異である多型
[5]  MDR1遺伝子の多型であって、変異型アレルが機能欠損型変異である多型
[6]  アキシチニブ投与量
[7]  前治療有無
上記入力部で入力した[1]〜[7]の情報に基づいてアキシチニブの予測薬物動態パラメータを算出し、当該予測薬物動態パラメータに基づいて、アキシチニブによる抗腫瘍活性及び/又は副作用を判定する演算部とを有する、アキシチニブの治療効果予測装置。
【請求項２３】
上記[1]はUGT1A1＊6及び/又はUGT1A1＊28又はこれらに対して連鎖不平衡となる多型であることを特徴とする請求項２２記載のアキシチニブの治療効果予測装置。
【請求項２４】
上記[2]はUGT1A7＊12又はこれに対して連鎖不平衡となる多型であることを特徴とする請求項２２記載のアキシチニブの治療効果予測装置。
【請求項２５】
上記[3]はUGT1A9＊1b又はこれに対して連鎖不平衡となる多型であることを特徴とする請求項２２記載のアキシチニブの治療効果予測装置。
【請求項２６】
上記[4]はBCRP遺伝子におけるC421A多型又はこれに対して連鎖不平衡となる多型であることを特徴とする請求項２２記載のアキシチニブの治療効果予測装置。
【請求項２７】
上記[5]はMDR1遺伝子におけるT1236C多型又はこれに対して連鎖不平衡となる多型であることを特徴とする請求項２２記載のアキシチニブの治療効果予測装置。
【請求項２８】
上記演算部は、上記予測薬物動態パラメータとして全身クリアランスの予測値（予測CLtot値）を算出することを特徴とする請求項２２記載のアキシチニブの治療効果予測装置。
【請求項２９】
上記演算部は、上記[1]の多型を有する場合に予測CLtot値が高くなり、上記[2]の多型を有する場合に予測CLtot値が低くなり、上記[3]の多型を有する場合に予測CLtot値が高くなり、上記[4]の多型を有する場合に予測CLtot値が低くなり、上記[5]の多型を有する場合に予測CLtot値が低くなる予測式を使用して上記全身クリアランスの予測値（予測CLtot値）を算出することを特徴とする請求項２８記載のアキシチニブの治療効果予測装置。
【請求項３０】
上記予測式は、上記[1]の多型を有する場合に算出される予測CLtot値への影響が最も高く、次に上記[2]の多型を有する場合に算出される予測CLtot値への影響が高く、次に上記[7]の前治療が有る場合に算出される予測CLtot値への影響が高くなる式であることを特徴とする請求項２９記載のアキシチニブの治療効果予測装置。
【請求項３１】
上記演算部は、アキシチニブによる抗腫瘍活性として腫瘍縮小効果を判定することを特徴とする請求項２２記載のアキシニチブによる治療効果予測装置。
【請求項３２】
上記演算部は、アキシチニブによる副作用としてタンパク尿を判定することを特徴とする請求項２２記載のアキシニチブによる治療効果予測装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、根治切除不能又は転移性の腎細胞癌等に適用可能なアキシチニブの薬物動態判定方法及びアキシチニブの薬物動態に基づくアキシニチブによる治療効果予測方法に関する。
【背景技術】
【０００２】
アキシチニブは、血管内皮増殖因子受容体（Vascular Endothelial Growth Factor Receptor: VEGFR）-1、-2、-3を選択的に阻害する経口チロシンキナーゼ阻害剤である。なお、VEGFR-1、-2、-3は、腫瘍増殖、癌の浸潤（腫瘍伸展）及び転移に関与していると考えられている。アキシチニブについては、これらVEGFR-1、-2、-3を選択的に阻害することで血管新生及びリンパ管新生を抑制し、腫瘍の増殖及び転移を抑制し、抗腫瘍活性を示すと考えられている。
【０００３】
現在、アキシチニブは、根治切除不能又は転移性の腎細胞癌に対する分子標的薬として使用されている。アキシチニブは、このような進行性腎細胞癌に対する治療効果の高い薬剤であるが、効果（腫瘍縮小率）や副作用の予測、及び至適投与量の設定が困難である。臨床においては10mg/日経口投与が基準開始量であるが、下痢脱水、血圧上昇などの有害事象の有無に応じて2〜14mg/日の範囲で投与量を変動する場合がある。臨床的には、投与量を減量しながらも腫瘍が完全消失した症例もあれば、投与量を増大しても血圧の変動なく治療効果も見られない例、基準量で投与していても突然脳出血をきたした症例などが報告されている。
【０００４】
アキシチニブについては、その血中濃度（薬物動態解析パラメータであるAUC: area under the curve）と治療効果及び有害事象とが相関することが報告されている（非特許文献１）。ところが、AUC等の薬物動態解析パラメータは、簡便に測定するためのキットやシステムがなく、臨床の現場において測定が非常に困難である。また、仮に測定システムが有ったとしても、最低6〜7回の採血が必要な薬物動態解析をすべての投与患者に対して実施することは実臨床では極めて困難である。
【０００５】
薬物の血中濃度を規定する要素として、一般的には、[1]投与量、投与速度及び投与回数、[2]患者体表面積、年齢及び性別、[3]体内の代謝速度の3要素に大別できる。なかでも[3]を規定する要素として、代謝酵素の発現量が挙げられる。たとえば、代謝酵素のなかには遺伝多型により発現量が大きく異なり、その結果、代謝速度が大きく異なることが知られている。アキシチニブの代謝酵素としては、CYP3A4及びUGT1Aが報告されている（非特許文献２）。従来、これらの代謝酵素にも遺伝多型が報告されている。たとえば、遺伝子多型UGT1A*28は、プロモーター領域内のTATA boxに関する多型であり、TAの繰り返し配列（通常6回）が7回になると薬剤代謝速度が低下し、薬剤の体内蓄積による毒性が増加することが、イリノテカンなどの薬剤で証明されている。
【０００６】
しかし、アキシチニブについては、上述した代謝酵素の遺伝子多型と、血中濃度及び毒性との関係は不明であり、また、イリノテカンに関する上記知見を当てはめられる根拠も無かった。
【０００７】
ところで、特許文献１には、腫瘍の組織試料におけるＣＤ６８発現レベルに基づいて、当該腫瘍のアキシチニブ感受性を評価する方法が開示されている。また、アキシチニブとは異なる抗癌剤であるが分子標的薬エルロチニブについては、例えば特許文献２に開示されるように、ヒトABCB1（Homo sapiens ATP-binding cassette, sub-family B (MDR/TAP), member 1）（P-糖蛋白質)遺伝子の多型に基づいてエルロチニブの副作用又は薬効の程度を判定する方法が知られている。詳細には、特許文献２には、ABCB1遺伝子の所定の遺伝子型が、エルロチニブの血中濃度と毒性をダイレクトに予測する指標となることが開示されている。
【０００８】
ところが、特許文献２に開示された技術は、上皮増殖因子受容体チロシンキナーゼ阻害薬であるエルロチニブに関する技術であり、また、エルロニチブの治療効果を直接予測するモデルではないこと、複雑な薬物動態予測をABCB1多型のみで行っていることなどから、アキシチニブに適用できる技術ではない。
【先行技術文献】
【特許文献】
【０００９】
【特許文献１】特開2015-210268号公報
【特許文献２】特開2010-263810号公報
【非特許文献】
【００１０】
【非特許文献１】Expert Opin. Drug Metab. Toxicol.(2012)8(2)
【非特許文献２】J. Clin Pharmacol.(2013)53:491-504
【発明の概要】
【発明が解決しようとする課題】
【００１１】
そこで、本発明は、簡便且つ高精度にアキシチニブの薬物動態を判定し、また、判定したアキシチニブの薬物動態に基づいてアキシニチブによる治療効果を予測する方法を提供するとともに、これら方法を適用した判定装置を提供することを目的とする。
【課題を解決するための手段】
【００１２】
本発明者らは、上記課題を解決すべく鋭意検討を行った結果、被験者の特定の多型に関する遺伝子型に基づいてアキシチニブの薬物動態を高精度に予測できること、また、予測したアキシチニブの薬物動態に基づいてアキシチニブの治療効果を高精度に予測できることを見出し、本発明を完成するに至った。本発明は以下を包含する。
【００１３】
（１）被験者に関する下記[1]〜[7]の情報を取得する工程と、
[1]  UGT1A1遺伝子の多型であって、変異型アレルが機能欠損型変異である多型
[2]  UGT1A7遺伝子の多型であって、変異型アレルが機能欠損型変異である多型
[3]  UGT1A9遺伝子の多型であって、変異型アレルが機能獲得型変異である多型
[4]  BCRP遺伝子の多型であって、変異型アレルが機能欠損型変異である多型
[5]  MDR1遺伝子の多型であって、変異型アレルが機能欠損型変異である多型
[6]  アキシチニブ投与量
[7]  前治療有無
上記[1]〜[7]の情報に基づいてアキシチニブの予測薬物動態パラメータを算出する工程とを含む、アキシチニブの薬物動態判定方法。
（２）上記[1]はUGT1A1＊6及び/又はUGT1A1＊28又はこれらに対して連鎖不平衡となる多型であることを特徴とする（１）記載のアキシチニブの薬物動態判定方法。
（３）上記[2]はUGT1A7＊12又はこれに対して連鎖不平衡となる多型であることを特徴とする（１）記載のアキシチニブの薬物動態判定方法。
（４）上記[3]はUGT1A9＊1b又はこれに対して連鎖不平衡となる多型であることを特徴とする（１）記載のアキシチニブの薬物動態判定方法。
（５）上記[4]はBCRP遺伝子におけるC421A多型又はこれに対して連鎖不平衡となる多型であることを特徴とする（１）記載のアキシチニブの薬物動態判定方法。
（６）上記[5]はMDR1遺伝子におけるT1236C多型又はこれに対して連鎖不平衡となる多型であることを特徴とする（１）記載のアキシチニブの薬物動態判定方法。
（７）上記予測薬物動態パラメータは、全身クリアランスの予測値（予測CLtot値）であることを特徴とする（１）記載のアキシチニブの薬物動態判定方法。
（８）上記全身クリアランスの予測値（予測CLtot値）を算出する際に、上記[1]の多型を有する場合に予測CLtot値が高くなり、上記[2]の多型を有する場合に予測CLtot値が低くなり、上記[3]の多型を有する場合に予測CLtot値が高くなり、上記[4]の多型を有する場合に予測CLtot値が低くなり、上記[5]の多型を有する場合に予測CLtot値が低くなる予測式を使用することを特徴とする（７）記載のアキシチニブの薬物動態判定方法。
（９）上記予測式は、上記[1]の多型を有する場合に算出される予測CLtot値への影響が最も高く、次に上記[2]の多型を有する場合に算出される予測CLtot値への影響が高く、次に上記[7]の前治療が有る場合に算出される予測CLtot値への影響が高くなる式であることを特徴とする（８）記載のアキシチニブの薬物動態判定方法。
（１０）上記（１）乃至（９）いずれかに記載のアキシチニブの薬物動態判定方法により算出したアキシチニブの予測薬物動態パラメータに基づいて、アキシチニブによる抗腫瘍活性及び/又は副作用を判定する工程を含む、アキシニチブによる治療効果予測方法。
（１１）アキシチニブによる抗腫瘍活性は、腫瘍縮小効果であることを特徴とする（１０）記載のアキシニチブによる治療効果予測方法。
（１２）アキシチニブによる副作用は、タンパク尿であることを特徴とする（１０）記載のアキシニチブによる治療効果予測方法。
（１３）被験者に関する下記[1]〜[7]の情報を入力する入力部と、
[1]  UGT1A1遺伝子の多型であって、変異型アレルが機能欠損型変異である多型
[2]  UGT1A7遺伝子の多型であって、変異型アレルが機能欠損型変異である多型
[3]  UGT1A9遺伝子の多型であって、変異型アレルが機能獲得型変異である多型
[4]  BCRP遺伝子の多型であって、変異型アレルが機能欠損型変異である多型
[5]  MDR1遺伝子の多型であって、変異型アレルが機能欠損型変異である多型
[6]  アキシチニブ投与量
[7]  前治療有無
上記入力部で入力した[1]〜[7]の情報に基づいてアキシチニブの予測薬物動態パラメータを算出する演算部とを有する、アキシチニブの薬物動態判定装置。
（１４）上記[1]はUGT1A1＊6及び/又はUGT1A1＊28又はこれらに対して連鎖不平衡となる多型であることを特徴とする（１３）記載のアキシチニブの薬物動態判定装置。
（１５）上記[2]はUGT1A7＊12又はこれに対して連鎖不平衡となる多型であることを特徴とする（１３）記載のアキシチニブの薬物動態判定装置。
（１６）上記[3]はUGT1A9＊1b又はこれに対して連鎖不平衡となる多型であることを特徴とする（１３）記載のアキシチニブの薬物動態判定装置。
（１７）上記[4]はBCRP遺伝子におけるC421A多型又はこれに対して連鎖不平衡となる多型であることを特徴とする（１３）記載のアキシチニブの薬物動態判定装置。
（１８）上記[5]はMDR1遺伝子におけるT1236C多型又はこれに対して連鎖不平衡となる多型であることを特徴とする（１３）記載のアキシチニブの薬物動態判定装置。
（１９）上記演算部は、上記予測薬物動態パラメータとして全身クリアランスの予測値（予測CLtot値）を算出することを特徴とする（１３）記載のアキシチニブの薬物動態判定装置。
（２０）上記演算部は、上記[1]の多型を有する場合に予測CLtot値が高くなり、上記[2]の多型を有する場合に予測CLtot値が低くなり、上記[3]の多型を有する場合に予測CLtot値が高くなり、上記[4]の多型を有する場合に予測CLtot値が低くなり、上記[5]の多型を有する場合に予測CLtot値が低くなる予測式を使用して上記全身クリアランスの予測値（予測CLtot値）を算出することを特徴とする（１９）記載のアキシチニブの薬物動態判定装置。
（２１）上記予測式は、上記[1]の多型を有する場合に算出される予測CLtot値への影響が最も高く、次に上記[2]の多型を有する場合に算出される予測CLtot値への影響が高く、次に上記[7]の前治療が有る場合に算出される予測CLtot値への影響が高くなる式であることを特徴とする（２０）記載のアキシチニブの薬物動態判定装置。
（２２）被験者に関する下記[1]〜[7]の情報を入力する入力部と、
[1]  UGT1A1遺伝子の多型であって、変異型アレルが機能欠損型変異である多型
[2]  UGT1A7遺伝子の多型であって、変異型アレルが機能欠損型変異である多型
[3]  UGT1A9遺伝子の多型であって、変異型アレルが機能獲得型変異である多型
[4]  BCRP遺伝子の多型であって、変異型アレルが機能欠損型変異である多型
[5]  MDR1遺伝子の多型であって、変異型アレルが機能欠損型変異である多型
[6]  アキシチニブ投与量
[7]  前治療有無
上記入力部で入力した[1]〜[7]の情報に基づいてアキシチニブの予測薬物動態パラメータを算出し、当該予測薬物動態パラメータに基づいて、アキシチニブによる抗腫瘍活性及び/又は副作用を判定する演算部とを有する、アキシチニブの治療効果予測装置。
（２３）上記[1]はUGT1A1＊6及び/又はUGT1A1＊28又はこれらに対して連鎖不平衡となる多型であることを特徴とする（２２）記載のアキシチニブの治療効果予測装置。
（２４）上記[2]はUGT1A7＊12又はこれに対して連鎖不平衡となる多型であることを特徴とする（２２）記載のアキシチニブの治療効果予測装置。
（２５）上記[3]はUGT1A9＊1b又はこれに対して連鎖不平衡となる多型であることを特徴とする（２２）記載のアキシチニブの治療効果予測装置。
（２６）上記[4]はBCRP遺伝子におけるC421A多型又はこれに対して連鎖不平衡となる多型であることを特徴とする（２２）記載のアキシチニブの治療効果予測装置。
（２７）上記[5]はMDR1遺伝子におけるT1236C多型又はこれに対して連鎖不平衡となる多型であることを特徴とする（２２）記載のアキシチニブの治療効果予測装置。
（２８）上記演算部は、上記予測薬物動態パラメータとして全身クリアランスの予測値（予測CLtot値）を算出することを特徴とする（２２）記載のアキシチニブの治療効果予測装置。
（２９）上記演算部は、上記[1]の多型を有する場合に予測CLtot値が高くなり、上記[2]の多型を有する場合に予測CLtot値が低くなり、上記[3]の多型を有する場合に予測CLtot値が高くなり、上記[4]の多型を有する場合に予測CLtot値が低くなり、上記[5]の多型を有する場合に予測CLtot値が低くなる予測式を使用して上記全身クリアランスの予測値（予測CLtot値）を算出することを特徴とする（２８）記載のアキシチニブの治療効果予測装置。
（３０）上記予測式は、上記[1]の多型を有する場合に算出される予測CLtot値への影響が最も高く、次に上記[2]の多型を有する場合に算出される予測CLtot値への影響が高く、次に上記[7]の前治療が有る場合に算出される予測CLtot値への影響が高くなる式であることを特徴とする（２９）記載のアキシチニブの治療効果予測装置。
（３１）上記演算部は、アキシチニブによる抗腫瘍活性として腫瘍縮小効果を判定することを特徴とする（２２）記載のアキシニチブによる治療効果予測装置。
（３２）上記演算部は、アキシチニブによる副作用としてタンパク尿を判定することを特徴とする（２２）記載のアキシニチブによる治療効果予測装置。
【発明の効果】
【００１４】
本発明に係るアキシチニブの薬物動態判定方法及び装置によれば、アキシチニブの薬物動態パラメータを非常に簡便な方法により高精度に予測することができる。
【００１５】
そして、本発明に係るアキシチニブによる治療効果判定方法及び装置によれば、予測したアキシチニブの薬物動態パラメータに基づいて、アキシチニブによる抗腫瘍効果及び/又は副作用を高精度に予測することができる。
【００１６】
したがって、本発明によれば、アキシチニブを必要とする患者に対して、アキシチニブの適切な投与量を決定するため有効な支援情報を提供することができる。
【図面の簡単な説明】
【００１７】
【図１】本発明を適用したアキシチニブの薬物動態判定装置を模式的に示す概略構成図である。
【図２】本発明を適用したアキシチニブによる治療効果判定装置を模式的に示す概略構成図である。
【図３】実施例で使用した40例について、所定の遺伝子多型の遺伝子型解析結果を一覧として示す図である。
【図４】実施例で使用した40例について、患者の背景因子、腫瘍縮小率、副作用の有無、及び薬物動態パラメータを解析した結果を一覧として示す図である。
【図５】実施例で使用した40例について、遺伝子多型の情報、患者の背景因子及び実測CLtot値を用いて、解析変数減少法（Backward Selection）により決定係数を計算した結果の一部を示す図である。
【図６】実施例１で特定した8因子を用いて算出した予測CLtot値と実測CLtot値とを比較した結果を示す特性図である。
【図７】実測CLtot値と腫瘍の最大縮小率とを比較した結果を示す特性図である。
【図８】実施例１で特定した8因子を用いて算出した予測CLtot値と腫瘍の最大縮小率とを比較した結果を示す特性図である。
【図９】実施例１で特定した8因子を用いて算出した予測CLtot値及び実測CLtot値とタンパク尿のグレードとを比較した結果を示す特性図である。
【発明を実施するための形態】
【００１８】
以下、本発明を詳細に説明する。
【００１９】
[概要]
本発明では、アキシチニブを投与する患者における所定の背景因子と当該患者における所定の遺伝子多型とに基づいて、アキシチニブの薬物動態を予測する。具体的には、アキシチニブに関する薬物動態パラメータを予測値（予測薬物動態パラメータ）として算出する。算出された予測薬物動態パラメータは、アキシチニブによる治療効果の判定に利用することができる。予測薬物動態パラメータに基づくアキシチニブによる治療効果を判定することによって、患者毎にアキシチニブの適切な投与量を決定することができる。
【００２０】
背景因子及び遺伝子多型については、それぞれ問診及び遺伝子型解析により情報として患者毎に取得することができる。或いは、背景因子については、過去の治療記録等に基づいて取得することもできるため、本発明を適用してアキシチニブの薬物動態を予測するに際して、必ずしも問診を要するものではない。また、遺伝子多型については、患者から採取した生体サンプル（血液、細胞等）を用いた各種の解析方法で取得することができる。或いは、遺伝子多型については、他の診察や治療に際して行った遺伝子型解析の結果を利用することもできるし、個人の遺伝子型データを格納したデータベース等があれば当該データベースから読み出して利用することもできる。したがって、遺伝子多型についても、本発明を適用してアキシチニブの薬物動態を予測するに際して、生体サンプルの採取及び遺伝子型解析を行う必要はない。
【００２１】
また、所定の背景因子と当該患者における所定の遺伝子多型は、予測薬物動態パラメータを目的変数としたとき、共に説明変数（従属変数）ということができる。言い換えると、所定の背景因子と当該患者における所定の遺伝子多型を説明変数とした予測式によって、目的変数である予測薬物動態パラメータを算出することができる。
【００２２】
ここで、アキシチニブとは、置換インダゾール誘導体であり、血管内皮増殖因子受容体VEGFR-1, -2及び-3を阻害するチロシンキナーゼ阻害剤（TKI）である。アキシチニブの分子式はC22H18N4OSであり、分子量が386.47であり、化学名がN-メチル-2-({3-[(1E)-2-(ピリジン-2-イル)エテン-1-イル]-1H-インダゾール-6-イル}スルファニル)ベンザミドである。アキシチニブについては、主としてCYP3A4/5によって代謝され、一部はCYP1A2、CYP2C19及びUGT1A1によって肝臓で代謝されることが知られている。アキシチニブは、効能及び効果を根治切除不能又は転移性の腎細胞癌とし、商品名：インライタとして販売されている。
【００２３】
ここで、アキシチニブの薬物動態とは、特に、アキシチニブの代謝及び排泄に関する体内動態を意味している。したがって、アキシチニブの薬物動態を示すパラメータとしては、特に限定されず、アキシチニブの代謝及び排泄に関するパラメータ、例えば全身クリアランス（CLtot）、血中濃度−時間曲線下面積（AUC）及び消失半減期を挙げることができる。
【００２４】
全身クリアランスは、時間あたりの薬物代謝・排泄量を体積に換算した値として定義される。すなわち、全身クリアランスは、単位を[容積]/[時間]とする値であり、例えば（mL/min）を単位として表される。また、全身クリアランスは、肝臓における代謝及び腎臓における排泄による薬物の消失量から求められる消失速度を、当該薬物の血中濃度で割ることで算出することができる。
【００２５】
血中濃度−時間曲線下面積とは、血中濃度を縦軸とし、時間を横軸としたグラフ（血中濃度−時間曲線）について、当該血中濃度を時間で積分した値である。つまり、血中濃度−時間曲線下面積とは、体内で利用された総薬物量を示す値である。また、上述した全身クリアランス（CLtot）は、薬物の投与量を血中濃度−時間曲線下面積（AUC）で割った値となる。
【００２６】
半減期とは、薬物が一次消失速度で減少する場合において所定の濃度から半分の濃度になるまでに要する時間を意味する。
【００２７】
[背景因子及び遺伝子多型]
アキシチニブの薬物動態を予測する際に使用する患者における所定の背景因子と当該患者における所定の遺伝子多型とについて詳述する。
【００２８】
アキシチニブの薬物動態を予測する際に使用する患者における所定の背景因子とは、当該患者に対するアキシチニブの投与量と、前治療有無に関する情報とを含む意味である。アキシチニブの投与量については、投与量そのものを説明変数としても良いし、投与量を所定の範囲に分けて当該範囲毎に所定の値を説明変数として割り当てても良い。投与量を所定の範囲に分けて説明変数を設定する際には、予測式の立て方にもよるが、投与量が多くなるにつれて高い数値を説明変数とすることができる。
【００２９】
また、前治療有無に関する情報については、予測式の立て方にもよるが、例えば、前治療無しの患者を０とし、前治療有りの患者を１とする説明変数を設定することができる。ここで、前治療とは、癌に対する化学療法及び免疫療法を含む意味である。化学療法としては、ソラフェニブ、スニチニブ、エベロリムス、テムシロリムス及びパゾパニブ等のアキシチニブを除く分子標的薬の投与が挙げられる、免疫療法としてはインターフェロン（INF）及びインターロイキン-2の投与が挙げられる。
【００３０】
アキシチニブの薬物動態を予測する際に使用する患者における遺伝子多型とは、以下の[1]〜[5]に規定される遺伝子多型である。
[1]  UGT1A1遺伝子の多型であって、変異型アレルが機能欠損型変異である多型
[2]  UGT1A7遺伝子の多型であって、変異型アレルが機能欠損型変異である多型
[3]  UGT1A9遺伝子の多型であって、変異型アレルが機能獲得型変異である多型
[4]  BCRP遺伝子の多型であって、変異型アレルが機能欠損型変異である多型
[5]  MDR1遺伝子の多型であって、変異型アレルが機能欠損型変異である多型
【００３１】
上記[1]の規定に含まれる多型としては、UGT1A1＊6及びUGT1A1＊28並びにこれらに対して連鎖不平衡となる多型を挙げることができる。
【００３２】
UGT1A1の遺伝子多型の１つであるUGT1A1＊6は、エキソン１に含まれる211番目が、多数を占める野生型（UGT1A1*1）ではグアニン（G）であるのに対し、アデニン（A）になっている遺伝子多型である。この１塩基の相違によりUGT1A1タンパク質における71番目のアミノ酸がグリシンからアルギニンに変異し、変異型のUGT1A1酵素は、野生型と比較して酵素活性が低下することとなる（機能欠損型変異）。
【００３３】
UGT1A1の遺伝子多型の１つであるUGT1A1＊28はプロモーター領域のTATAボックスにおけるTA配列の繰り返しが、多数を占める野生型（UGT1A1*1）では6回であるのに対し、7回の繰り返しになっている遺伝子多型であり、このTA2塩基挿入という違いにより遺伝子の発現量が低下し、結果としてUGT活性が低下する（機能欠損型変異）。
【００３４】
上記[2]の規定に含まれる多型としては、UGT1A7＊12及びこれに対して連鎖不平衡となる多型を挙げることができる。UGT1A7の遺伝子多型の１つであるUGT1A7＊12は、プロモーター領域の-57番目、コーディング領域の622番目及び760番目が、多数を占める野生型（UGT1A7*1）ではそれぞれチミン、チミン及びシトシンであるのに対し、それぞれグアニン、シトシン及びチミンになっている遺伝子多型である。プロモーター領域の１塩基の相違により遺伝子の発現量が低下し、また、コーディング領域の2塩基の相違によりUGT1A7タンパク質における208番目のアミノ酸がトリプトファンからアルギニンに変異し、且つ、254番のアミノ酸がアルギニンからストップコドンに変異することとなる。この変異型のUGT1A7酵素は、野生型と比較して発現量及び酵素活性がともに低下することとなる（機能欠損型変異）。なお、以下の説明において、UGT1A7＊12のうち、プロモーター領域の-57番目に位置する多型についてはUGT1A7＊12(-57)と称する。
【００３５】
上記[3]の規定に含まれる多型としては、UGT1A9＊1b及びこれに対して連鎖不平衡となる多型を挙げることができる。UGT1A9の遺伝子多型の１つであるUGT1A9＊1bはプロモーター領域の-118番目から始まる(dT)が、野生型が9回（(dT)9）であるのに対し、10回になっている遺伝子多型（(dT)10）であり、このチミンの１塩基挿入という違いにより遺伝子の発現量が上昇し、結果としてUGT活性が上昇する（機能獲得型変異）。
【００３６】
上記[4]の規定に含まれる多型としては、BCRP遺伝子におけるC421A多型及びこれに対して連鎖不平衡となる多型を挙げることができる。BCRPの遺伝子多型の１つであるC421A多型は、コーディング領域の421番目のシトシンがアデニンになっている遺伝子多型である。この一塩基の相違により、BCRPタンパク質における141番目のグルタミンがリシンに変異し、基質薬物を補足して排出するといった輸送機能が著しく低下する（機能欠損型変異）。
【００３７】
上記[5]の規定に含まれる多型としては、MDR1遺伝子におけるC1236T多型及びこれに対して連鎖不平衡となる多型を挙げることができる。MDR1の遺伝子多型の１つであるT1236C多型は、コーディング領域の1236番目のシトシンがチミンになっている遺伝子多型である。この一塩基の相違はMDR1タンパク質におけるアミノ酸置換を伴わないが、mRNAの安定性に影響を与え、マイナーアリルであるTを持つ遺伝子型ではmRNA量が低下する（機能欠損型変異）。
【００３８】
ところで、上記[1]〜[5]の規定における、連鎖不平衡となる多型とは、特に限定されないが、連鎖不平衡係数D’が0.8以上、より好ましくは0.95以上、さらに好ましくは0.99以上、最も好ましくは１である連鎖不平衡状態にある多型とすることができる。連鎖不平衡とは、２つの対立遺伝子がそれぞれ独立に遺伝する場合よりも大きな頻度で互いに連鎖して遺伝することをいう。連鎖不平衡係数D’は、２つの遺伝子多型について第一の遺伝子多型の各アレルを（A，a）、第二の遺伝子多型の各アレルを（B，b）とし、４つのハプロタイプ（AB，Ab，aB，ab）の各頻度をそれぞれPAB，PAb，PaB，Pabとすると、下記式により得られる。
【００３９】
D’＝（PABPab-PAbPaB）/Min[（PAB+PaB）（PaB+Pab），（PAB+PAb）（PAb+Pab）]
[式中、Min[（PAB+PaB）（PaB+Pab），（PAB+PAb）（PAb+Pab）]は、（PAB+PaB）（PaB+Pab）と（PAB+PAb）（PAb+Pab）とのうち、値の小さい方をとることを意味する。]
【００４０】
なお、上記[1]〜[5]に規定された遺伝子多型の有無については、予測式の立て方にもよるが、野生型アレルをホモで有する場合には０、当該遺伝子多型をヘテロで有する場合には１、そして、当該遺伝子多型をホモで有する場合には２といった説明変数を設定することができる。ただし、上記[1]〜[5]に規定された遺伝子多型について設定する説明変数は、この具体例に限定されず、上記の順に０、１及び３といった非等間隔に設定しても良い。
【００４１】
[予測式]
上述したように、患者における所定の背景因子と当該患者における所定の遺伝子多型とを説明変数として、アキシチニブの薬物動態パラメータを目的変数とする予測式によって、患者毎の予測薬物動態パラメータを計算する。予測式については、何ら限定されるものではなく、実際の患者における上記背景因子及び上記遺伝子多型並びに実測したCLtot値を用いた回帰分析によって得ることができる。回帰分析には、特に限定されないが、例えばべき乗モデル、指数モデル、漸近指数モデル、ロジスティク成長モデル及びゴンペルツ成長モデル等の非線形重回帰分析手法を使用することができる。
【００４２】
このようにして得られた予測式は、一例として実施例に示されるが、上述した患者における所定の背景因子と当該患者における所定の遺伝子多型から予測全身クリアランス（予測CLtot値）を算出することができる。この場合、予測式は、上記[1]〜[5]の遺伝子多型についてそれぞれ独立した項、アキシチニブの投与量に関する項及び前治療有無に関する項を含み、各項に係る係数が設定される。
【００４３】
特に、予測式は、上記[1]の遺伝子多型を有する場合に予測CLtot値が高くなり、上記[2]の多型を有する場合に予測CLtot値が低くなり、上記[3]の多型を有する場合に予測CLtot値が高くなり、上記[4]の多型を有する場合に予測CLtot値が低くなり、上記[5]の多型を有する場合に予測CLtot値が低くなる予測式である。また、この場合、予測式は、アキシチニブの投与量が多量であるほど予測CLtot値が高くなり、前治療有りの場合に予測CLtot値が高くなるものとする。
【００４４】
また、予測式で算出される予測CLtot値は、上記[1]〜[5]の遺伝子多型のうち、上記[1]の多型を有する場合に算出される予測CLtot値への影響が最も高く、次に上記[2]の多型を有する場合に算出される予測CLtot値への影響が高く、次に上記[7]の前治療が有る場合に算出される予測CLtot値への影響が高くなる式とすることが好ましい。ここで、算出される予測CLtot値への影響は、係数の絶対値の大きさで調整することができる。例えば、上記[1]の多型に関する項に係る係数の絶対値を、上記[2]の多型に関する項に係る係数の絶対値より大きくすることで、算出される予測CLtot値に対して上記[1]の多型を有する場合の影響が上記[2]の多型を有する場合の影響より相対的に大きくすることができる。
【００４５】
なお、薬物動態パラメータを予測する予測式としては、上記[1]〜[5]に規定する遺伝子多型に関する項、アキシチニブ投与量に関する項及び前治療有無に関する項に加えて、他の遺伝子多型に関する項や他の背景要因に関する項を含むように設定しても良い。他の遺伝子多型としては、例えばUGT1A1遺伝子、UGT1A7遺伝子又はUGT1A9遺伝子に含まれる他の遺伝子多型を挙げることができる。また、他の背景要因としては、例えば、患者年齢、性別、病理組織診断結果、既往症等を挙げることができる。
【００４６】
さらに、予測式で算出される予測CLtot値といった予測薬物動態パラメータには信頼区間を推定することができる。信頼区間としては、四分位範囲とすることもできるし、正規分布を仮定して95%信頼区間とすることもできる。予測式によって算出された予測薬物動態パラメータが信頼区間から外れる場合には、予測不能と判断することができる。
【００４７】
[アキシチニブによる治療効果判定]
上述のように、患者における所定の背景因子と当該患者における所定の遺伝子多型に基づいて当該患者についてアキシチニブの薬物動態パラメータを予測することができる。この予測された薬物動態パラメータに基づいて、当該患者に関するアキシチニブによる治療効果を判定することができる。
【００４８】
ここで、アキシチニブによる治療効果とは、治療対象の癌を改善する効果及びアキシチニブ投与による副作用の両方を含む意味である。すなわち、予測された薬物動態パラメータに基づいて、当該患者に対してアキシチニブを投与した場合の癌改善効果及び副作用の発生について高確度で判断することができる。
【００４９】
アキシチニブによる治療対象の癌としては、例えば、腎細胞癌、特に根治切除不能又は転移性の腎細胞癌を挙げることができる。ただし、本発明においてはこれに限定されず、アキシチニブによる治療対象の癌として線維肉腫、悪性線維性組織球腫 、脂肪肉腫、横紋筋肉腫、平滑筋肉腫、血管肉腫、カポジ肉腫、リンパ管肉腫、滑膜肉腫、軟骨肉腫、骨肉腫などの肉腫、脳腫瘍、頭頚部癌、乳癌、肺癌、食道癌、胃癌、十二指腸癌、虫垂癌、大腸癌、直腸癌、肝癌、膵癌、胆嚢癌、胆管癌、肛門癌、尿管癌、膀胱癌、前立腺癌、陰茎癌、精巣癌、子宮癌、卵巣癌、外陰癌、膣癌、皮膚癌などの癌腫、さらには白血病や悪性リンパ腫などが挙げられる。
【００５０】
したがって、患者について予測したアキシチニブの薬物動態パラメータに基づいて、これら癌に対するアキシチニブの効果を予測することができる。アキシチニブによる癌の改善効果としては、特に限定されないが、例えば腫瘍縮小率で評価することができる。腫瘍縮小率は、固形癌の治療効果判定のための新ガイドライン（RESICT version1.1）に基づき、「ベースライン評価時の全標的病変の径の和（非リンパ節病変では長径、リンパ節病変では短径）を、ベースライン径和(a) として算出し、治療後の全標的病変の径の和から引いた値(b)をベースライン径和で割った値(a)に100をかけたもの」、すなわち式：100 x [(a)-(b)]/(a)にて算出することができる。また、アキシチニブによる癌の改善効果としては、その他にもMRI画像やCT画像による腫瘍の大きさ数及び発生部位等に基づいて判断することもできる。
【００５１】
また、アキシチニブ投与による副作用としては、例えば、タンパク尿を挙げることができる。ただし、本発明においてはこれに限定されず、アキシチニブの副作用としては、下痢、高血圧、疲労、悪心、食欲減退、発声障害、手足症候群、甲状腺機能低下症、無力症、嘔吐、体重減少、粘膜の炎症、口内炎、発疹、便秘、頭痛、皮膚乾燥、味覚異常、TSH増加、悪心、AST（GOT）増加、ALT（GPT）増加、鼻出血、関節痛、ALP増加、腹痛、LDH増加、倦怠感、咳嗽、胸痛、血小板数減少及び浮腫を挙げることができる。また、アキシチニブ投与による重大な副作用としては、高血圧、高血圧クリーゼ、動脈血栓塞栓症、静脈血栓塞栓症、出血、消化管穿孔並びに瘻孔形成、甲状腺機能障害、創傷治癒遅延、可逆性後白質脳症症候群、肝機能障害及び心不全を挙げることができる。
【００５２】
したがって、患者について予測したアキシチニブの薬物動態パラメータに基づいて、これらアキシチニブによる副作用の発生の有無及び程度を予測することができる。
【００５３】
以上のように、予測された薬物動態パラメータに基づいて、患者毎にアキシチニブによる治療効果を判定することができる。この判定結果を用いることによって、患者毎にアキシチニブによる治療方針を最適化することができる。アキシチニブによる治療方針とは、アキシチニブの投与量及び投与タイミングを含む意味である。特に、この判定結果を用いることによって、患者毎にアキシチニブの投与量を最適化することが好ましい。
【００５４】
すなわち、予測された薬物動態パラメータにより、アキシチニブの代謝・排泄が遅いと判断される場合（例えば予測CLtot値が低値である場合）、副作用の発生が予測できるためアキシチニブ投与量を低くすることができる。逆に予測された薬物動態パラメータにより、アキシチニブの代謝・排泄が早いと判断される場合（例えば予測CLtot値が高値である場合）、副作用の可能性が低く且つ抗腫瘍活性も低いと予測できるためアキシチニブ投与量を高くすることができる。また、予測された薬物動態パラメータが中央値近傍（例えば25〜50パーセンタイルの範囲）であれば、通常のアキシチニブの投与量とすることができる。
【００５５】
なお、通常のアキシチニブの投与量とは、成人には1回5mgを1日2回経口投与（10mg/day）である。したがって、上述した例において、アキシチニブの投与量を低くするとは、10mg/day未満の量、例えば7mg/dayとすることになる。また、上述した例において、アキシチニブの投与量を高くするとは、10mg/dayより多い量、例えば14mg/dayとすることになる。
【００５６】
[薬物動態判定装置、治療効果判定装置]
本発明に係るアキシチニブの薬物動態判定装置は、上述したように、患者における所定の背景因子と当該患者における所定の遺伝子多型に基づいて当該患者についてアキシチニブの薬物動態パラメータを予測する装置である。
【００５７】
薬物動態判定装置は、例えば、図１に示すように、患者における所定の背景因子と当該患者における所定の遺伝子多型に関する情報を入力する入力部１と、入力部１で入力した情報に基づいてアキシチニブの予測薬物動態パラメータを算出する演算部２とを備えている。また、薬物動態判定装置は、入力部１で入力した情報や、演算部２で算出した予測薬物動態パラメータを記憶する記憶部３を備えていても良い。また、薬物動態判定装置は、演算部２で算出した予測薬物動態パラメータ、記憶部３に記憶されて予測薬物動態パラメータを出力する出力部４を備えていても良い。
【００５８】
なお、薬物動態判定装置は、いわゆるコンピュータにより実現することができる。すなわち薬物動態判定装置における入力部１は、コンピュータにおけるキーボードやマウス等の入力装置、或いは外部記憶装置（公的データベース等）との入力インターフェースを意味する。また、演算部２は、入力された情報に基づいて上述した予測式を用いて予測薬物動態パラメータを算出する演算機能を有するCPU等を意味する。さらに記憶部３は、コンピュータ内部のメモリやハードディスク、或いは外部記憶装置とすることができる。なお、記憶部３は、複数の薬物動態判定装置が共有することができる。
【００５９】
以上のように構成された薬物動態判定装置では、入力部１が患者における所定の背景因子と当該患者における所定の遺伝子多型を入力する。これら情報は操作者が入力しても良いし、例えば、遺伝子多型解析装置から入力しても良い。後者の場合、薬物動態判定装置は、遺伝子型解析装置と接続され、当該遺伝子型解析装置から出力された遺伝子型に関する情報が入力されることとなる。入力された背景因子及び/又は遺伝子多型に関する情報は、演算部２が記憶部３に記憶することもできる。
【００６０】
演算部２は、入力部１で入力された又は記憶部３に記憶された背景因子及び/又は遺伝子多型に関する情報と、予測薬物動態パラメータを算出するための予測式とを用いて予測薬物動態パラメータを算出する。ここで、予測式は、記憶部３に記憶されていても良いし、入力部１から入力しても良い。演算部２で算出された予測薬物動態パラメータは、例えばディスプレイやプリンタ等の出力部４に出力することができる。なお、演算部２で算出された予測薬物動態パラメータは、図１には図示しないが、インターネットやイントラネット等のネットワークを介して他のコンピュータやタブレット等の装置に出力することもできる。
【００６１】
また、演算部２は、算出した予測薬物動態パラメータが予め設定した上記信頼区間内に入るか外れるか判断し、予測薬物動態パラメータが信頼区間を外れる場合には予測不能と出力することができる。
【００６２】
一方、本発明に係るアキシチニブによる治療効果判定装置は、上述したように、患者における所定の背景因子と当該患者における所定の遺伝子多型に基づいて当該患者についてアキシチニブの薬物動態パラメータを予測し、予測された薬物動態パラメータに基づいてアキシチニブによる治療効果を判定する装置である。
【００６３】
治療効果判定装置は、例えば、図２に示すように、患者における所定の背景因子と当該患者における所定の遺伝子多型に関する情報を入力する入力部１と、入力部１で入力した情報に基づいてアキシチニブの予測薬物動態パラメータを算出し、当該予測薬物動態パラメータに基づいてアキシチニブによる治療効果を判定する演算部２とを備えている。また、治療効果判定装置は、入力部１で入力した情報や、演算部２で算出した予測薬物動態パラメータ及び/又は治療効果を記憶する記憶部３を備えていても良い。また、治療効果判定装置は、演算部２で算出した予測薬物動態パラメータ及び/又は治療効果、記憶部３に記憶されて予測薬物動態パラメータ及び/又は治療効果を出力する出力部４を備えていても良い。
【００６４】
なお、図２に示した治療効果判定装置において、図１に示した薬物動態判定装置と同一の構成は同一の符号を付することによってその詳細な説明を省略する。特に図２に示した治療効果判定装置では、演算部３は、予測された薬物動態パラメータに基づいてアキシチニブによる治療効果として、予測された腫瘍縮小率及び/又は副作用発生の可能性を判定する。また、治療効果判定装置は、これら治療効果の予測に基づいてアキシチニブを用いた治療方針（投与量及び/又は投与タイミング）を導き出しても良い。
【実施例】
【００６５】
以下、実施例により本発明をより詳細に説明するが、本発明の技術的範囲は以下の実施例に限定されるものではない。
【００６６】
〔実施例１〕
本実施例では、アキシチニブの治療効果や副作用との相関を認めた全身クリアランス（CLtot；薬物動態解析パラメータの一つ）を予測するモデルを、関連する遺伝子多型情報と患者背景因子を変数とした回帰モデルを用いて作成した。
【００６７】
実験方法
（１）アキシチニブ投与
ａ）対象患者
対象患者は、山口大学医学部付属病院  泌尿器科において20歳以上で腎細胞癌と診断されアキシチニブ治療を導入予定で「インフォームド・コンセントのための説明文書」の内容を理解し、「同意文書」に自署することにより自由意思により本研究に参加同意する進行性腎細胞癌患者である。
ｂ）アキシチニブ投与法
1) アキシチニブ標準投与開始量10mg（1回投与量5mg）経口投与をおこなった。
2) 臨床的に増量のパラメータとされる高血圧が出現しなかった場合、投与量を14mg（1回投与量7mg）に増量した。
3) Grade 2の有害事象（NCI CTCAE ver.4.0の有害事象の評価規準による）が出現した場合、休薬後投与量10mgで再開し、Grade3の有害事象が出現した場合、休薬後、投与量を7mg（一回投与量：朝5mg, 夕2mg）に減量再開した。
【００６８】
（２）遺伝多型解析
患者末梢血（10ml）単核球成分よりDNAを抽出（約1μL）し、2つの容器に分注した。一方のDNAをジーンチップ（東洋鋼鈑社より供与）基盤チップ上の鋳型DNAにハイブリダイズさせ、専用リーダーで蛍光強度を定量することにより、UGT1A1*6、UGT1A1*28、UGT1A1*60、UGT1A7*12、UGT1A7*2、UGT1A9*1b、UGT1A1*93、UGT1A4*1b、UGT1A4*3、UGT1A5*3、UGT1A5*6、BCRP（421C>A）、MDR1（1236T>C）、MDR1（2677G>T/A）及びMDR1（3435C>T）の遺伝子型を決定した。図３に遺伝子型を解析した結果を示した。
【００６９】
（３）薬物動態濃度解析
同意が得られたアキシチニブ投与予定患者に対し、薬物投与量が変更される毎に薬物動態濃度解析を行った。採血は、薬剤投与開始8日以降（投与中断があった場合は再開日より8日以降、すなわち連続する7日間の投与後）に行った。静脈用留置カテーテルに採血ルートを確保後、アキシチニブ投与前、投与後30分、1、2、3、4、8、12時間に採血（5ml）を行い、アキシチニブ血中濃度測定を行った。採血ルートは採血後ヘパリン注入（0.5 ml）にて複数回使用可能な状態とし、採血前には2mlの血液を廃棄後、サンプル用の採血を行った。測定は液体クロマトグラフィ法（LC MS/MS）で測定し、拡張最少二乗法により最高血清中濃度（Cmax）、最高血清中濃度到達時間（Tmax）、血清中濃度曲線下面積（AUC）、みかけの血清中消失半減期（T1/2）を算出し、extensive metabolizer及びslow metabolizerを決定した。
【００７０】
（４）癌縮小率
癌縮小率の計算方法
縮小率は、固形癌の治療効果判定のための新ガイドライン（RESICT version1.1）に基づき、「ベースライン評価時の全標的病変の径の和（非リンパ節病変では長径、リンパ節病変では短径）を、ベースライン径和(a) として算出し、治療後の全標的病変の径の和から引いた値(b)をベースライン径和で割った値(a)に100をかけたもの」とした。つまり100 x [(a)-(b)]/(a)で計算した。
【００７１】
（５）回帰モデルの構築
上記（３）及び（４）の結果を纏めて図４に示した。図３及び図４に示した結果を用いて、UGT1Aに関する6因子、CYP3Aに関する1因子、ABC関連遺伝子に関する4因子、年齢、性別、前治療の有無、投与量を加えた15因子でCLtot値を推定する回帰モデル（exponential regression model）を構築した。回帰モデルには、下記式で表される指数関数モデルを用いた。
【００７２】
【数１】
2018014938.tif 000003
【００７３】
ここで、遺伝子型、前治療の有無、投与量等の上記15因子をxiとし、CLtot値をyとして、図３及び４に示した40例に対して式のyと予測値（右辺）の二乗誤差が最小となるようにb0〜bnを求めた。なお、本実施例では、遺伝子型について野生型、ヘテロ型及び変異型をそれぞれ0、1及び2とおいた。
【００７４】
解析結果
全15因子を候補として，変数減少法（Backward Selection）により決定係数が最大となり因子数が最小となる組合せを探索した。図５に示すように、全因子を用いたときの決定係数R2＝0.48であったのに対して、UGT1A1遺伝子におけるUGT1A1*6及びUGT1A1*6*28、UGT1A7遺伝子におけるUGT1A7*12、UGT1A9遺伝子におけるUGT1A9*1b、BCRP遺伝子における421C>Aの多型、MDR1遺伝子における1236T>Cの多型、アキシチニブ投与量及び前治療の有無の8因子を用いたときの決定係数はR2=0.63であった。
【００７５】
また、図５に示すように、これら8因子に対してUGT1A7*2、年齢、又は性別を含めて決定係数を算出しても、これら8因子を用いた場合の決定係数と同等であることが判った。さらに、これら8因子のなかから１因子（MDR1遺伝子における1236T>Cの多型）減らした7因子で決定係数を算出すると、8因子を用いた場合の決定係数から0.04低下し、さらに2因子（UGT1A9*1bとMDR1遺伝子における1236T>Cの多型）減らした6因子で決定係数を算出すると、8因子を用いた場合の決定係数から0.12も低下することが明らかとなった。
【００７６】
なお、図５において、各因子の下に記載した数値は、上記式における係数を意味し、当該係数で上記式を規定することで予測全身クリアランス（予測CLtot値）を算出するための予測式となる。したがって、各因子に係る係数の絶対値は，予測CLtot値の予測における影響の強さを表している。係数が正であれば因子の数値が高いほど予測CLtot値が高い傾向にあり、係数が負であれば因子の数値が低いほど予測CLtot値が低い傾向にある。係数が0に近い（絶対値が小さい）ほど、予測における影響が小さいことを表している。
【００７７】
40例の患者について、それぞれ8因子を用いて予測CLtot値を算出した場合の、実測したCLtot値と比較した結果を図６に示した。図６からも、これら8因子を用いて算出した予測CLtot値は、実測したCLtot値と高い相関関係にあることが理解できる。
【００７８】
〔実施例２〕
本実施例では、実施例１で作製した予測CLtot値を算出するモデルを実際の臨床例へ適用した例を説明する。
【００７９】
先ず、実測したCLtot値と最大腫瘍縮小率との関係を図７に示した。図７に示したように、CLtot（実測値）と最大縮小率との間には相関があることが明らかとなった。そこで、予測CLtot値と最大縮小率との相関を調べることとした。実施例１で構築した予測モデルにより算出した予測CLtot値と最大腫瘍縮小率との関係を図８に示した。図８に示すように、実施例１で構築した予測モデルにより算出した予測CLtot値のうち、中央値以下の低値群は、腫瘍縮小（中央値：13.5%以上）と有意の相関を認め、その予測確率は感度、特異度、陽性的中率とも76.5%であった。
【００８０】
次に、予測CLtot値及び実測Cltot値とタンパク尿のグレードとの相関を調べた。結果を図９に示した。図９に示すように、予測CLtot値は副作用（タンパク尿）と有意の相関が認められた。
【００８１】
本実施例の結果より、実施例１で構築した予測モデルにより算出した予測CLtot値は、腎細胞癌に対するアキシチニブ治療効果（腫瘍改善効果及び副作用）の予測に有用であることが明らかとなった。
【図１】
2018014938.tif 000004
【図２】
2018014938.tif 000005
【図３】
2018014938.tif 000006
【図４】
2018014938.tif 000007
【図５】
2018014938.tif 000008
【図６】
2018014938.tif 000009
【図７】
2018014938.tif 000010
【図８】
2018014938.tif 000011
【図９】
2018014938.tif 000012