| Literature DB >> 32799335 |
Jing Li1, Jun Jiang1, Jianmei Wu1, Xun Bao1, Nader Sanai2.
Abstract
A better understanding of the human central nervous system (CNS) pharmacokinetics is critical to the selection of the right drug and refinement of dosing regimen for more effective treatment of primary and metastatic brain cancer. Using the physiologically-based pharmacokinetic (PBPK) modeling approach, we systematically compared the CNS pharmacokinetics of three cyclin D-cyclin dependent kinase 4 and 6 (CDK4/6) inhibitors (ribociclib, palbociclib, and abemaciclib) in patients with cancer. A PBPK model platform was developed and verified for predicting plasma and CNS pharmacokinetics. Target engagement ratio (TER), defined as the ratio of the average steady-state unbound drug brain concentration to the in vitro half-maximal inhibitory concentration (IC50 ) for CDK4/6 inhibition, was used as a crude predictor of efficacy. As compared with ribociclib and palbociclib, abemaciclib penetrated into the human brain to a larger extent, but at a slower rate, and was retained in the brain longer. Following the standard dosing regimens, the predicted CDK4/6 TERs were 26/5.2 for abemaciclib, 2.4/0.62 for ribociclib, and 0.36/0.27 for palbociclib. Simulations suggested that abemaciclib achieved comparable TERs following twice daily or daily dosing; ribociclib may sufficiently inhibit both CDK4 and CDK6 at the maximum tolerated dose; whereas, palbociclib achieved TERs < 0.5 even at a dose 50% higher than the standard dose. In conclusion, the PBPK modeling, supported by available preclinical and clinical evidence, suggests that abemaciclib is the best CDK4/6 inhibitor for brain cancer treatment, whereas palbociclib is not recommended. The model refined dosing regimen is 300 mg daily on a 4-weeks-on schedule for abemaciclib, and 900 mg daily on a 3-weeks-on/1-week-off schedule for ribociclib.Entities:
Year: 2020 PMID: 32799335 PMCID: PMC7854954 DOI: 10.1002/cpt.2021
Source DB: PubMed Journal: Clin Pharmacol Ther ISSN: 0009-9236 Impact factor: 6.875
Figure 1Model structure of the whole‐body physiologically‐based pharmacokinetic model integrated with a 4‐compartment permeability‐limited brain (4Brain) model, which has 4 compartments representing the brain blood, brain mass, cranial cerebrospinal fluid (CSF), and spinal CSF. The 4Brain model assumes: (1) drug transport across the blood‐brain barrier (BBB) is governed by bidirectional passive permeability (PSB) and ABCB1‐mediated and ABCG2‐mediated active efflux clearance (CLefflux,BBB); (2) drug transport across the blood‐CSF barrier was controlled by bidirectional PSB and ABCB1‐mediated active influx clearance (CLuptake,CSF); (3) drug transport between the brain mass and cranial CSF is diffusive transport at a bi‐directional clearance (PSB); (4) fluid balance is maintained by the circulation of CSF between spinal and cranial compartments and reabsorbed into the brain blood; (5) the cerebral blood flow rate (Q Brain) links the 4Brain model to whole‐body model; (6) only unbound and unionized drug can passively pass through all barriers, while transporters act upon unbound drug (including both unionized and ionized species); (7) metabolism in brain mass is negligible; and (8) all compartments are well‐stirred with defined volumes. Flow rates are described by the CSF secretion rate (Q BCSFB), bulk flow rate from brain mass to cranial CSF (Q bulk), CSF flow rate out of cranial and spinal compartments (Q sink), CSF shuttle flow rate between cranial and spinal compartments (Q Sin and Q Sout), and water transfer rate from the brain blood to brain mass (Q BBB).
Drug‐specific parameters in the whole‐body‐4Brain PBPK models
| Ribociclib | Palbociclib | Abemaciclib | Comments/reference | |
|---|---|---|---|---|
| Physicochemical | ||||
| MW, g/mol | 434.6 | 447.5 | 506.6 | PubChem |
| LogP | 2.38 | 2.77 | 4.25 | ChemAxon |
| PKa1, strongest acidic | 11.59 | 11.34 | 10.27 | ChemAxon |
| PKa2, strongest basic | 8.87 | 8.86 | 7.94 | ChemAxon |
| B/P | 1.5 | 1.9 | 1.2 | Investigator brochure |
|
| 0.125 | 0.126 | 0.06 | Experimental determined |
| Absorption | ||||
| Fa | 0.9 | 0.5 | 0.6 | Assigned |
| Ka | 0.8 | 0.367 | 0.197 | Observed |
| Lag time, hours | 0.5 | 0.658 | 0.5 | Observed |
| fugut | 0.0061 | 0.007 | 0.004 | Simcyp predicted |
|
| 9.61 | 14.48 | 11.08 | Simcyp predicted |
| Distribution | ||||
|
| 12.19 | 16.31 | 7.86 | Predicted by Method 3, |
| Elimination | ||||
| HLMs: | ||||
| Vmax, pmol/min/mg | 313.7 | 518.5 | 2003 | Experimental determined from HLM |
|
| 13.4 | 78.1 | 126.2 | Experimental determined from HLM |
| CLint, µL/min/mg | 23.4 | 6.64 | 15.9 | Experimental determined from HLM |
| fuinc | 0.35 | 0.2 | 0.2 | Experimental determined from HLM |
| HIMs: | ||||
| Vmax, pmol/min/mg | 104.8 | 280.7 | 494.3 | Experimental determined from HIM |
|
| 28.8 | 70.4 | 115.1 | Experimental determined from HIM |
| CLint, µL/min/mg | 3.64 | 3.99 | 4.29 | Experimental determined from HIM |
| fumic | 0.35 | 0.2 | 0.2 | Experimental determined |
| CLR, L/h | 7.5 | 10.5 | 1.5 | Assigned |
| 4Brain model | ||||
| BBB | ||||
| PSB, L/h | 135 | 307 | 25 | Estimated by Eq. 1 |
| fu,br | 0.044 | 0.015 | 0.006 | Experimental determined from patient brain tumor tissue |
| CLABCB1,vitro, µL/min/mg | 21770 | 26160 | 800 | Experimental determined Eq. 2 |
| CLABCG2,vitro, µL/min/mg | 0 | 8580 | 90 | Experimental determined Eq. 2 |
| ABCB1 RAF at BBB | 87.3 | 87.3 | 87.3 | Determined based on ABCB1 abundance in normal human brain microvessels by Eq. 3 |
| ABCB1 RAF at BBB | 5.88 | 5.88 | 5.88 | Determined based on ABCB1 abundance in human glioblastoma microvessels by Eq. 3 |
| ABCG2 RAF at BBB | 125.3 | 125.3 | 125.3 | Determined based on ABCG2 abundance in normal human brain microvessels by Eq. 3 |
| ABCG2 RAF at BBB | 55.47 | 55.47 | 55.47 | Determined based on ABCG2 abundance in human glioblastoma microvessels by Eq. 3 |
| Blood‐CSF barrier | ||||
| PSC, L/h | 67.5 | 150 | 13 | Assumed to be half of PSB |
| fucsf | 1 | 1 | 1 | Assigned given low CSF protein concentration |
| CLABCB1,vitro, µL/min/mg | 21770 | 26160 | 800 | Experimental determined Eq. 2 |
| ABCB1 RAF at BCCSF | 8.7 | 8.7 | 8.7 | Assigned and validated by observed ribociclib CSF data |
| Brain‐cranial CSF barrier | ||||
| PSE, L/h | 300 | 300 | 300 | Assigned assuming no barrier function |
BBB, blood‐brain barrier; B/P, blood‐to‐plasma partition ratio; CLABCB1,vitro and CLABCB1,vitro, ABCB1‐mediated and ABCG2‐mediated in vitro efflux clearance, respectively; CLint, in vitro intrinsic metabolic clearance; CLR, renal clearance; f u,br, fraction unbound drug in brain tissue; f u,p, fraction of unbound drug in plasma; fucsf, fraction unbound in cerebrospinal fluid; fugut, fraction of unbound drug in enterocytes; fumic, fraction of unbound in microsomal incubation; K m, substrate concentration at which half of Vmax is achieved; logP, logarithm of the neutral species octanol‐to‐buffer partition ratio; MW, molecule weight; PKa, acid dissociation constant; PBPK, physiologically‐based pharmacokinetic; PSB, passive permeability‐surface area product at the BBB; PSC, passive permeability‐surface area product at the blood‐cranial CSF barrier; PSE, passive permeability‐surface area product at the brain‐cranial CSF barrier; Q gut, gut blood flow; RAF, in vivo‐in vitro relative activity factor; Vmax, maximum metabolic rate; V ss, volume of distribution at steady‐state.
(Eq. 1), where P app,A→B is the apparent permeability determined from MDCKII cell monolayer; SA is the human brain microvasculature surface area (mean, 20 m2); and λ is unionization efficiency.
(Eq. 2), where CLefflux,vitro (µL/min/mg) is the in vitro efflux transporter‐mediated intrinsic clearance; NER is the net efflux ratio determined from MDCKII‐ABCB1or MDCKII–ABCG2; P app,A‐B is the apparent passive permeability determined from MDCKII; SA is the filter surface area (0.143 cm2) in a 96‐well transwell; λ is unionization efficiency; and Procell is the cell membrane protein amount in a 96‐well transwell.
(Eq. 3), where RAF is the relative activity factor; BMvPGB is the milligrams of brain microvessels per gram brain; BW is the average human brain weight; abundance in vivo or in vitro represents the ABCB1/ABCG2 transporter protein expression level in human brain microvessels or in MDCKII‐ABCB1 and –ABCG2 cells, respectively.
Transcellular permeability and efflux ratios of the three CDK4/6 inhibitors, determined from the parental MDCKII, MDCKII‐ABCB1, and MDCKII‐ABCG2 cell monolayers
| Ribociclib | Palbociclib | Abemaciclib | |
|---|---|---|---|
| Parental MDCKII | |||
| Apparent permeability ( | 7.17 ± 1.85 | 9.38 ± 1.16 | 9.09 ± 2.21 |
| Net efflux ratio | 1.02 ± 0.21 | 1.55 ± 0.17 | 1.03 ± 0.30 |
| Unionization efficiency ( | 0.033 | 0.034 | 0.224 |
| Intrinsic permeability ( | 218.6 | 280.0 | 40.6 |
| MDCKII‐ABCB1 | |||
| Apparent permeability ( | 1.77 ± 0.82 | 5.56 ± 1.10 | 5.65 ± 0.86 |
| Net efflux ratio | 10.5 ± 3.30 | 5.92 ± 1.04 | 2.84 ± 0.49 |
| MDCKII‐ABCG2 | |||
| Apparent permeability | 5.90 ± 0.91 | 7.73 ± 0.25 | 7.98 ± 0.42 |
| Net efflux ratio | 1.02 ± 0.20 | 2.64 ± 1.05 | 1.21 ± 0.11 |
Data are expressed as the mean ± SD from at least three independent experiments (with duplicates in each experiment).
Bidirectional permeability experiments were performed at pH 7.4 in both apical and basolateral chambers.
Net efflux ratio was the efflux ratio in the absence of an ABCB1 inhibitor (elacridar) or ABCG2 inhibitor (Ko145) divided by the efflux ratio in the presence of the inhibitor.
Unionization efficiency (λ) is the ratio of unionized form to total drug (the sum of unionized and ionized forms), where the unionized‐to‐ionized ratio is calculated based on Henderson‐Hasselbalch equation: .
Intrinsic permeability is the transcellular permeability of the unbound and unionized drug, estimated as the mean P app,A‐B normalized by λ.
Figure 2The pH‐dependent apparent permeability of the three cyclin D‐cyclin dependent kinase 4 and 6 inhibitors across the MDCKII, MDCKII‐ABCB1, and MDCKII‐ABCG2 cell monolayers. Bi‐directional permeability experiments were performed at a fixed apical pH (7.4) and varying basolateral pH (7.4, 7.0, 6.5, and 6.0). The apparent permeability in the apical‐to‐basolateral (P app,A‐B) and basolateral‐to‐apical (P app,B‐A) directions and efflux ratio (ER) are expressed as the mean fold change relative to those determined in the control (at pH 7.4 in both apical and basolateral chambers).
Figure 3Physiologically‐based pharmacokinetic model‐simulated and clinically observed plasma and central nervous system (CNS) pharmacokinetics of ribociclib. (a) Simulated and observed unbound ribociclib plasma concentration‐time profiles. (b) Simulated and observed ribociclib cerebrospinal fluid concentration‐time profiles. (c) Simulated unbound ribociclib brain concentration‐time profiles in the human normal brain (with brain pH 7.12 and ABCB1 abundance 3.38 pmol/mg at the BBB). (d) Simulated unbound ribociclib concentration‐time profiles in brain tumors with pH 6.8 and ABCB1 abundance 0.14 pmol/mg at the BBB. (e) Simulated unbound ribociclib concentration‐time profiles in brain tumors with pH 6.8 and loss of ABCB1 expression at the BBB. (f) Simulated unbound ribociclib concentration‐time profiles in brain tumors with pH 6.5 and loss of ABCB1 expression at the BBB. Simulations of 10 trials with 10 subjects in each were performed in the Simcyp virtual cancer patient population following 5‐day ribociclib treatment at a daily oral dose of 900 mg. Observed clinical plasma and CNS pharmacokinetic data were obtained from 12 glioblastoma patients treated with ribociclib (900 mg q.d. for 5 days).
PBPK model‐predicted unbound drug exposure in the plasma, normal brain, cranial CSF, and spinal CSF at the steady‐state following the standard or modified dosing regimens
| Ribociclib | Palbociclib | Abemaciclib | ||||
|---|---|---|---|---|---|---|
| Standard (600 mg q.d., 21 days) | Alternative (900 mg q.d., 21 days) | Standard (125 mg q.d., 21 days) | Alternative (187 mg q.d., 21 days) | Standard (150 mg b.i.d., 28 days) | Alternative (300 mg q.d., 28 days) | |
| Plasma | ||||||
| Tss,max, hours | 1.32 | 1.32 | 3.84 | 3.84 | 2.91 | 3.84 |
| Css,max, µmol/L | 0.298 | 0.448 | 0.038 | 0.057 | 0.055 | 0.062 |
| Css,min, µmol/L | 0.092 | 0.139 | 0.022 | 0.033 | 0.045 | 0.040 |
| Css,ave, µmol/L | 0.176 | 0.264 | 0.029 | 0.044 | 0.050 | 0.050 |
| AUC24h, µmol/L*h | 4.291 | 6.440 | 0.779 | 1.167 | 1.683 | 1.326 |
| Brain | ||||||
| Tss,max, hours | 3.84 | 3.84 | 3.84 | 3.84 | 6.27 | 10.56 |
| Css,max, µmol/L | 0.039 | 0.059 | 0.005 | 0.007 | 0.052 | 0.053 |
| Css,min, µmol/L | 0.014 | 0.020 | 0.003 | 0.004 | 0.051 | 0.050 |
| Css,ave, µmol/L | 0.024 | 0.036 | 0.004 | 0.005 | 0.052 | 0.052 |
| AUC24h, µmol/L*h | 0.624 | 0.935 | 0.096 | 0.143 | 1.746 | 1.391 |
| Brain | 0.14 | 0.14 | 0.12 | 0.12 | 1.04 | 1.05 |
| TER for CDK4 | 2.4 | 3.6 | 0.36 | 0.45 | 26 | 26 |
| TER for CDK6 | 0.62 | 0.92 | 0.27 | 0.33 | 5.2 | 5.2 |
| Cranial CSF | ||||||
| Tss,max, hours | 1.32 | 1.32 | 3.84 | 3.84 | 6.27 | 10.56 |
| Css,max, µmol/L | 0.475 | 0.713 | 0.068 | 0.102 | 0.039 | 0.040 |
| Css,min, µmol/L | 0.149 | 0.223 | 0.038 | 0.058 | 0.038 | 0.036 |
| Css,ave, µmol/L | 0.281 | 0.421 | 0.052 | 0.078 | 0.038 | 0.038 |
| AUC24h, µmol/L*h | 6.890 | 10.33 | 1.385 | 2.078 | 1.311 | 1.045 |
| CSF | 1.60 | 1.60 | 1.78 | 1.78 | 0.78 | 0.79 |
| Spinal CSF | ||||||
| Tss,max, hours | 3.84 | 3.84 | 6.36 | 6.36 | 6.27 | 10.56 |
| Css,max, µmol/L | 0.404 | 0.606 | 0.063 | 0.094 | 0.039 | 0.039 |
| Css,min, µmol/L | 0.159 | 0.238 | 0.040 | 0.060 | 0.038 | 0.037 |
| Css,ave, µmol/L | 0.263 | 0.394 | 0.051 | 0.076 | 0.038 | 0.038 |
| AUC24h, µmol/L*h | 6.742 | 10.11 | 1.379 | 2.069 | 1.293 | 1.030 |
| CSF | 1.57 | 1.57 | 1.77 | 1.77 | 0.77 | 0.78 |
AUC24h, area under the concentration time curve during 24 hours at the steady‐state; CDK, cyclin D‐cyclin dependent kinase; CSF, cerebrospinal fluid; Css,ave, average steady‐state concentration; Css,max, maximum steady‐state concentration; Css,min, trough steady‐state concentration; PBPK, physiologically‐based pharmacokinetic; TER, target engagement ratio.
Simulations of 10 trials with 10 subjects in each trial were performed in the Simcyp cancer patient population. Data are presented as the population mean values.
Brain K p,uu is estimated as the AUC24h ratio of unbound brain to unbound plasma at the steady‐state.
TER is calculated as the ratio of the average steady‐state unbound brain concentrations to the in vitro IC50 for inhibiting CDK4/6 enzymes.
CSF K p,uu is estimated as the AUC24h ratio of CSF to unbound plasma at the steady‐state.
Figure 4Physiologically‐based pharmacokinetic model‐simulated plasma and central nervous system pharmacokinetics of the three cyclin D‐cyclin dependent kinase 4 and 6 inhibitors following varying dosing regimens. (a) Simulated population mean plasma, CSF, and brain concentration time profiles of unbound ribociclib following a single oral dose (600 mg). (b) Simulated population mean plasma, CSF, and brain concentration time profiles of unbound palbociclib following a single oral dose (125 mg). (c) Simulated population mean plasma, CSF, and brain concentration time profiles of unbound abemaciclib following a single oral dose (150 mg). (d, e) Simulated unbound ribociclib brain concentration time profiles following the standard dosing regimen (600 mg q.d. for 3‐weeks) and modified dosing regimen (900 mg q.d. for 3 weeks). (f, g) Simulated unbound palbociclib brain concentration time profiles following the standard dosing regimen (125 mg q.d. for 3‐weeks) and modified dosing regimen (187 mg q.d. for 3 weeks). (h, i) Simulated unbound abemaciclib brain concentration time profiles following the standard dosing regimen (150 mg b.i.d. for 4‐weeks) and modified dosing regimen (300 mg q.d. for 4 weeks). Simulations of 10 trials with 10 subjects in each trial were performed in the Simcyp virtual cancer patient population.