Optimising the combination dosing strategy of abemaciclib and vemurafenib in BRAF-mutated melanoma xenograft tumours.

BACKGROUND: Resistance to BRAF inhibition is a major cause of treatment failure for BRAF-mutated metastatic melanoma patients. Abemaciclib, a cyclin-dependent kinase 4 and 6 inhibitor, overcomes this resistance in xenograft tumours and offers a promising drug combination. The present work aims to characterise the quantitative pharmacology of the abemaciclib/vemurafenib combination using a semimechanistic pharmacokinetic/pharmacodynamic modelling approach and to identify an optimum dosing regimen for potential clinical evaluation.
METHODS: A PK/biomarker model was developed to connect abemaciclib/vemurafenib concentrations to changes in MAPK and cell cycle pathway biomarkers in A375 BRAF-mutated melanoma xenografts. Resultant tumour growth inhibition was described by relating (i) MAPK pathway inhibition to apoptosis, (ii) mitotic cell density to tumour growth and, under resistant conditions, (iii) retinoblastoma protein inhibition to cell survival.
RESULTS: The model successfully described vemurafenib/abemaciclib-mediated changes in MAPK pathway and cell cycle biomarkers. Initial tumour shrinkage by vemurafenib, acquisition of resistance and subsequent abemaciclib-mediated efficacy were successfully captured and externally validated. Model simulations illustrate the benefit of intermittent vemurafenib therapy over continuous treatment, and indicate that continuous abemaciclib in combination with intermittent vemurafenib offers the potential for considerable tumour regression.
CONCLUSIONS: The quantitative pharmacology of the abemaciclib/vemurafenib combination was successfully characterised and an optimised, clinically-relevant dosing strategy was identified.

The recent approval of the BRAF inhibitors, vemurafenib and dabrafenib, represents a major improvement in the prognosis of patients with BRAF-mutated (V600E/K) metastatic melanoma. In a pivotal phase III registration trial, vemurafenib demonstrated significantly improved progression-free survival (PFS) over dacarbazine (5.3 vs 1.6 months) (Chapman ), with a similar improvement also observed for dabrafenib (Hauschild ). However, the majority of patients who initially responded to either vemurafenib or dabrafenib developed resistance and ultimately relapsed after 6–12 months of treatment (Chapman ; Hauschild ). Although combination with an MEK inhibitor further improved PFS (9.4 vs 5.8 months), it was unable to prevent acquisition of resistance to BRAF inhibition (Flaherty ). As such, resistance to vemurafenib/dabrafenib remains an area of significant unmet medical need for patients with BRAF-mutated metastatic melanoma.

Tumour biopsies obtained from relapsed patients typically exhibit MAPK reactivation (Trunzer ; Shi ), as signified by elevation of ERK1/2 phosphorylation (pERK) (Trunzer ) and Cyclin D1 (Smalley ). This mode of resistance to BRAF inhibition has been successfully replicated in preclinical xenograft studies in mice with human BRAF-mutant melanoma cell lines (Das Thakur ; Yadav ). Using these resistant cell lines, it was shown that although vemurafenib still inhibits pERK, the MAPK pathway is upregulated such that maximum vemurafenib-mediated inhibition cannot sufficiently suppress pERK to induce apoptosis (Das Thakur ). Furthermore, while a period of chronic vemurafenib therapy induces pERK elevation, temporary cessation of treatment allows pERK to return to near-normal levels and subsequent reintroduction of vemurafenib becomes efficacious once more (Das Thakur ). In the clinic, patients receiving vemurafenib on an intermittent schedule have demonstrated response to treatment (Dooley ; Koop ), with some evidence of reduced toxicity (Dooley ). Furthermore, those rechallenged with vemurafenib after previously progressing on treatment exhibited further response to BRAF inhibition (Seghers ; Romano ). However, all clinical evidence of the benefit of an intermittent schedule for BRAF inhibitors to date is restricted to isolated case studies.

While intermittent therapy may delay resistance to vemurafenib, using drug combinations to prevent or exploit MAPK pathway reactivation may further improve clinical efficacy over BRAF inhibition alone. The potential combination with a cyclin-dependent kinase 4 and 6 (CDK4/6) inhibitor, such as abemaciclib or palbociclib, is of particular interest (Sheppard and McArthur, 2013; Yadav ). Abemaciclib recently demonstrated a level of efficacy in vemurafenib-resistant A375 xenograft tumours which was far greater than that expected in treatment-naïve tumours (Yadav ). Using western blotting, Yadav noted increased cleaved PARP levels in vemurafenib-resistant cells treated with retinoblastoma protein (Rb) siRINA, suggesting increased dependency on Rb expression for cell survival. Chronic dosing with a CDK4/6 inhibitor such as abemaciclib suppresses phosphorylation of Rb (pRb) (Gelbert ; Tate ), which in turn downregulates Total Rb (Yadav ) as a result of the inherent autoregulation of the E2F1-Rb pathway (Shan ). The authors therefore proposed that abemaciclib-mediated inhibition of Total Rb in vemurafenib-resistant tumours may provide additional clinical benefit over the currently approved therapies for patients with BRAF-mutated metastatic melanoma (Yadav , 2014).

