A Translational Quantitative Systems Pharmacology Model for CD3 Bispecific Molecules: Application to Quantify T Cell-Mediated Tumor Cell Killing by P-Cadherin LP DART®.

CD3 bispecific antibody constructs recruit cytolytic T cells to kill tumor cells, offering a potent approach to treat cancer. T cell activation is driven by the formation of a trimolecular complex (trimer) between drugs, T cells, and tumor cells, mimicking an immune synapse. A translational quantitative systems pharmacology (QSP) model is proposed for CD3 bispecific molecules capable of predicting trimer concentration and linking it to tumor cell killing. The model was used to quantify the pharmacokinetic (PK)/pharmacodynamic (PD) relationship of a CD3 bispecific targeting P-cadherin (PF-06671008). It describes the disposition of PF-06671008 in the central compartment and tumor in mouse xenograft models, including binding to target and T cells in the tumor to form the trimer. The model incorporates T cell distribution to the tumor, proliferation, and contraction. PK/PD parameters were estimated for PF-06671008 and a tumor stasis concentration (TSC) was calculated as an estimate of minimum efficacious trimer concentration. TSC values ranged from 0.0092 to 0.064 pM across mouse tumor models. The model was translated to the clinic and used to predict the disposition of PF-06671008 in patients, including the impact of binding to soluble P-cadherin. The predicted terminal half-life of PF-06671008 in the clinic was approximately 1 day, and P-cadherin expression and number of T cells in the tumor were shown to be sensitive parameters impacting clinical efficacy. A translational QSP model is presented for CD3 bispecific molecules, which integrates in silico, in vitro and in vivo data in a mechanistic framework, to quantify and predict efficacy across species.

Introduction

Immunotherapy, which recruits a patient’s own immune system to kill cancer cells, has begun to revolutionize cancer treatment (1). Within the class of immune-oncology therapies are the bispecific immune cell re-targeting molecules (2). These are typically recombinant bispecific antibodies, or antibody fragments, with one binding domain targeting a specific tumor antigen of choice and the other domain targeting CD3 on T cells. Because CD3 serves as the signaling component of the T cell receptor (TCR) complex, these CD3 bispecific molecules enable T cells to circumvent the need for the interaction between the TCR and antigen presented by major histocompatibility complex (MHC) class I molecules. This expands the repertoire of T cells able to recognize the tumor and stimulate them to act as effector cells (3). Similar to the standard immune synapse formation, once a threshold of bispecific-mediated molecular interactions has been reached, CD3 signals the T cell to initiate a cytotoxic response toward the adjacent tumor cell expressing the specific antigen. Cytotoxicity is mediated by the release of cytotoxic granules containing perforin and granzymes by the T cell. Perforin is a pore-forming protein enabling entry of granzymes, and the granzymes trigger a caspase cascade that leads to apoptosis. Activation of T cells leads to the transient release of cytokines and T cell proliferation, recruitment, and infiltration into the tumor environment, which drives serial killing of tumor cells.