While intermittent vemurafenib and combined vemurafenib/abemaciclib therapy have proved promising in preclinical studies (Das Thakur ; Yadav ), the optimal dosing schedule that could maximise its clinical efficacy remains to be identified. The objective of this study was therefore to explore dosing strategies for the combination of abemaciclib and vemurafenib by means of a semimechanistic PK/PD model, which captures the quantitative pharmacology of abemaciclib overcoming resistance to vemurafenib. To this end, an existing semimechanistic PK/PD model for abemaciclib (Tate ) was further developed to include vemurafenib-mediated inhibition of the MAPK pathway and tumour shrinkage, acquisition of resistance through MAPK pathway reactivation and subsequent abemaciclib-mediated efficacy as a result of Total Rb inhibition. This model, once externally validated, was used to identify the combination treatment schedule yielding the best predicted response in A375 xenograft tumours. Model simulations not only supported existing preclinical studies reporting the benefit of intermittent vemurafenib therapy over continuous treatment (Das Thakur ) but also indicated that daily administration of abemaciclib in combination with intermittent vemurafenib therapy offers the potential for considerable tumour regression. To our knowledge, this is the first PK/PD modelling approach to formalise acquisition of resistance to vemurafenib in BRAF-mutated melanoma and to mechanistically evaluate dose schedule optimisation for combination therapies in oncology.

Materials and methods

In vivo experiments

Biomarker modulation in colo-205 xenograft tumours

The in vivo studies evaluating abemaciclib-mediated cell cycle arrest in colo-205 xenograft tumours used in this work were described previously (Tate ). In addition to phospho-(ser780)-Rb (pRb), topoisomerase II α (TopoIIα) and phospho-(ser10)-histone H3 (pHH3), tumour lysates from selected studies were also used to evaluate inhibition of Total Rb after abemaciclib. In the first study, mice received a single oral 50 mg kg−1 dose and tumour tissue was collected at 1, 6, 24, 36 and 48 h; in the second study, mice received 50 mg kg−1 orally once a day for 56 days, with tumour tissue collections at 28 and 56 days (n=5 per time point). A summary of the available experimental data is provided in Supplementary Table 1.

Biomarker modulation in A375 xenograft tumours

A series of in vivo experiments were performed to evaluate the dose–response and time course of abemaciclib- or vemurafenib-mediated modulation of biomarkers in mice bearing A375 xenograft tumours. First, the time course of abemaciclib-mediated target engagement and cell cycle arrest was assessed after a single 45 or 90 mg kg−1 oral dose, with blood and tumour tissue collection occurring at 1, 6, 24, 36 and 48 h postdose. Second, the time course of vemurafenib-mediated biomarker modulation was evaluated at a single 15 mg kg−1 oral dose, where blood and tumour tissue were collected at 1, 2, 6, 12, 24 and 48 h postdose. Third, the dose response of vemurafenib was assessed after a single oral dose of 7.5, 30 or 60 mg kg−1, with blood and tumour tissue collection at 4, 8 and 24 h postdose (n=5 per time point). All tumour tissue was evaluated for the cell cycle biomarkers, pRb, TopoIIα, pHH3 and Total Rb. In addition, tumour tissue obtained in the vemurafenib studies was also assayed for phospho-(ser217/221)-mitogen-activated protein kinase kinase (pMEK), phospho-(thr202/tyr204)-ERK (pERK) and Cyclin D1.

Tumour growth inhibition in A375 xenograft tumours

First, the in vivo efficacy of the abemaciclib and vemurafenib combination was evaluated. For 21 days, mice bearing A375 xenograft tumours received either the combination or monotherapy of vemurafenib two times a day at 10 mg kg−1 and abemaciclib daily at 45 mg kg−1 (n=8 per dose level). Treated tumour growth was compared with tumour growth in mice receiving one of three vehicle formulations (n=24). Tumour size was measured every 3–4 days for ∼6 weeks after the first dose.

Second, the in vivo efficacy of abemaciclib in mice bearing vemurafenib-resistant A375 xenograft tumours was evaluated. For 48 days, mice received vemurafenib at 15 mg kg−1 two times a day to allow resistance to develop; mice then received either 90 mg kg−1 abemaciclib daily or continued on vemurafenib therapy for 28 days (n=8 per group). An additional eight mice received a vehicle formulation as a control. Tumour size was measured every 3–4 days for ∼13 weeks after the first dose.

All animal studies were performed in accordance with American Association for Laboratory Animal Care institutional guidelines, and all protocols were approved by the Eli Lilly and Company Institutional Animal Care and Use Committee.

Development of the integrated semimechanistic PK/PD model

While the melanoma cell line, A375, was the focus of the present work, there is wealth of data available for the preclinical tool cell line, colo-205 (Gelbert ; Tate ; Yadav ). As the existing PK/PD model was shown to accurately predict biomarker response to abemaciclib in either cell line (Tate ), this allowed the authors to capitalise on the available data and develop the cell cycle model using Total Rb response data obtained in colo-205. As such, the following describes initial model development using colo-205 biomarker data, before validated translation to A375 to incorporate biomarker response to vemurafenib and efficacy.

PK/biomarker model for abemaciclib

A semimechanistic PK/PD model was previously developed to describe abemaciclib target engagement and resulting cell cycle arrest in mice bearing colo-205 xenograft tumours (Tate ). Briefly, the model comprises four transit compartments, each corresponding to early G1, late G1, S and G2/M cell cycle phases in a sequential, open-loop manner. The biomarkers, pRb, TopoIIα and pHH3, represent the change in cell phase density in late G1, S and G2/M phases, respectively (Gurley ; DeCaprio ; Goswami ; Kitagawa ), where abemaciclib inhibits production of pRb and thus induces cell cycle arrest in the G1 phase (Gelbert ). This model was further developed to include inhibition of Total Rb caused by abemaciclib-mediated pRb inhibition, as a result of autoregulation of the E2F1-Rb pathway (Shan ) (see Equations 9, 10, 11, 12, 13). The model was refitted with all colo-205 biomarker data previously used for model development, with the addition of Total Rb data obtained after single or multiple doses of 50 mg kg−1 abemaciclib. The rate constants for the cell cycle model were then recalculated using cell phase density (Huang ; Lam ), as shown previously (Tate ), to predict abemaciclib-mediated inhibition of pRb, TopoIIα, pHH3 and Total Rb in A375 xenograft tumours.

PK model for vemurafenib

The plasma–concentration time course of vemurafenib in mice bearing xenograft tumours was described by a one-compartment model with saturable absorption and linear clearance. The PK model is given by the following differential equations:

where AGut is the amount of drug in the gut, CVemurafenib is the concentration of vemurafenib in the plasma, Vmax and Km describe the saturable absorption, CL denotes drug clearance from plasma and V is the volume of distribution. Initial conditions for each PK model compartment were set to zero, except AGut, which was set equal to the dose.

Biomarker model for abemaciclib and vemurafenib

The biomarker model for vemurafenib was developed using A375 xenograft tumour data to describe drug-induced changes in the MAPK pathway. An indirect response model was used to relate vemurafenib plasma concentrations to inhibition of pMEK (Yang ), which was in turn directly related to inhibition of pERK (Dhillon ; Yang ), where IC50,Vemurafenib is the potency of vemurafenib and kout,pMEK is the turnover rate constant:

A series of transit compartments then describe the impact of pERK inhibition on Cyclin D1 (Terada ; Stacey, 2003) and the downstream effect of Cyclin D1 inhibition on the cell cycle (Stacey, 2003), where kout,Cyclin D1 is the Cyclin D1 turnover rate constant:

Three transit compartments were included to account for the apparent delay between vemurafenib-mediated inhibition of Cyclin D1 and pRb inhibition:

Inhibition of pRb both directly by abemaciclib (Gelbert ) and indirectly by vemurafenib (Stacey, 2003; Tap ) was incorporated into the cell cycle model, where Imax is the maximum drug effect and IC50,Abemaciclib is the drug potency:

The inhibitory effect of pRb downregulation on the expression of Total Rb, caused by autoregulation of the Rb pathway through E2F-1 (Shan ), was included in the biomarker model, where kout,Total Rb is the Total Rb turnover rate constant:

Each zero-order input in the above equations is defined as the product of the turnover rate (plus elimination rate, where relevant) and the baseline level of each biomarker:

As described previously (Tate ), the rate constants, kR, kG1S, kSG2 and kMG1, drive the rate of transition of the cells through each of the cell cycle phases. The initial phase distributions of the colo-205 and A375 cell lines (Huang ; Lam ) were used to adjust baseline levels of pRb, TopoIIα and pHH3 and to calculate the rate constants:

The initial conditions for all other biomarkers, pMEK, pERK, Cyclin D1, Total Rb and all associated transit compartments, were set equal to the baseline level of 100, representing the percentage of control.

Effect of emerging resistance on the biomarker model

The emerging resistance of the tumour to vemurafenib by the upregulation of the MAPK pathway (Das Thakur ) was incorporated into the biomarker model by use of a modulator compartment. In the model, the chronic presence of vemurafenib drives upregulation of the modulator compartment, which in turn increases production of pMEK, pERK and Cyclin D1. Although inhibition of the MAPK pathway by vemurafenib still occurs, the peak inhibition level no longer falls below non-resistant baseline levels of each biomarker, in accordance with previous experimental observations (Das Thakur ). As a result, vemurafenib-mediated tumour shrinkage no longer occurs, and tumour growth resumes.

The modulator compartment (Mod) and its effect on pMEK is described as follows, where Scalar is the magnitude of vemurafenib effect on the modulator compartment and kout,Mod is the turnover rate:

The zero-order input for the modulator compartment is defined below:

Tumour model

The growth rate of A375 xenograft tumours in control mice was described using the model previously developed by Simeoni . Tumour growth was affected in one of three ways: (1) by inhibition of the cell cycle (EpHH3) (Tate ), (2) by apoptosis triggered by the downstream effect of pERK inhibition (EpERK) (Mandal ) and (3) by increasing dependence on Total Rb for survival under vemurafenib-resistant conditions (ETotal Rb) (Yadav )

The exponential and linear growth parameters are denoted by λ0 and λ1, and the cytotoxic effect rate constants are denoted by kdeath,pERK and kdeath,Total Rb. The effect of reduced mitosis on tumour growth, EpHH3, is described below, where γ is fixed to 0.241 (Tate ):

The apoptotic effect of pERK inhibition, EpERK, is given by:

As the resistance to vemurafenib therapy develops, the tumour becomes increasingly dependent on Total Rb expression for survival (Yadav ). The apoptotic effect of Total Rb under resistant conditions was thus incorporated as the product of Total Rb and Cyclin D1 expression (ETotalRb). Therefore, when Cyclin D1 is overexpressed as a result of resistance to vemurafenib (Smalley ; Yadav ), administration of abemaciclib causes cell apoptosis and tumour shrinkage (Yadav ). To prevent Total Rb-mediated tumour shrinkage occurring under non-resistant conditions, ETotalRb is set equal to zero if either Total Rb is above baseline or if Cyclin D1 is below baseline.

A schematic representation of the biomarker/tumour growth model is provided in Figure 1, showing the effect of vemurafenib under non-resistant conditions (A), and the acquisition of BRAF inhibitor resistance and efficacy mediated by abemaciclib under resistant conditions (B).

Figure 1

Schematic representation of the integrated semimechanistic PK/PD model describing (A) vemurafenib-meditated efficacy under non-resistant conditions through inhibition of the MAPK pathway and the cell cycle, and (B) acquisition of resistance to vemurafenib by the upregulation of the MAPK pathway and abemaciclib-mediated efficacy through cell cycle arrest and increased dependence on Total Rb for survival. Arrowheads denote positive relationships, flatheads denote negative relationships. Red arrows indicate route of efficacy; purple arrows indicate elevation of baseline biomarker levels as a result of vemurafenib resistance. Shaded boxes denote observed compartments. Parameters are defined in the text.

Model implementation and validation

The data were analysed in a sequential manner using NONMEM VII (ICON Development Solutions, Ellicott City, MD, USA). A population approach was used where longitudinal data were available. Model selection was based on goodness of fit and diagnostic plots, parameter estimate precision, the Akaike Information Criterion value and the Objective Function Value. Internal validation of the models was performed by graphical comparison of the raw data used to develop the model with the 5th, 50th and 95th percentiles of 1000 model simulations (visual predictive check (VPC)). In addition, external validation was performed by VPC of the model vs experimental data not included in the fitting process, where available.

Results

Biomarker response to abemaciclib in colo-205 and A375 xenograft tumour tissue