In 2014, blinatumomab (CD3-CD19) was the first CD3 bispecific construct approved in the USA for the treatment of resistant/refractory B cell acute lymphocytic leukemia (B-ALL) (4). Blinatumomab is also being investigated in phase 2 clinical trial in patients with resistant/refractory non-Hodgkin’s lymphoma (NHL) (5). The first generation bispecific T cell retargeting molecules such as blinatumomab are tandemly linked single-chain Fv (scFv) known as bispecific T cell engager (BiTE) molecules (2,3). These molecules are around 50 kDa and have a short circulating half-life (approx. 2 h) requiring constant infusion through the use of a pump to achieve a stable therapeutic exposure of the molecule (6). New generation CD3 bispecifics with a variety of formats are being tested in clinical trials. These include PF-06671008 which is a P-cadherin-specific LP DART: a molecule based on the DART® platform, but containing a human IgG1 Fc domain to extend the half-life (7). This bispecific targets CD3 and P-cadherin expressed on solid tumors. P-cadherin is a member of a family of molecules that mediate calcium-dependent cell-cell adhesion and has been reported to correlate with increased tumor cell motility and invasiveness when over-expressed (8–10). Upregulation of P-cadherin has been reported in breast, gastric, endometrial, colorectal, and pancreatic carcinomas and correlates with poor survival of breast cancer patients (11–14). In contrast, P-cadherin has low expression in normal tissues, making it an attractive target for immunotherapy (12). In preclinical studies, in vitro and in vivo data indicate that PF-06671008 is a highly potent molecule eliciting P-cadherin expression-dependent cytotoxic T cell activity across a range of tumor indications (15). In addition, PF-06671008 is stable and has desirable biophysical and PK properties with a half-life of 3.7–6 days in mouse (7,15). PF-06671008 is currently being investigated in phase 1 clinical trials in patients with advanced solid tumors with the potential to have P-cadherin expression (https://clinicaltrials.gov/ct2/show/NCT02659631).

In order to characterize the in vivo efficacy of PF-06671008 in tumor-bearing mice, a quantitative systems pharmacology (QSP) model was established. This model integrates the PK of PF-06671008, its binding to shed P-cadherin and circulating T cells in the systemic circulation, its biodisposition in the tumor and the formation of a trimolecular complex (trimer) with T cells, and P-cadherin expressing tumor cells in the tumor microenvironment (TME). The model incorporates T cell kinetics in the tumor including T cell proliferation and contraction. The concentration of the trimer in the tumor is used to drive efficacy in mouse using an optimized transduction model of tumor cell growth and killing. In this manuscript, we discuss the use of the model to characterize the underlying pharmacology in mouse, and translation of the preclinical efficacy data to the clinic by incorporation of predicted human PK and disease parameters. The quantitative translational framework for CD3 bispecific molecules presented here can aid in drug design, candidate selection, and clinical dosing regimen projection.

Materials and Methods

In Vivo Studies

All procedures in animals were approved by the Pfizer Institutional Animal Care and Use Committees and studies were performed according to established guidelines.

PF-06671008 Mouse PK Study

PF-06671008 was administered as a single intravenous (IV) dose of 0.05 or 0.5 mg/kg to HCT-116 tumor-bearing female NOD-scid IL-2rgnull (NSG) mice, (n = 3/time point/dose) with or without human peripheral blood mononuclear cell (PBMC) engraftment. Mice were injected with 5 × 106 HCT-116 cells in Matrigel subcutaneously in the dorsal left flank. When the tumors had grown to approximately 0.5 g in size (after 14 days), the mice were administered PF-06671008. Serum and tumor samples were collected at predetermined time points from 5 min to 240 h post-dose.

ELISA Assay to Quantify PF-06671008

PF-06671008 concentrations in mouse serum and tumor homogenate were determined using an enzyme-linked immunosorbant assay (96-well format) with colorimetric detection. Briefly, the capture protein was a polyclonal goat antibody recognizing the CD3 scFv domain and the detection antibody was a goat anti-human IgG-Biotin (Qualex), followed by HRP-streptavidin conjugate (Jackson ImmunoResearch, West Grove, PA). Optical density was measured on a spectrophotometer (molecular devices). The lower limit of quantitation (LLOQ) of the assay was 12.5 ng/mL for serum samples and 1.5 ng/mL for tumor samples. The minimal required dilution was 1:25 for serum and 1:6 for tumor.

Flow Cytometric Tumor-Infiltrating Lymphocyte Analysis

HCT-116 tumor-bearing mice (n = 3) engrafted with human PBMC and administered a single IV dose of 0.01, 0.05, or 0.5 mg/kg PF-06671008 were euthanized pre-dose and 24, 72, and 144 h following dosing to assess tumor-infiltrating human CD3+ lymphocytes. Tumor samples were collected into gentleMACS C tubes containing human tumor cell dissociation buffer (Miltenyi Biotech) and processed to single-cell suspensions using the manufacturer’s suggested protocol for soft human tumors using the gentleMACS tissue dissociator (Miltenyi Biotech). After subsequent washing steps and live cell counting (using a hemocytometer and trypan blue exclusion), 1 × 106 live cells from each sample were collected and stained with CD3 FITC (BD Pharmingen) for 30 min on ice. Samples were analyzed using LSRII with FACS Diva software (BD Pharmingen). Absolute numbers of CD3+ T cells per gram of tumor were then calculated using the number of CD3+ events and sample tumor weight.

PF-06671008 Mouse Xenograft Studies

Mouse xenograft studies were completed in human T cell engrafted (HCT-116) or adoptive transfer (HCT-116 or SUM-149) established tumor models. In the human T cell engrafted model, tumor cells (5 × 106 HCT-116) were implanted subcutaneously (SC) into the right flank of 6–8-week-old female NSG mice as a 0.2 mL bolus mixed with 4 mg/mL Cultrex basement membrane extract (Trevigen) in PBS. Seven days prior to randomization, mice were inoculated with 5 × 106 or 2.5 × 106 freshly isolated human PBMC as an intraperitoneal injection of 0.2 mL cell suspension in PBS. In addition to vehicle (PBS), dose levels of 0.01, 0.05, 0.1, 0.15, and 0.5 mg/kg PF-06671008 were administered for HCT-116 studies (n = 10/dose). The doses were administered IV as a q7d × 2 regimen.

For the T cell adoptive transfer established tumor model, 8- to 10-week-old NSG mice were inoculated with either 5 × 106 HCT-116 cells in the flank or 5 × 106 SUM-149 cells in the mammary fat pad in a total injection volume of 0.2 mL, 7 days prior to randomization. HCT-116 cells were suspended in PBS, while SUM-149 cells were suspended in growth media and mixed 1:1 with Matrigel basement membrane matrix (BD Biosciences, San Jose, CA). T cells, which had been isolated from PBMCs, were activated and expanded using Dynabeads Human T-Expander CD3/CD28 magnetic beads (Life Technologies) for 6–9 days, depending on the study, were harvested, and resuspended in PBS at 1 × 107 cells/mL for in vivo inoculation. An initial dose of PF-06671008 or vehicle was administered to mice on day 0 and on the following day, mice were inoculated with 0.5, 1, 2, 2.5, or 5 × 106 T cells/mice IV. In addition to vehicle, dose levels of 0.05, 0.15, and 0.5 mg/kg PF-06671008 were administered for HCT116 xenograft studies and 0.05, 0.15, and 0.5 mg/kg PF-06671008 for SUM149 xenograft studies (n = 10/dose). The doses were administered IV as a q7d × 3 or q7d × 5 regimen.

Tumor volume was measured using a digital Vernier caliper (Mitutoyo America, Aurora, IL), and volumes were calculated by the use of the modified ellipsoid formula ½ (width2 × length). Tumor measurements were collected twice weekly, with continuous health monitoring, until the animals had to be euthanized due to tumor burden or health concerns out to a maximum of 16 days (HCT-116 in the human T cell engrafted model), 42 days (SUM-149 adoptive transfer), or 65 days (HCT-116 adoptive transfer). For a full description of the mouse xenograft studies, please refer to reference (15).

PF-06671008 Cynomolgus Monkey PK Study

The PK of PF-06671008 in cynomolgus monkey was evaluated following IV bolus and SC administration at weekly escalating doses for 1 month. The IV doses administered were 1.1/3.3, 3.3/10, or 10/20 μg/kg/week, and the SC dose administered was 10/30 μg/kg/week. This study has been described previously (16).

Measurement of Soluble P-Cadherin

Baseline sPcad levels were measured in cynomolgus monkey, healthy volunteer, and cancer patient serum samples (Bioreclamation). sPcad levels were also measured in cynomolgus monkey after treatment with PF-06671008 (in-house samples). A qualified Meso Scale Discovery (MSD) human P-cadherin kit was used to measure soluble P-cadherin levels, as described previously (16).

Modeling of Mouse Tumor Growth Inhibition Data

A QSP model was constructed to describe the disposition of PF-06671008 and T cells in the central compartment and tumor of the xenograft mouse models (Fig. 1a). The model accounts for the binding of PF-06671008 to tumor cells and T cells in the extracellular space of the TME to form trimers. The trimers are assumed to drive tumor cell killing. Description of all the symbols and parameters used in the mouse equations are shown in Tables I and II.

Fig. 1

a Translational quantitative systems pharmacology model for CD3 bispecific molecules. Parameter descriptions and values are summarized in Tables I to III. The figure represents both mouse and human models, with the following exceptions: binding to sPcad was only included in the human model and T cell proliferation/ exhaustion in the tumor were only included in the mouse model. b Schematic of the bell-shaped concentration relationship which can be observed for CD3 bispecific molecules. Formation of trimers between drugs, T cells, and tumor cells, is required for efficacy. The QSP model predicts trimer concentration and links it to tumor cell killing

Table I

Model Variables and Terms Used in Equations

Variable	Definition	Unit
TotPcad t	Total Pcad in the tumor	nM
TotCD3 p	Total CD3 in the central compartment	nM
TotCD3 t	Total CD3 in the tumor	nM
Tcellsp	T cells in the central compartment	Cells/L
Tcellstm	T cells migrated from plasma to tumor, during 5-day lag time	Cells/L
Tcellst	T cells in the tumor	Cells/L
Dose	Dose of PF-06671008	nmols
C1	Concentration of PF-06671008 in central compartment	nM
C2	Concentration of PF-06671008 in the peripheral compartment	nM
C3	Concentration of PF-06671008 in the tumor	nM
DCD3p	Dimer of PF-06671008-CD3 in the central compartment	nM
DCD3t	Dimer of PF-06671008-CD3 in the tumor	nM
DPcadp	Dimer of PF-06671008-Pcad in the central compartment	nM
DPcadt	Dimer of PF-06671008-Pcad in the tumor	nM
Trimer	Trimer of PF-06771008-CD3-Pcad in the tumor	nM
w	Total tumor volume	mm3
M1, M2, M3, M4	Tumor volume in growth and three transduction compartments	nM

Table II

Mouse Model Parameters

	Parameter	Definition	Unit	Value (CV%)			Source
Binding	konCD3, koffCD3, Kd_CD3	Binding of PF-06671008 to CD3	1/nM/h, 1/h nM	1.72, 19.66 11.4			(7)
	konPcad, koffPcad, Kd_Pcad	Binding of PF-06671008 to P-cadherin		1.57, 0.74 0.47
Central/peripheral compartment	V 1	Volume of distribution in the central compartment	mL/kg	49.6 (9)			Estimated from mouse PK data. kel = CL/V1 k12 = CLd/V1 k21 = Cld/V2
	V 2	Volume of distribution in the peripheral compartment	mL/kg	60.7 (16)
	CL	Clearance	mL/h/kg	0.45 (12)
	CL d	Inter-compartmental clearance	mL/h/kg	4.95 (28)
	Omega CL	Inter individual variability in clearance	–	0.064 (41)
	a	Additive error	–	0.067 (32)
	b	Proportional error	–	0.207 (15)
	Tcells p 0	Number of T cells administered/mouse	Cells/L	0, 2.5e8, 5e8, 1e9, 1.25e9, 2.5e9			See “Methods” section
	CD3	CD3 expression on T cells	Receptors/cell	100,000			(17,18)
Tumor disposition of PF-06671008/T cells	P	Permeability of drug into tumor	μm/day	334			(19)
	D	Diffusivity of drug into tumor	cm2/day	0.022
	ε	Void fraction in tumor for drug	–	0.24
	Rcap	Capillary radius	μm	8
	Rkrogh	Average distance between 2 capillaries	μm	75
	Tlag	Lag time for T cell disposition/onset of T cell proliferation in the tumor	Day	5			Set empirically, using in-house data
	kelT, k12T, k21T	T cell redistribution from the central compartment to the tumor	1/day	2.51, 0.002, 0.0005
Tumor compartment	P rate	T cell proliferation rate	1/h	Function of dose. Eq. (9)			See “Methods” section
	k _exhaust c	Slope of T cell decline	1/h	0.0412			Interpolated from (20)
	Tumor cellst	Number of tumor cells	Cells/g of tumor	1e8			(21)
	mPcad	P-cadherin expression on tumor cells	Receptors/cell	28,706 (HCT-116) 17,500 (SUM-149)			(15)
	Rtumor	Tumor radius	cm	Calculated from w			Measured
	k int	P-cadherin internalization rate	Day−1	0.1728 (−)			Estimated from mouse tumor PK data. Represents 96 h half-life of internalization.
Tumor growth inhibition		Mouse tumor models		HCT116 a	HCT116 b	SUM149 b
	k g0	Exponential tumor growth rate	1/day	0.30 (−)	0.19 (3)	0.12 (3)	Estimated in mouse models from unperturbed tumor growth data
	k g	Linear tumor growth rate	mm3/day	105 (4)	123 (2)	74.3 (5)
	M max	Maximum tumor volume	mm3	3.8 × 103 (−)	6.0 × 103 (−)	5.8 × 103 (−)
	ψ	Switch between exponential and linear growth phases	–	20 (−)			Fixed based on (22)
	k max	Maximum killing rate	1/day	0.74 (7)	1.32 (7)	2.71 (14)	Estimated in mouse models
	kC 50	Concentration at half maximum kill rate	nM	1.0 × 10−4 (6)	6.9 × 10−5 (7)	2.0 × 10−4 (15)
	τ	Transduction time between tumor compartments	Day	4.78 (10)	3.99 (1)	2.25 (3)
	Omega k g0	Inter-individual variability in exponential growth rate		0.46 (14)	0.34 (11)	0.12 (25)
	Omega k g	Inter-individual variability in linear growth rate		0.35 (13)	0.16 (13)	0.16 (28)
	a	Additive error		5 (−)	60 (−)	60 (−)
	b	Proportional error		0.26 (3)	0.06 (6)	0.01(50)
	TSC (pM) (80% confidence interval)	Tumor static concentration of the trimer		0.064 (0.044, 0.096)	0.011 (0.0096, 0.013)	0.0092 (0.0071, 0.012)

T cell engrafted tumor model

T cell adoptive transfer tumor model

Onset of exhaustion of T cells set to 7 days after disposition in the tumor

MWt of PF-06671008 = 105 kDa, MWt of sPcad = 85 kDa

Modeling of PF-06671008 and T cells in central/peripheral compartments and distribution to the tumor: Following systemic administration to mouse, PF-06671008 is assumed to be able to distribute to a peripheral compartment, distribute to the tumor, bind to circulating T cells, or be cleared from the central compartment. In the mouse model, PF-06671008 does not bind to sPcad.

Mouse PK: The mouse serum concentration profiles in human PBMC-engrafted mice, following IV administration of PF-06671008 at 0.05 and 0.5 mg/kg, were described using a two-compartment model with linear elimination from the central compartment (Eqs. 1 and 2). C1, C2, and C3 are the concentrations of the drug, PF-06671008, in plasma, peripheral compartment, and tumor, respectively. kel is the elimination rate of PF-06671008 from the central compartment. k12 and k21 are the inter-compartmental rate constants describing the distribution of PF-06671008 between the central compartment and the peripheral compartment. These values were fixed in the subsequent TGI modeling. Distribution of free PF-06671008 to the extracellular environment of the tumor was characterized using tumor disposition equations (Eqs. 3 and 4) that have been described previously (19,25,26). Briefly, P is the rate of permeability and D is the diffusion of the drug, across and around the tumor blood vessels. R is the radius of the tumor blood capillary, R is the average distance between 2 capillaries, R is the radius of the tumor, and ε is the tumor void volume for the drug.

Binding to T cells: Binding of PF-06671008 to circulating T cells was determined from CD3 binding (k and k), the number of CD3 receptors per T cell (CD3) and number of T cells in the central compartment, or plasma (Tcells). These values were used to calculate total CD3 in the central compartment (TotCD3) (Eq. 5). Binding to CD3 (and P-cadherin) was determined using surface plasmon resonance (SPR) assays run on a Biacore instrument as described previously (7). The number of CD3 receptors per T cell was taken from literature data (17,18). The number of T cells administered per mouse was used to inform the initial number of T cells in the central compartment. See Eq. 6 for binding of PF-06671008 to CD3 in the central compartment. DCD3p and DCD3t are the concentration of drug-CD3 dimers in plasma and tumor, respectively.

T cell trafficking: Following administration of T cells to the mouse, T cells were assumed to be able to distribute to the tumor (Eqs. 7 and 8), bind to PF-06671008, or be cleared from the central compartment. k and k are the rate constants describing the distribution of T cells between the central compartment and the tumor. k is the elimination rate of the T cells from the central compartment. These parameters were determined from modeling of in-house PBMC data in tumor-bearing mice (not shown here). A lag time of 5 days was introduced to accommodate the disposition and start of the proliferation of T cells at the tumor site. This was informed from in-house immunohistochemistry data and was equivalent to the time of T cell observation in the tumor. Tcells are T cells which have migrated from the central compartment to TME during the 5 days.

Modeling of T cell proliferation and trimer formation in the TME

Table III

Predicted Human Parameters Used in Simulations

	Parameter	Definition	Unit	Value (CV%)	Source
Binding	konCD3, koffCD3, Kd_CD3	Binding of PF-06671008 to CD3	1/nM/h, 1/h nM	1.72, 19.66, 11.4	(7)
	konPcad, koffPcad, Kd_Pcad	Binding of PF-06671008 to P-cadherin		1.57, 0.74, 0.47
Central/peripheral compartment	V 1	Volume of distribution in central compartment	mL/kg	40.2	Allometrically scaled from cynomolgus monkey PK analysis (16) kel = CL/V1 k12 = CLd/V1 k21 = Cld/V2
	V 2	Volume of distribution in peripheral compartment	mL/kg	211
	CL	Clearance	mL/h/kg	4.61
	CL D	Inter-compartmental clearance	mL/h/kg	25.2
	sPcad	sPcadherin concentration in central compartment	nM	1.1 (0.4–4.1)	Measured in-house (= 92.7 ng/mL) Median value of healthy subjects and patient data
	k deg	sPcad degradation rate	1/h	0.15	Allometrically scaled from cynomolgus monkey PK analysis (0.31 1/h in cyno) (16)
	k degcx	sPcad-PF-06671008 complex degradation rate	1/h	0.115	Assumed to equal PF-06671008 elimination rate (CL/V1)
	Tcells p	T cell concentration in central compartment	Cells/μL	5000	(23)
	CD3	CD3 expression on T cells	Receptors/cell	100,000	(17,18)
Tumor disposition of PF-06671008/T cells	P	Permeability of drug into tumor	μm/d	334	(19)
	D	Diffusivity of drug into tumor	cm2/d	0.022
	ε	Void fraction in tumor for drug	–	0.24
	Rcap	Capillary radius	μm	8
	Rkrogh	Average distance between 2 capillaries	μm	75
	Tcells t a	Number of T cells in tumor	Cells/g of tumor	6.49e5	(24)
	Tumor cellst	Number of Tumor cells	Cells/g of tumor	1e8	(21)
	mPcad	P-cadherin expression on tumor cells	Receptors/cell	28,706	(15)
	Rtumor	Tumor radius	cm	1	Assumed
	k int	Internalization rate with PBMCs	Day−1	0.1728 (−)	Estimated from mouse tumor PK data. Represents 96 h half-life of internalization.

Assume no proliferation in tumor; ksyn = kdeg*sPcad. MWt of PF-06671008 = 105 kDa, MWt of sPcad = 85 kDa

T cell kinetics in the tumor: CD3+ cells/mg tumor were measured in HCT-116 tumor-bearing mice, engrafted with human PBMCs, following administration of PF-06671008 at 0.01, 0.05, or 0.5 mg/kg (described above). This data was used to determine the proliferation rate of T cells in the tumor. The relationship between CD3+ cells/mg tumor with time at each dose level was described using an exponential function. The slope of each line represents the rate of proliferation of CD3+ cells and was plotted versus PF-06671008 dose. An empirical model was then used to describe the CD3+ proliferation rates (Prate) as a function of dose (Eq. 9). Please see Supplementary Material for additional information and plots. T cells migrating into the TME during the 5-day lag time (Tcells) undergo proliferation for 7 days (Eq. 10). Following proliferation, T cells undergo contraction which was characterized using mono-exponential decline (k) (Eq. 11). The time (7 days) and the rate of decline (0.0412 1/h) were estimated from literature data (20). It was assumed that T cell proliferation was only taking place in the tumor environment and that proliferation and contraction rates were the same in the human T cell engraftment and adoptive transfer mouse tumor models.

Trimer formation: In the TME, PF-06671008 can bind to P-cadherin on tumor cells or to CD3 on T cells to form dimers, or both tumor cells and T cells to form the active trimers. The binding constants between drug and P-cadherin are k and k and the binding constants between drug and CD3 are k and k. In addition to binding affinity values, trimer formation was a function of P-cadherin receptors per tumor cell (mPcad), number of tumor cells (Tumor), CD3 receptors per T cell (CD3), and number of T cells in the TME (Tcellst). These values were used to calculate total Pcad (TotPcad) and total CD3 (TotCD3) in the TME (Eqs. 12 and 13). P-cadherin receptor expression in HCT-116 and SUM-149 tumor cell lines was determined by phycoerythrin (PE) labeling of anti-P-cadherin mAb and flow cytometry to determine the number of PE-labeled antibodies bound per cell. This study has been described previously (15). Internalization rate of drug bound to P-cadherin (k) in the tumor was determined from the mouse PK study, completed in the presence of PBMCs. The number of tumor cells was determined from xenograft data. The number of CD3 receptors/T cell was taken from published data (17,18). See Eqs. 14–16 for binding of PF-06671008 to P-cadherin and CD3 to form dimers and trimers.

Tumor growth inhibition: The mouse xenograft PK/PD relationship was established by relating mouse PF-06671008 trimer concentration in the TME to measured xenograft tumor volume data using an optimized cell distribution transduction model (27). The presented model is a modified version of the model by Simeoni et al. (22). Briefly, the unperturbed tumor growth was fitted first using individual animal growth data from the vehicle control group, using a logistic model describing linear (k) and exponential (k) growth. The measured initial tumor volume in each animal was used to inform the initial conditions (M1). M1–M4 are the tumor volume in the growth compartment and three transduction compartments, respectively. w is the total tumor volume (mm3). The inter-individual variability of the growth parameters and the maximum tumor volume (M) obtained from the unperturbed growth model were then fixed in the simultaneous estimation of growth and drug effect parameters from the complete tumor volume data set. Tumor cell killing was driven by the concentration of the trimolecular complex (Trimer). τ is the transduction time, k is the maximum kill rate, kc is the concentration of the trimer in the tumor at half the maximal kill rate, and ψ is the constant for switching from exponential to linear growth patterns. Equations 17–22 describe the tumor growth inhibition modeling.

Determination of tumor static concentration (TSC): TSC is the concentration of trimers at which tumor growth and death rate are equal, and is defined as the minimal efficacious concentration (Ceff). This PK/PD-derived parameter combines growth information and drug effect, providing insight into the efficacy of PF-06671008 in mouse xenograft models. TSC was used as a translational factor for extrapolation of xenograft data to the clinic. See Eq. 23 for TSC calculation.

An 80% confidence interval on TSC was calculated using parametric bootstrap by resampling from the estimated parameters using a log-normal distribution.

Modeling

All modeling was performed using Monolix software v4.3.3 (Paris, France). The quality of the model fitting was assessed using the following:

Diagnostic plots: (a) plots of observations versus population/individual predictions and comparison with line of unity, (b) plots of weighted residuals versus time/concentration and check for systematic deviation from zero, (c) visual predictive checks of observations and predictions for all individuals at each dose level to check for goodness of fit.

Diagnostic criteria: (a) reasonable precision of the parameter estimates (RSE/CV%), (b) lack of correlation between model predicted parameters (< 0.95), (c) lack of shrinkage (η−) as a check for model over-parameterization (< 40%), (d) reduction in objective function values and/or Aikake and Schwarz criterion for model comparison.

Translation of the Model to Human

Prediction of Human PK

Human PK parameters were predicted from cynomolgus monkey PK parameters using a two-compartmental PK model which incorporates binding to sPcad (Table III and Fig. 1a). PK parameters were scaled from monkey to a human using allometric exponents of 0.9 for clearance, 1 for the volume of distribution, and − 0.25 for absorption rate. These exponents were selected as they have been previously identified as optimal for monoclonal antibodies (28). The degradation rate of sPcad (k) was scaled from monkey to human using an exponent of − 0.25. The degradation rate of the PF-06671008-sPcad complex (k) was assumed to be the same as PF-06671008 elimination rate.

Prediction of Clinical PK/PD

The QSP model used to describe the PK/PD relationship in the mouse was translated to a human using the physiological parameters and assumptions described in Table III. An important difference from the mouse model is that PF-06671008 binds to circulating target (sPcad) to form drug-P-cadherin (DPcadp) dimers in the central compartment in the human model. The additional model equations are shown in Eqs. 24–26. In addition, T cell proliferation/contraction kinetics were not included in the human model. Instead, a “steady-state” number of T cells in the tumor are assumed (T). All model simulations were completed using Berkeley-Madonna v8.3.18.

Sensitivity Analyses

Local sensitivity analyses were performed to assess the sensitivity of the QSP model to P-cadherin receptor expression on tumor cells, and to tumor T cell (effector) to tumor cell ratio (E:T), as these are potentially variable parameters in patients. P-cadherin receptor numbers of 1000, 3000, 10,000, and 28,706 were used for simulations with the human model. These values represented the range of P-cadherin expression measured across human tumor cell lines (15). The nominal value of E:T used in the model was 1:150, which is thought to be representative of a solid tumor (21,24). In the sensitivity analysis, E:T ratios of tenfold lower (1:15) and tenfold higher (1:1500) than the nominal value were investigated in the human model. For quantitative comparison, sensitivity was represented as predicted tumor trimer concentration at each expression level, or E:T ratio, following an IV dose of 0.1 μg/kg PF-06671008 QW to cancer patients.

Results

Serum and Tumor PK of PF-06671008 in Mouse

PK profiles of PF-06671008 in PBMC engrafted and non-PBMC engrafted HCT-116 tumor-bearing mice following single-dose IV administration at 0.05 and 0.5 mg/kg are shown in Fig. 2a. The area under the curve (AUC) of PF-06671008 in serum was dose proportional between 0.05 and 0.5 mg/kg and similar between PBMC engrafted and non-PBMC-engrafted mice (Supplemental Fig. 1). In contrast, the tumor AUC from the study with PBMC engraftment was more than fivefold higher than the study without PBMCs (Fig. 2b). This was attributed to a reduction in the internalization of PF-06671008 bound to P-cadherin on tumor cells in the presence of PBMCs.

Fig. 2

a Serum and b tumor PK profiles of PF-06671008 in PBMC engrafted and non-PBMC engrafted HCT-116 tumor-bearing mice following single-dose intravenous administration at 0.05 and 0.5 mg/kg

The serum PK in the PBMC-engrafted mice was used for PK modeling. The estimated serum PK parameter estimates for PF-06671008 are shown in Table II, and the goodness of fit plots are shown in Supplemental Fig. 1. The tumor internalization rate in the presence of PBMCs was used in the TGI PK/PD modeling (Table II).

Tumor T Cell Kinetics

HCT-116 tumor-bearing mice engrafted with PBMCs and administered PF-06671008 showed dose-dependent increases of tumor-infiltrating/proliferating CD3+ lymphocytes (TILs) over time (Fig. 3). The relationship with time was transformed to calculate a proliferation rate of CD3+ cells as a function of the dose which was used to describe tumor T cell kinetics in the QSP model (Supplementary Fig. 2).

Fig. 3

PF-06671008-induced tumor T cell proliferation in mice bearing HCT-116 tumors with human PBMC engraftment. Number of CD3+ cells/mg of the tumor (with standard deviations) is plotted against time following IV administration of control and PF-06671008 at 10 μg/kg, 50 μg/kg, and 500 μg/kg

PK/PD Relationship of PF-06671008 in Mouse Xenograft Models

The QSP model (Fig. 1a) was used to fit the tumor growth inhibition data obtained from the HCT-116 and SUM-149 mouse xenograft studies. The tumor trimer concentration was used as a driver of tumor cell killing. Estimated model parameters with percent coefficient of variation (CV) and calculated tumor trimer TSCs with 80% confidence intervals are shown in Table II. Parameters were estimated with good precision as assessed by the percentage CV for all cell lines. The goodness of fit and model performance were assessed using the goodness of fit plots (population prediction, individual prediction, and visual predictive check) that are shown in Supplemental Fig. 3 for HCT-116 in the T cell engrafted model and for HCT-116 and SUM-149 in T cell adoptive transfer experimental model. Overall, the median response and variability of all cell lines were described well by the mechanistic model. The calculated population median TSCs were 0.064, 0.011, and 0.0092 pM for HCT-116 in T cell engrafted model and for HCT-116 and SUM-149 in T cell adoptive transfer experimental models, respectively. The Ceff for tumor stasis is defined as the geometric mean of the TSCs in three mouse xenograft models and was calculated to be 0.028 pM trimer concentration in the tumor.

Serum P-Cadherin Concentrations Across Species

The concentrations of sPcad in serum samples from cynomolgus monkey, healthy humans, and cancer patients are shown in Table IV. There was no difference in sPcad levels in the serum of healthy human volunteers and cancer patients. Higher variability was observed in the lung and colorectal cancer samples compared to samples from breast cancer patients or healthy humans. Levels of sPcad in cynomolgus monkeys were similar to those in human.

Table IV

Concentration of Soluble P-Cadherin in Cynomolgus Monkey and Human Serum

Species	Disease state	n	Soluble P-cadherin concentration
			Median (ng/mL)	Range (ng/mL)
Cynomolgus monkeya	Healthy	32	47	29–273
Cynomolgus monkeyb	Healthy	4	68	57–74
Humana	Healthy	40	90	45–150
Humana	Breast cancer patients	23	78	32–190
Humana	Colon cancer patients	31	102	36–328
Humana	Lung cancer patients	25	102	65–320

Samples from Bioreclamation (Westbury, NY)

Samples from in-house studies

Clinical PK Predictions for PF-06671008

The predicted human PK parameters for PF-06671008 are shown in Table III. The predicted human CL and Vss were 4.6 mL/h/kg and 251 mL/kg, respectively, and the terminal half-life was predicted to be approximately 1 day.

Clinical PK/PD Predictions for PF-06671008 and Sensitivity to P-Cadherin Expression on Tumor Cells and T Cell Number

To translate the QSP model from mouse to human, the predicted human PK was incorporated along with assumptions and parameters describing the human physiology (Table III). The model simulated serum PK and tumor trimer concentrations following IV infusion of PF-06671008 at 0.01, 0.1, and 1 μg/kg QW to cancer patients are shown in Fig. 4 a and b, respectively. Expression levels of P-cadherin on tumor cells are expected to vary across patients. To investigate the potential impact on tumor trimer concentrations, a sensitivity analysis was performed varying P-cadherin receptor numbers from 1000 to 28,706 (HCT-116). Predicted tumor trimer concentration increases with increasing receptor expression (Fig. 5a) suggesting that P-cadherin expression is a sensitive parameter. Tumor immune status is also likely to vary across patients. The nominal E:T ratio in the model is assumed to be low (1:150) in a solid tumor (21,24). To investigate the potential impact of tumor T cell number on tumor trimer concentrations, a sensitivity analysis was performed varying E:T ratio from 1:15 to 1:1500 and assuming a constant number of tumor cells. Predicted tumor trimer concentration correlates with E:T ratio (Fig. 5b), suggesting T cells in the tumor is a sensitive parameter.

Fig. 4

Model simulated a serum PK and b tumor trimer concentrations following IV infusion of PF-06671008 at 0.01, 0.1, and 1 μg/kg QW to cancer patients

Fig. 5

Model simulated tumor trimer concentrations at a different P-cadherin receptor expression values (1000–28,706 receptors/cell) and b different E:T ratios (1:1500–1:15) following IV infusion of PF-06671008 at 0.1 μg/kg QW to cancer patients

Discussion

Complex Exposure-Response Relationships for CD3 Bispecific Molecules

Bispecific antibodies are emerging as a leading class of biotherapeutic drugs in oncology, with the potential to enhance efficacy, increase tumor targeting and reduce systemic toxicity compared to their monospecific counterparts. These formats can vary in their molecular weight, PK, and ability to support immune effector functions. Perhaps more significantly, they can also vary in geometry a number of antigen binding sites, and the intrinsic affinity of individual arms (29). As a result of this complexity, dose-response relationships for bispecific antibodies can be non-intuitive and difficult to rationalize.

An additional complexity emerges for the CD3 bispecific T cell retargeting modality, where efficacy is driven by the formation of a trimer between the drug, T cell, and tumor cell. A bell-shaped concentration-response relationship can be observed (Fig. 1b), which is a well-described phenomenon for ternary complexes (30–33). When concentrations of antibodies are low, conditions favor the formation of trimers, with an optimal antibody concentration needed for trimer formation. However, as concentrations increase further, antibodies will be in excess and the equilibrium will shift to the formation of dimers between antibodies and T cell, or antibodies and tumor cell. This results in a decrease of response as dimers cannot trigger cytotoxicity. Since trimer concentration is a function of drug Kd values, tumor antigen expression, CD3 expression, and E:T ratio, a single drug concentration could potentially result in different trimer concentrations. Therefore, interpretation of response by drug exposure alone can be misleading. For the CD3 bispecific molecule discussed in this manuscript (PF-06671008), a bell-shaped dose-response relationship was not observed in mouse xenograft studies. This is probably because there was high P-cadherin expression on the tumor cell lines studied and good infiltration of T cells. In addition, PF-06671008 is a potent drug with low-Kd values for P-cadherin and CD3. As a result sufficient, trimer concentrations were achieved at each dose for efficacy. The bell-shaped relationship has been confirmed for other CD3 bispecific molecules in-house, where target expression is lower and/or affinity is weaker. It has also been observed in the literature from modeling of in vitro and in vivo experimental data (34,35).

Translational QSP Model for CD3 Bispecific Molecules

QSP models, which map out the causal path between drug administration and effect in a mechanistic framework, can be a useful tool to deconvolute complex mechanisms (36). Some examples of the use of mechanistic PK/PD models to quantify and understand the system dynamics of CD3 bispecific molecules are emerging in the literature. For example, Jiang et al. (34) proposed a cell-killing model based on target cell-biologic-effector cell complex formation and used it to describe and predict in vitro cytotoxicity data for multiple T cell redirecting bispecific antibodies under different experimental conditions. Campagne et al. (37) developed a PK/PD model for a bispecific CD123/CD3 DART molecule in non-human primates. The model describes DART molecule binding to peripheral CD3 expressing cells and CD123+ cells, T cell trafficking, activation, and expansion, and resulting peripheral depletion of CD123 cells.

In this manuscript, a translational QSP model is proposed for CD3 bispecific T cell retargeting molecules, capable of predicting trimer formation and linking it to tumor cell killing in in vivo efficacy models. In addition, the mechanistic nature of the model enables the integration of patient data/parameters and subsequent clinical predictions. The model consists of 3 parts describing the central, tumor, and effect compartments (Fig. 1a). The first part includes the bispecific antibody PK, binding to circulating T cells, and binding to the soluble target (when applicable) in the central compartment. The second part describes the distribution of the antibodies to the tumor compartment using mechanistic tumor penetration equations, and parameters calculated based on the drug’s molecular weight and tumor size (19,26,38). If the model is being used for a liquid tumor, these drug exchange tumor penetration parameters can simply be removed, as liquid tumors are assumed to provide less of a diffusion barrier than solid tumors, and equilibrium can be assumed between drug concentration in the central compartment and tumor interstitium. In the tumor compartment, the model incorporates binding of the drug to CD3 on T cells and the specific antigen on tumor cells to form inactive dimers and ultimately the active trimers. In the mouse model, a simple description of T cell expansion and contraction is included, constructed using mouse TIL data and published information. For translation of this model to human, data on T cell kinetics was not available and instead a baseline concentration of T cells was assumed with no proliferation.

In the third part of the model, the trimer concentration is used as the basis for quantifying tumor volume reduction using a tumor growth inhibition model. The model used is a transduction model describing tumor cell growth and tumor cell killing (as a function of the tumor trimer concentration). The model parameters from each mouse study can be used to calculate a secondary parameter called the TSC. This is the concentration of trimer where the tumor is neither growing nor regressing and can be considered as the minimum concentration of tumor trimer required for efficacy. The TSC is a useful parameter which can be used as a pharmacodynamic index to rank compounds, or to understand the difference in compound potency across mouse xenograft models, or as the denominator in therapeutic index calculations.

Application of the QSP Model to Quantify PK/PD Relationship for PF-06671008 in Mouse Xenograft Models

The model was used to quantify the preclinical PK/PD relationship of a CD3 bispecific molecule targeting P-cadherin (PF-06671008). To implement the model, the first step was to collect drug and system parameters describing the mechanism of action in the mouse. To calculate trimer concentration in the tumor, receptor expression of P-cadherin was determined for the HCT-116 and SUM-149 human tumor cell lines used in the mouse xenograft experiments. P-cadherin receptor expression in both cell lines (28,706 for HCT-116 and 17,500 for SUM-149) was lower than the expression of CD3 on T cells (100,000 (17,18)). This is typical for CD3 bispecific molecules as an expression of most tumor targets is less than 100,000 and as a result, tumor antigen receptor expression can be limiting and a key driver of efficacy. This was exemplified for a carcinoembryonic antigen T cell bispecific (CEA-TCB) for the treatment of solid tumors. CEA-TCB activity was found to be strongly correlated with CEA expression, with a higher potency observed in highly CEA-expressing tumor cells, with a threshold of 10,000 CEA-binding sites/cell (39). Target affinity data was also required to calculate trimer concentration. PF-06671008 binds to P-cadherin with a Kd of 0.47 nM and CD3 with a Kd of 11.4 nM (7). Binding to the tumor target antigen is often more potent than binding to CD3 on T cells in order to target the CD3 bispecific toward the tumor and away from peripheral tissues (35). In addition, strong binding to CD3 has been shown to drive more rapid clearance of an anti-CD3/anti-CLL1 bispecific in preclinical in vivo models (40).

The QSP model was used to integrate the mouse PK for PF-06671008 with the TGI data and to calculate TSCs in T cell engrafted (HCT-116) and T cell adoptive transfer (HCT-116 and SUM-149) established mouse tumor models. TSC values were very similar in the adoptive transfer model for both the SUM-149 and HCT-116 tumor cell lines (0.0092 and 0.011 pM respectively, with overlapping 80% confidence intervals). In contrast, a sixfold higher TSC value was obtained in the T cell engrafted versus T cell adoptive transfer model with the same cell line (HCT116, 0.064pM), and the respective 80% confidence intervals do not overlap. This is probably due to differences in T cell engraftment between the two mouse tumor models. In the T cell engrafted model, the T cells are administered as freshly isolated human PBMCs, 7 days prior to drug administration. In contrast, in the adoptive transfer model, activated T cells are given 1 day post-drug treatment. There are also other factors which can result in different TSCs including initial tumor size and differences in tumor growth rates.

Translation of the Model to the Clinic

The first step in translation to human was the prediction of the clinical PK parameters. For PF-06671008, circulating soluble target can act as a sink for the drug and reduce free drug exposure by forming complexes with PF-06671008. The reduction of free sPcad concentrations in cynomolgus monkey following dosing of PF-06671008 has been reported previously (16). The human PK of PF-06671008 was predicted from cynomolgus monkey PK using a two-compartmental PK model which incorporates binding to sPcad. Levels of sPcad were measured in healthy volunteers and in breast, colon, and lung cancer patients and the median concentration in cancer patients was used in the human model.

The next step in the clinical translation process was the incorporation of human systems parameters into the QSP model. These parameters are summarized in Table III and include T cell concentration in the circulation and tumors, tumor cell concentration, and typical tumor volumes in cancer patients. Values for all of these parameters were obtained from the literature. CD3 receptor expression was kept the same as the mouse model (which used human T cells or PBMCs). P-cadherin expression of 28,706 was used in the clinical simulations. This was the value from the HCT-116 cell line and represents a medium-high level expression of P-cadherin measured across human tumor cell lines used in in vitro cytotoxicity experiments (874–37,582 (15)).

The model simulated serum PK and tumor trimer concentrations following IV infusion of PF-06671008 at 0.01, 0.1, and 1 μg/kg QW to cancer patients are shown in Fig. 4 a and b, respectively. In human, the terminal half-life of PF-06671008 was predicted to be approximately 1 day. The concentration of trimer in the tumor, which is the more relevant concentration for efficacy, accumulates slowly (Cmax approx. 2 days post first dose) and persists for longer (Fig. 4b). This is due to slow diffusion of the drug into the tumor and formation of a more stable trimer which is retained within the TME. Since receptor expression of tumor target was known to be a key parameter, a sensitivity analysis was completed using the human model with P-cadherin expression varying from 1000 to 28,706 receptors/cell. This analysis confirmed that P-cadherin receptor expression was a sensitive parameter and that concentration of trimer formed in the tumor correlates with expression level (Fig. 5a). This has an impact on predicted clinical efficacy with a higher dose required for efficacy in patients with lower P-cadherin expression. In addition, the T cell number in the tumor was found to be a sensitive parameter (Fig. 5b), with a higher predicted concentration of trimer in the tumor with increasing E:T ratio. High doses of PF-06671008 were also simulated, to check to see where the bell-shaped relationship might be observed. At doses of > 1.8 mg/kg, a reduction in tumor trimer concentration is predicted with increasing dose levels (Supplementary Fig. 4). However, at these doses, the predicted trimer concentrations in the tumor are high enough that good responses would be expected (assuming the doses would be tolerated). A translational flow diagram describing the steps taken to translate CD3 bispecific drugs from preclinical TGI data in the mouse to human is shown in Fig. 6.

Fig. 6

QSP model-based strategy for translating preclinical data for CD3 bispecific compounds to the clinic. “Biomeasures” can be defined as system-dependent parameters. TAA is a tumor-associated antigen

The translational QSP model described for CD3 bispecific compounds can be used to drive decision-making at different stages of the drug discovery and development continuum. At early stages, the model can be used to provide guidance on the compound selection, by predicting optimal Kd values for CD3 and the tumor antigen. This can be achieved by modeling of in vitro data, using a reduced version of the model without the PK (central and peripheral) compartments. For example, the model was previously used to describe the in vitro exposure-response of PF-06671008 in cytotoxicity assays and was able to simultaneously describe the kinetics of tumor and T cells at various E:T ratios (16). Once a lead compound has been selected, the model can be used to predict clinical doses and regimens and to optimize efficient clinical study design (41). A precision medicine approach could be adopted, whereby parameters in the model such as immune cell numbers, or tumor target expression levels, are tailored to individual characteristics of patients. This could result in a recommendation of different doses for different patients. The model has also been used to predict clinical starting dose for PF-06671008 using a minimal biological effect level approach (MABEL), which is recommended for CD3 bispecific constructs due to their immune agonistic activity following target engagement (16,42). A recent analysis by the FDA concluded that receptor occupancy-based methods were not advised for CD3 bispecifics. The QSP modeling approach is more suitable to determine MABEL as efficacy is driven by drug bound to both T cells and tumor cell, rather receptor occupancy of either target singly. It is also independent of E:T ratio or other experimental specificities.

The model in its current state is very useful for a range of tasks from optimization of drug design to clinical dose predictions. However, opportunities exist to improve the model. For example, the current model includes an empirical description of T cell activation/proliferation in mouse, constructed based on TIL analysis across dose and time. A more mechanistic model could be developed by collection and characterization of more tumor lymphocyte kinetic data across species. In addition, the model is based upon a “well-mixed” hypothesis in which tumor target and T cells are assumed to be homogeneously distributed throughout the tumor environment with equal opportunity for trimer formation. However, tumors are known to be a complex environment with a heterogeneous distribution of T cells and tumor cells expressing target. Future versions of the model will take this into account.

Conclusion

The mechanistic PK/PD model and translational framework described for CD3 bispecific molecules provide a holistic solution for quantitative decision-making throughout the drug discovery and development process. In this manuscript, the use of the model to characterize the in vivo PK/PD relationship of a P-cadherin/CD3 bispecific construct (PF-06671008) across mouse efficacy models is described. The model can also be translated to the clinic for human PK/PD predictions and sensitivity analysis to determine important parameters driving efficacy. The model can be applied at early stages to aid in the design of CD3 bispecific constructs and to select molecules with optimal properties.

(DOCX 362 kb)