The cell cycle model previously developed for colo-205 xenograft tumours (Tate ) was extended to include downregulation of Total Rb as a result of abemaciclib-mediated pRb inhibition (Shan ). The refitted model successfully described the slow turnover rate of Total Rb and its modest inhibition after single and daily oral doses of 50 mg kg−1 abemaciclib (Supplementary Figure 1). The model parameters related to pRb, TopoIIα and pHH3 remain largely unchanged from the previous modelling effort (Tate ), and as such, the model also successfully describes the concentration-dependent time course of pRb, TopoIIα and pHH3 (Table 1 and Supplementary Figure 1, additional data not shown).

Table 1

Parameter estimates for the integrated semimechanistic PK/PD model for vemurafenib and abemaciclib in mice bearing colo-205 or A375 xenograft tumours

Parameter	Colo-205 xenograft tumours	A375 xenograft tumours (model estimation)	A375 xenograft tumours (external validation)
Biomarker model
Imax	0.907 (1.7)	*	*
IC50,Abemaciclib (ng ml−1)	4.99 (21)	*	*
IC50,Vemurafenib (ng ml−1)	—	16 800 (7.7)	**
kout,pMEK (h−1)	—	1.76 (31)	**
kout,Cyclin D1 (h−1)	—	0.175 (5.9)	**
γCyclin D1	—	5 (fixed)	**
kR (h−1)	0.154 (6.3)	*	*
kG1S (h−1)	0.631a	0.618a	**
kSG2 (h−1)	0.0947a	0.445a	**
kMG1 (h−1)	0.272a	1.98a	**
kel (h−1)	0.0463 (29)	*	*
γpRb	0.458 (16)	*	*
kout,Total Rb (h−1)	0.0524 (38)	*	*
Tumour model
λ0 (h−1)	—	0.00290 (9.2)	0.00350 (5.4)
λ1 (mg h−1)	—	5.16 (41)	6.31 (14)
ω0 (mg)	—	26.6 (19)	114 (2.9)
kdeath,pERK (h−1)	—	0.0128 (29)	**
kout,modulator (h−1)	—	0.00191 (64)	**
Scalar	—	0.0000616 (26)	**
kdeath,Rb (h−1)	—	0.00990 (41)	**

For all fitted parameters, standard error of the estimate (%) is given within parentheses.

Estimates of interindividual variability and residual error are provided in Supplementary Table 3.

—, Not used in estimation or simulation.

*Simulation performed using parameter estimate obtained from colo-205 xenograft tumour data.

**Simulation performed using parameter estimate obtained from A375 xenograft tumour data.

Calculated from literature baseline values of cell density in the respective cell line.

The model was validated through accurate prediction of PK (data not shown) and biomarker response after 45 or 90 mg kg−1 of abemaciclib in A375 xenograft tumours by altering the initial conditions of the cell cycle biomarker model, and assuming equivalent potency to colo-205 xenograft tumours (Supplementary Figure 2), as shown previously (Tate ).

PK and biomarker response to vemurafenib in A375 xenograft tumour tissue

The disposition of orally administered vemurafenib in xenograft tumour-bearing mice was best described by a one-compartment PK model with saturable absorption, with a Michaelis constant (Km) corresponding to a 19.9 mg kg−1 dose (Table 1, Supplementary Figure 3 and Supplementary Table 2). Notably, neither vemurafenib nor abemaciclib exhibited altered pharmacokinetics after coadministration (data not shown).

The PK model was then linked to a biomarker model describing the concentration-dependent effects of vemurafenib on pMEK and pERK (Dhillon , Yang ), and the resulting inhibition of Cyclin D1 formation (Terada ; Stacey, 2003) (Figure 2, Supplementary Figures 4–6). Inhibition of MEK phosphorylation by vemurafenib was potent and rapid, with an IC50,Vemurafenib of 16.8 μg ml−1 (Table 1). It was assumed (through model testing) that vemurafenib was able to completely block pMEK synthesis (i.e. Imax=1). Although pERK is downstream of pMEK, the effect of vemurafenib-mediated inhibition of pMEK on pERK was found to be instantaneous. Cyclin D1 exhibits a much longer turnover half-life (4 h), resulting in maximum inhibition at 10 h after a 15 mg kg−1 dose of vemurafenib, compared with 2 h for upstream markers, pMEK and pERK (Figure 2).

Figure 2

Visual predictive check of the biomarker model in A375 xenograft-bearing mice following a single 15 mg kg−1 oral dose of vemurafenib. The circles denote observed pMEK, pERK, Cyclin D1, pRb, TopoIIα, pHH3 and Total Rb data in treated tumours, expressed as a percentage of the control value observed in the vehicle group. The solid and dotted lines represent the median, and the 5th and 95th percentiles of 1000 individual model simulations, respectively.

By relating Cyclin D1 to the pRb, the effect of vemurafenib-mediated Cyclin D1 inhibition on cell cycle arrest (i.e. pRb, TopoIIα, pHH3 and Total Rb) (Stacey, 2003) was also captured (Table 1, Figure 2, Supplementary Figures 4–6). The model indicates that cell cycle progression is highly sensitive to Cyclin D1 expression; the same vemurafenib dose of 15 mg kg−1 inhibits Cyclin D1 by 33%, which in turn inhibits pRb by 50% (Figure 2). Although initial inhibition of pHH3 appears modest (Figure 2), the model predicts that 15 mg kg−1 every 12 h inhibits pHH3 expression by 90% (under non-resistant conditions).

Inhibition of A375 xenograft tumour growth by vemurafenib acquisition of its resistance and subsequent efficacy by abemaciclib

The untreated growth of vehicle-treated A375 xenograft tumours was best described by the model previously developed by Simeoni (Figure 3A and Table 1). Vemurafenib-mediated tumour growth inhibition under non-resistant conditions was successfully captured by a model structure incorporating both the apoptotic effect of pERK inhibition (Mandal ) and the cell growth inhibitory effect of pHH3 inhibition (Tate ). This model indicates that 98% of antitumour activity is driven by pERK-related apoptosis, whereas the contribution of Cyclin D1-related cell cycle arrest is minimal (2%).

Figure 3

Visual predictive check of the tumour growth model in A375 xenograft-bearing mice for two studies: (A) after treatment with vehicle (top left) or with vemurafenib 15 mg kg−1 twice a day for 48 days, subsequently continuing with vemurafenib treatment (top right) or switching to 90 mg kg−1 abemaciclib every day for a further 28 days (bottom centre); and (B) after treatment with vehicle (top left) or with vemurafenib 10 mg kg−1 twice a day (top right) or abemaciclib 45 mg kg−1 every day alone (bottom left) or in combination for 21 days (bottom right). The circles denote observed A375 tumour size data. The solid and dotted lines represent the median, and the 5th and 95th percentiles of 1000 individual model simulations, respectively.

Acquisition of resistance as a result of continuous vemurafenib treatment was best described by incorporating an upstream modulator compartment to induce overexpression of pMEK, pERK and Cyclin D1 (Das Thakur ) (Figure 3A). Time to resistance was estimated to be ∼30 days, whereupon pERK and Cyclin D1 levels are elevated nearly twofold over baseline. Moreover, sensitivity to abemaciclib monotherapy during vemurafenib resistance was successfully described by conditioning tumour cell survival to Total Rb when Cyclin D1 is overexpressed. The model indicates that under resistant conditions, 68% of abemaciclib-mediated tumour shrinkage is driven by sensitivity to Total Rb, whereas the remaining 32% is attributed to cell cycle arrest.

Model validation: predicting tumour response in a second in vivo study

Model validation was performed by reestimating the growth of untreated A375 xenograft tumours in a second in vivo study and simulating the drug effect (Tate ) elicited by vemurafenib and abemaciclib alone and in combination (Figure 3B and Table 1). As with the first in vivo study, untreated tumour growth was best described by the model previously developed by Simeoni . The inhibitory effect of vemurafenib (10 mg kg−1 twice a day) and abemaciclib (45 mg kg−1 every day) alone or in combination for 21 days was successfully predicted using the integrated semimechanistic tumour model. Owing to the short duration of the study, resistance to vemurafenib was neither present in the observed data nor in the model prediction.

Model predictions: identifying an optimum dosing schedule for the combination of abemaciclib and vemurafenib

Simulations of the model were performed to assess the optimum dosing schedule for vemurafenib and abemaciclib, providing insight into potentially efficacious clinical scenarios. Dosing regimens were constrained to those considered most clinically comparable (15 mg kg−1 twice a day for vemurafenib and 50 mg kg−1 every day for abemaciclib) (Grippo ; Tate ), whereas scheduling was varied across a range of monotherapy and combotherapy, and continuous and intermittent options. In accordance with experimental data (Das Thakur ; Yadav ), an intermittent schedule of vemurafenib (15 mg kg−1 twice a day for 2 weeks on a 3-week cycle) was predicted to provide a more durable tumour response over continuous therapy (15 mg kg−1 twice a day) (Figure 4C vs Figure 4A). However, the onset of resistance was only delayed by intermittent vemurafenib therapy, eventually leading to tumour regrowth. The combination of continuous abemaciclib (50 mg kg−1 every day) with continuous vemurafenib (15 mg kg−1 twice a day) was also predicted to provide superior benefit over continuous vemurafenib alone (Figure 4B vs Figure 4A), in accordance with previous observations (Yadav ). However, the magnitude of response was comparable to that obtained with intermittent vemurafenib monotherapy, whereas resistance occurred more rapidly. Finally, continuous abemaciclib (50 mg kg−1 every day) combined with intermittent vemurafenib (15 mg kg−1 twice a day for 2 weeks on a 3-week cycle) was predicted to maximise tumour response over either intermittent vemurafenib monotherapy or the combination of abemaciclib with continuous vemurafenib (Figure 4D). This optimised combination was predicted by the model to achieve complete tumour regression and prevent tumour regrowth throughout the duration of treatment.

Figure 4

Simulations of the PK/biomarker/tumour growth model for four long-term dosing scenarios: (A) vemurafenib 15 mg kg−1 twice a day given continuously, (B) vemurafenib 15 mg kg−1 twice a day and abemaciclib 50 mg kg−1 every day in combination, both given continuously, (C) vemurafenib 15 mg kg−1 twice a day given on an intermittent 2 weeks on/1 week off schedule, (D) intermittent vemurafenib 15 mg kg−1 twice a day in combination with continuous abemaciclib 50 mg kg−1 every day.

Discussion

In this study, an integrated semimechanistic PK/PD model was established describing efficacy and developing resistance to vemurafenib in A375 xenograft tumours, and heightened sensitivity to abemaciclib therapy arising as a result of continuous vemurafenib treatment. The purpose of developing such a model was (1) to characterise the quantitative pharmacology of abemaciclib-mediated efficacy in vemurafenib-resistant BRAF-mutated melanoma, and (2) to identify an optimal dosing strategy for the abemaciclib/vemurafenib combination. To our knowledge, this is the first PK/PD modelling approach to mechanistically evaluate dose schedule optimisation for combination therapies in oncology.

A semimechanistic PK/PD model for abemaciclib in colo-205 xenograft tumours was previously shown to predict biomarker modulation (pRb, TopoIIα and pHH3) and efficacy in A375 xenograft tumours (Tate ). To describe abemaciclib-mediated inhibition of Total Rb, thought to be the mechanistic basis for sensitivity to CDK4/6 inhibition in vemurafenib-resistant tumours (Yadav ), this PK/PD model was expanded in the current work to describe the downregulation of Total Rb as a result of pRb inhibition, caused by the inherent autoregulation of the E2F1-Rb pathway (Shan ). To exploit the wealth of in house biomarker data available in colo-205 xenograft tumours (Gelbert ; Tate ; Yadav ), the dynamics of Total Rb inhibition were estimated using colo-205 data and, as before (Tate ), the model was translated to A375 xenograft tumours by recalibration of the cell cycle. An updated, more intensively sampled data set provided a further and more robust validation analysis for biomarker response to abemaciclib in A375 xenograft tumours (Supplementary Figure 2).

Once the biomarker model for abemaciclib was successfully updated, the semimechanistic PK/PD model was further expanded to successfully describe vemurafenib-mediated inhibition of the MAPK pathway modulation (pMEK, pERK and Cyclin D1) and cell cycle arrest (pRb, TopoIIα, pHH3 and Total Rb). Tumour shrinkage was related to inhibition of pERK and pHH3, each describing efficacy as a result of apoptosis or cell cycle arrest, respectively (Lee ). Model simulations reveal the majority of vemurafenib-mediated efficacy under non-resistant conditions is through pERK-related apoptosis (98%), whereas cell cycle arrest has only a minimal effect on tumour growth inhibition. Internal–external validation of the model was achieved using a previously established technique (Tate ), whereby uncontrolled tumour growth is refitted and the extent of tumour shrinkage as a result of drug treatment is accurately predicted (Figure 3B).

Although vemurafenib is initially efficacious, continuous treatment leads to resistance and tumour regrowth as a result of MAPK pathway reactivation (Trunzer ). Resistance to vemurafenib was incorporated into the PK/PD model by use of a modulator compartment to mechanistically drive emerging upregulation of the MAPK pathway markers, pMEK, pERK and Cyclin D1. This modulator compartment can be considered an empirical approximation of the underlying genetic processes leading to resistance, which for A375 cells includes FGF-mediated ERK activation and/or p61-BRAF splicing (Yadav ). As modulator upregulation is related to vemurafenib concentrations, resistance remains until cessation of vemurafenib treatment, at which point the biomarkers of the MAPK pathway gradually return to their original baseline levels and resistance is effectively reversed. This model structure is based on the growing evidence that vemurafenib resistance is reversible: altered EGFR expression levels as a result of vemurafenib treatment are counterselected in its absence (Sun ), intermittent dosing delays tumour progression in preclinical models of BRAF-mutated melanoma (Das Thakur ; Abdel-Wahab ) and patients who have previously developed resistance to treatment have been successfully rechallenged with BRAF inhibition (Seghers ; Romano ). Most notably, resistance-related overexpression of pERK returns to near-normal levels during a treatment holiday (Das Thakur ), with a recovery rate that appears comparable to model estimates (half-life ∼6 days). Consequently, by incorporating resistance to vemurafenib as a dynamic, reversible phenomenon, the PK/PD model was able to predict the additional benefit mediated by intermittent therapy over continuous treatment (Figure 4C). However, research in this area remains limited, and further work is required to investigate the apparent loss of resistance upon cessation of vemurafenib treatment and the clinical utility of an intermittent dosing schedule.

The CDK4/6 inhibitor, abemaciclib, was recently found to elicit a level of efficacy in vemurafenib-resistant A375 xenograft tumours, which was far greater than would be expected in treatment-naïve tumours (Yadav ). This amplified sensitivity to abemaciclib is thought to be caused by an increased dependency on Rb expression for cell survival (Yadav ). Resistant cells were also shown to be sensitive to Cyclin D1 expression, which is substantially elevated after continuous exposure to vemurafenib (Smalley ; Yadav ). As such, Cyclin D1 was incorporated as a marker of resistance in the PK/PD model, thus conditioning the model to increased sensitivity to inhibition of Total Rb. Despite only moderate levels of inhibition (∼40% at 50 mg kg−1 every day), Total Rb is sufficiently reduced by abemaciclib to induce apoptosis in vemurafenib-resistant tumours. Indeed, while abemaciclib efficacy in non-resistant tumours is wholly attributed to cell cycle arrest (50 mg kg−1 every day) (Tate ), the majority of tumour shrinkage in resistant tumours is mediated through sensitivity to Total Rb inhibition (68%), rendering cell cycle arrest a minor route (32%). Various dosing scenarios were explored using the PK/PD model; the strategy whereby vemurafenib is dosed intermittently during continuous abemaciclib treatment is predicted to be more effective than a continuous combination. This approach synergistically exploits both increased sensitivity to abemaciclib under vemurafenib-resistant conditions and improved sensitivity to BRAF inhibition after its temporary withdrawal. While the treatment schedule investigated in the current modelling work involves more frequent vemurafenib dosing (21-day cycle; 2 weeks on, 1 week off) than previous in vivo studies (42-day cycle; 4 weeks on, 2 weeks off; Das Thakur ), the dosing intensity remains comparable (2 : 1). The higher frequency of treatment in the current study was designed to be amenable to common clinical practice; model simulations of the 42-day schedule are only marginally less effective than the 21-day cycle (data not shown). Furthermore, reducing the dosing intensity of vemurafenib (e.g. 1 week on, 2 weeks off) did not considerably worsen the predicted level of efficacy (data not shown), and may offer an opportunity to lessen the side effects associated with BRAF inhibition.

As the clinical doses for each inhibitor are well established, the PK/PD model was designed to investigate optimum dosing schedules rather than to explore alternative dose levels. As such, the in vivo efficacy data for vemurafenib covered a limited dose range of 10–15 mg kg−1 twice a day. Although these doses exhibit exposure levels below clinical concentrations (Grippo ), the sensitivity of A375 xenograft tumours to vemurafenib requires the use of lower doses (Wong ). In contrast, the 50 mg kg−1 every day dose of abemaciclib in mouse achieves exposure levels that have been observed in patients demonstrating clinical efficacy (Tate ). It is therefore considered that, within the constraints of xenograft sensitivity, the doses used in the current work are relevant to current clinical practice.

To conclude, previous xenograft studies have demonstrated the additional benefit of abemaciclib either in combination with vemurafenib in treatment-naïve human melanoma tumours or as a secondary monotherapy treatment in tumours that have already developed resistance to vemurafenib (Yadav ; Yadav ). In the current work, a PK/PD model developed to optimise the treatment schedules indicates that intermittent therapy with vemurafenib in combination with continuous abemaciclib treatment could lead to complete response in the median tumour volume (Figure 4D). This combination scenario allows resistance to wane during the break in vemurafenib dosing, whereas abemaciclib yields additional efficacy in resistant cells through inhibition of Total Rb (Figure 4D). Additional in vivo studies are required to fully assess the extent of tumour shrinkage mediated by intermittent vemurafenib and continuous abemaciclib, and to further validate the predictive ability of the combination PK/PD model.