Linear pharmacokinetic parameters for monoclonal antibodies are similar within a species and across different pharmacological targets: A comparison between human, cynomolgus monkey and hFcRn Tg32 transgenic mouse using a population-modeling approach.

The linear pharmacokinetics (PK) of therapeutic monoclonal antibodies (mAbs) can be considered a class property with values that are similar to endogenous IgG. Knowledge of these parameters across species could be used to avoid unnecessary in vivo PK studies and to enable early PK predictions and pharmacokinetic/pharmacodynamic (PK/PD) simulations. In this work, population-pharmacokinetic (popPK) modeling was used to determine a single set of 'typical' popPK parameters describing the linear PK of mAbs in human, cynomolgus monkey and transgenic mice expressing the human neonatal Fc receptor (hFcRn Tg32), using a rich dataset of 27 mAbs. Non-linear PK was excluded from the datasets and a 2-compartment model was applied to describe mAb disposition. Typical human popPK estimates compared well with data from comparator mAbs with linear PK in the clinic. Outliers with higher than typical clearance were found to have non-specific interactions in an affinity-capture self-interaction nanoparticle spectroscopy assay, offering a potential tool to screen out these mAbs at an early stage. Translational strategies were investigated for prediction of human linear PK of mAbs, including use of typical human popPK parameters and allometric exponents from cynomolgus monkey and Tg32 mouse. Each method gave good prediction of human PK with parameters predicted within 2-fold. These strategies offer alternative options to the use of cynomolgus monkeys for human PK predictions of linear mAbs, based on in silico methods (typical human popPK parameters) or using a rodent species (Tg32 mouse), and call into question the value of completing extensive in vivo preclinical PK to inform linear mAb PK.

Introduction

Therapeutic monoclonal antibodies (mAbs) have become an important class of drugs, with over 40 antibody-based therapies approved by the US FDA across multiple indications and many more in clinical trials. Advances in antibody engineering have enabled rapid progress from the first generation of highly immunogenic murine and chimeric antibodies to better tolerated humanized and fully human mAbs. Recently, the variety of antibody-like modalities has evolved further to include Fc-fusion proteins, antibody-drug conjugates and bispecific antibody products. In addition to their exquisite specificity and potency, mAbs are successful therapeutics due to their long pharmacokinetic (PK) half-life. Low clearance of mAbs from the systemic circulation enables them to be administered less frequently than their peptide or small molecule counterparts, which is more convenient for the treatment of chronic diseases.

The PK properties of mAbs are a function of their large size (150 kDa), relative polarity, Fc-receptor binding and specific binding to target antigens. The primary elimination route for mAbs is cellular uptake followed by proteolytic degradation. There are two distinct catabolic pathways for mAbs. The first is a non-specific, linear (first-order) clearance (CL) pathway mediated by fluid phase pinocytosis or unspecific fluid phase endocytosis. This common pathway shared by endogenous IgG and therapeutic mAbs operates independently of the specific interaction between a mAb and its pharmacological target. This pathway is not easily saturated at therapeutic doses, and tends to result in linear CL. FcRn functions as a salvage receptor to protect IgG from rapid intracellular catabolism, and is responsible for the long half-life of endogenous IgG and exogenous IgG-based therapeutic proteins.

The second catabolic pathway is a non-linear (target mediated) CL pathway mediated by the specific interaction between the Fab region of the antibody and its pharmacological target. This pathway is often referred to as target-mediated drug disposition (TMDD). Following binding of the mAb to its target on the cell surface, the mAb-antigen complex is internalized and then trafficked via the endosomes to the lysosomes where the complex is degraded. When the target binding is saturated, the relative importance of target binding to overall disposition is diminished and the mAb is eliminated by first order process. For mAbs exhibiting this pathway, disposition depends upon the concentration and distribution of the mAb, along with target receptor expression, internalization and turnover rates. Certain mAbs to soluble targets can also undergo TMDD driven by binding of 2 or more mAbs to form multimeric complexes that are rapidly eliminated by phagocytosis. mAbs cleared primarily by TMDD will have dose-dependent non-linear elimination. For these mAbs, PK is distinctly different from the catabolism of endogenous IgG, with higher CL and shorter half-life values at lower doses.

The rate and extent of mAb distribution is very slow and depends upon extravasation in tissue, distribution within the particular tissue and degradation. The convective CL of mAbs from the tissue is thought to be more efficient than the process of convective extravasation, thereby maintaining relatively low mAb concentrations in the interstitial fluid. , As a result of this, mAbs often have small apparent volumes of distribution. ,

The ability to predict the PK of a drug prior to first-in-human studies is of utmost importance to reduce attrition in Phase 1. Preclinical testing of mAbs often occurs in a rodent species and non-human primate (NHP) to understand efficacy and toxicity prior to human dosing. Cynomolgus monkey is the preferred strain of NHP for preclinical PK and toxicological studies due to a high genetic similarity with humans, and therefore greater likelihood of target antigen sequence homology, comparable binding affinities for cynomolgus monkey vs human FcRn, and similar tissue cross-reactivity profiles. Cynomolgus monkey is also the preferred species for predicting the PK of mAbs with linear CL in humans: several groups have reported the successful use of fixed allometric exponents to predict CL and volume of distribution of mAbs in human from data in cynomolgus monkey.[11-14] For mAbs that exhibit non-linear CL due to TMDD, scaling of PK is more challenging. In order to take into account the kinetics of mAb binding to its target, a mechanistic TMDD model is required, with proper exploration of species differences in target expression and binding. ,[15-17]

In this study, the objective was to complete a comprehensive analysis of the linear PK of therapeutic mAbs using popPK methods. The dataset was composed of 27 Pfizer mAbs where PK had been generated in at least one of three species: human, cynomolgus monkeys and human FcRn Tg32 homozygous transgenic mice. Given that the linear, Fc-mediated elimination is a common pathway shared by both endogenous IgG and therapeutic IgG of mAbs, it was theorized that a single set of ‘typical’ linear PK parameters could be estimated for each species, describing the linear PK of all mAbs in the dataset. PopPK is an ideal technique to investigate this because it can separately estimate variability, including variability between mAbs, between individuals in a given mAb dataset and also random, non-specific error. The resulting ‘typical’ linear PK parameter estimates could be used for designing PK/PD and toxicology studies, dose predictions and viability assessments. For tractable mAbs with linear PK in preclinical species, it could forfeit the requirement for allometric scaling for clinical PK predictions; instead the typical PK parameters could be used as a substitute until clinical PK is obtained. In addition to cynomolgus monkey, we also studied the potential of a human FcRn Tg32 transgenic mouse model to predict linear PK of mAbs in human. The availability of a rodent model to accurately estimate human PK of mAbs would enable earlier predictions, before cynomolgus monkey data is routinely available. In addition to the single species analysis, a combined analysis was performed on the entire dataset of human, cynomolgus monkey and hFcRn Tg32 transgenic mice and used to estimate allometric co-efficients between species. While our primary focus was to explore PK predictions utilizing a variety of species, we also investigated the value of an in vitro assay measuring self-association in predicting the CL of mAbs in the dataset. The analysis herein provides robust strategies for predicting linear human PK of mAbs, which could improve throughput for lead drug candidate selection, and potentially increase the overall success while decreasing the time for non-clinical development of mAbs.

Results

mAb PK dataset and selection of linear dose range

Properties of the mAbs included in this study are summarized in Table 1. Of the 27 mAbs analyzed, 12 were IgG1 and 15 were IgG2; 16 were specific for soluble ligand targets, 9 bound to membrane targets and 2 had both membrane and soluble targets; 16 of the mAbs were fully human, 10 were humanized and 1 was from a human phage display library. All of the mAbs had similar binding Kd values to FcRn and had wild type sequences for the FcRn binding region. For 18 of these mAbs, clinical PK data were available; cynomolgus monkey PK data were available for 23 mAbs; and 11 mAbs had PK in hFcRn Tg32 transgenic mouse. As discussed in the materials and methods section, non-linear data were removed from the dataset to enable analysis of linear PK only. This was completed via a combination of visual inspection and application of an algorithm applied to a linear regression of dose and AUC to test for deviation of the slope from 1 (Fig. 1). The remaining linear dose range and number of dose levels included in the popPK analysis for each mAb are summarized in Table 1.

Table 1.

Monoclonal antibodies in this study.

mAb number	Immunoglobin type	Human or Humanized	Type of mAb target	AC-SINS Score	Species	Linear dose range (mg/kg)	Dose levels excluded (mg/kg)	Number of dose levels included
mAb1	IgG2	Human	membrane bound	—	Cyno	2	None	1
mAb2	IgG2	Humanized	soluble	—	Cyno	10– 100	None	3
mAb3	IgG2	Humanized	soluble	3	Cyno, Tg32 Mice	1–10, 5	None, None	2, 1
mAb4	IgG2	Human	membrane bound	—	Cyno	1–10	0.01–0.3	2
mAb5	IgG2	Human	membrane bound	—	Human, Cyno	1–15, 5–50	0.5, None	8, 2
mAb6	IgG2	Human	membrane bound	—	Human	0.1–15	None	6
mAb7	IgG2	Human	membrane bound	—	Human, Cyno	0.2–20, 30–100	None, 3–5	8, 2
mAb8	IgG1	Human	membrane bound	—	Cyno	3–30	0.3	2
mAb9	IgG1	Humanized	soluble	—	Human, Cyno	1.5–5, 1–10	0.5, None	2, 2
mAb10	IgG2	Humanized	soluble	—	Human, Cyno	1–24, 1–100	0.3, 0.1	6, 3
mAb11	IgG2	Human	soluble	—	Cyno	1– 5	None	2
mAb12	IgG1	Human	soluble	—	Human, Cyno	5–600, 2	None	6, 1
mAb13	IgG1	Synthetic	soluble	—	Cyno	2	None	1
mAb14	IgG2	Humanized	soluble	—	Human, Cyno	1–18, 3–100	0.3, 0.1–0.3	6, 3
mAb15	IgG2	Humanized	soluble	—	Human	1–6	None	3
mAb16	IgG2	Humanized	soluble	—	Human, Cyno	0.1–10, 0.8–16.5	0.01–0.03, None	4, 3
mAb17	IgG2	Humanized	soluble	—	Human, Cyno	1–10, 10–200	0.1–0.3, None	4, 3
mAb18	IgG1	Human	membrane bound	24	Human, Cyno, Tg32 Mice	30–120 mg, 2–10, 5	None, None, None	3, 2, 1
mAb19	IgG1	Human	membrane bound	17	Human, Cyno, Tg32 Mice	1–8, 0.5–5, 5	0.5, None, None	4, 2, 1
mAb20	IgG2	Human	membrane bound and soluble	0	Human, Cyno, Tg32 Mice	0.3–10, 1–100, 5	0.03– 0.1, None, None	3, 3, 1
mAb21	IgG1	Human	membrane bound and soluble	6	Human, Cyno, Tg32 Mice	10–800 mg, 10–100, 5	1– 3 mg, None, None	4, 2, 1
mAb22	IgG2	Human	soluble	1	Human, Cyno, Tg32 Mice	22–700 mg, 0.5–5	7 mg, None, None	6, 2
mAb23	IgG1	Humanized	soluble	0	Human, Cyno, Tg32 Mice	3, 1–100, 5	None, None, None	1, 2, 1
mAb24	IgG1	Humanized	soluble	1	Human, Cyno, Tg32 Mice	3, 1–100, 5	None, None, None	1, 2, 1
mAb25	IgG1	Human	soluble	10	Human, Cyno, Tg32 Mice	0.3– 11.25, 50, 5	None, None, None	5, 1, 1
mAb26	IgG1	Human	membrane bound	21	Tg32 Mice	5	None	1
mAb27	IgG1	Human	soluble	2	Tg32 Mice	5	None	1

Figure 1.

(a) 2-compartment PK model and (b) algorithm to test for linearity of PK data.

Population PK analysis across species

The concentration versus time relationships in Tg32 mouse, cynomolgus monkey and human mAb datasets were described using a 2-compartmental PK model (Fig. 1) with inter-individual variability (IIV) on CL and volume of the central compartment (V1). Residual error was determined per compound and covariance was estimated between CL and V1. This model adequately captured the PK of all compounds in each species. The PK data for the mAbs in each dataset could be described using a single set of popPK parameters. The parameter estimates for each species are shown in Table 2, with 95% confidence intervals. The popPK estimate of human CL (0.15 mL/h/kg (0.14–0.16)) was in the same range as the CL of endogenous IgG (0.125 mL/h/kg). , CL was lowest in humans (0.15 mL/h/kg (0.14–0.16)), followed by cynomolgus monkey (0.27 mL/h/kg (0.24–0.30)) and then Tg32 mouse (0.35 mL/h/kg (0.28– 0.41)). Population estimates of the volumes of distribution in the central (V1) and peripheral (V2) compartments were typically small and approximated plasma volume. The inter-compartmental clearance parameter, Q, varied the most across species and was estimated to be 0.27 mL/h/kg (0.25– 0.30) in human, 1.00 mL/h/kg (0.8–1.20) in cynomolgus monkey and 4.40 mL/h/kg (3.17– 5.62) in Tg32 mouse. Inter-individual variability in V1 was low and moderate-high for CL, which may be expected from the methodology chosen in which PK parameters were estimated for all compounds combined, instead of estimating separate parameters for each compound. Covariance between CL and V1 was low, representing low correlation between the random effects on the parameter estimates. The relative standard error (RSE) was lower than 35% for all parameters representing low uncertainty in parameter estimation.

Table 2.

Single species popPK parameter estimates for mAbs with linear CL.

		Human (n=18)		Cynomolgus Monkey (n=23)		Tg32 hFcRn Transgenic Mouse (n=11)
Parameter	Unit	Value (95%-CI)	RSE (%)	Value (95%-CI)	RSE (%)	Value (95%-CI)	RSE (%)
CL	mL/h/kg	0.15 (0.14–0.16)	3.03	0.27 (0.24–0.30)	5.22	0.35 (0.28–0.41)	9.06
V1	mL/kg	46.31 (45.14–47.48)	1.29	39.29 (37.16–41.41)	2.76	59.28 (54.65–63.90)	3.98
Q	mL/h/kg	0.27 (0.25–0.30)	5.12	1.00 (0.80–1.20)	10.33	4.40 (3.17–5.62)	14.24
V2	mL/kg	31.47 (28.63–34.31)	4.60	27.56 (24.83–30.29)	5.05	60.54 (52.80–68.29)	6.53
IIV CL	—	0.48 (0.43–0.53)	5.29	0.38 (0.31–0.46)	9.48	0.41 (0.29–0.53)	15.34
COV CL-V1	—	0.09 (0.07–0.10)	11.48	0.09 (0.05–0.13)	23.28	0.11 (0.04–0.18)	33.55
IIV V1	—	0.09 (0.08–0.11)	6.66	0.10 (0.07–0.13)	14.90	0.12 (0.07–0.18)	23.91

CL: clearance from the central compartment, V1: volume of the central compartment, Q: inter-compartment distribution clearance, V2: volume of the peripheral compartment. RSE(%): relative standard error, calculated as: standard error of estimation / estimated value x 100%. 95% CI: 95% confidence interval calculated from the standard error. IIV: inter-individual variability. COV CL-V1: covariance between clearance and volume. Residual errors per compound are shown in Supplementary Table 1.

Results from a jackknife analysis showed consistent estimation of PK parameters without significant influence from removal of one mAb from the dataset at a time. PK parameters from the jackknife analysis with 95% confidence intervals are shown in Supplementary Fig. 1. Residual error per compound and diagnostic plots of observed concentration versus individual and population predictions in each species are included in Supplementary Tables 1 and 2 and Supplementary Figs. 2 and 3.

Distribution of clearance and volume of distribution of the central compartment

The distribution of CL estimates and the individual estimates of CL (with IIV) for each mAb across species are shown in Fig. 2. This plot shows a typical log normal distribution of individual CL estimates. In general, CL is lowest in human, followed by cynomolgus monkey and then Tg32 mouse. The distribution of population and individual estimates of V1 for each mAb across species are shown in Fig. 3. The estimates of V1 are normally distributed in cynomolgus monkey and human and center on plasma volume (30– 50 mL/kg). It should be noted that more variability was observed in individual estimates of V1 in Tg32 mouse compared with cynomolgus monkey or human. This may have been due to restricted time points in the initial phase of PK sample collection in Tg32 mouse.

Figure 2.

Distribution of population and individual mAb estimates (with variability) of clearance (CL) in the combined human, cynomolgus monkey and hFcRn Tg32 mouse dataset.

Figure 3.

Distribution of population and individual mAb estimates (with variability) of volume of distribution of the central compartment (V1) in the combined human, cynomolgus monkey and hFcRn Tg32 mouse dataset.

Comparison of human population PK parameter estimates with population PK of literature mAbs

In order to put the results into context, popPK parameters determined for the current mAb dataset were compared with therapeutic mAbs with linear CL in clinical studies. 5 fully human, humanized or chimeric therapeutic mAbs were identified with linear CL in humans and 2-compartment popPK analyses published in the literature. Population estimates of CL, Q, V1 and V2 for the Pfizer mAbs and bevacizumab, infliximab (2 different indications), pertuzumab, rituximab and trastuzumab in human are shown in Table 3. The parameter estimates reported for the comparator clinical mAbs are very close to the population parameter estimates generated from the analysis of 18 clinical mAbs described herein. An outlier was infliximab for ulcerative colitis (UC), which had a reported mean CL of 0.24 mL/h/kg and Q of 4.25 mL/h /kg, which deviates significantly from the population CL estimates of 0.15 mL/h /kg (0.14– 0.16) and Q estimates of 0.27 mL/h /kg (0.25– 0.3), estimated in this analysis. Fig. 4 shows the concentration versus time profiles for the median and 95% prediction interval with observed dose-normalized concentration-time data for the 18 clinical mAbs in this study and the 5 clinical therapeutic mAbs (bevacizumab, infliximab, pertuzumab, rituximab and trastuzumab).

Table 3.

Comparison of typical popPK parameter estimates in human with test set of clinical mAbs with linear PK.

		Human	Literature clinical mAbs
Parameter	Unit	Value (95% CI)	Bevacizumab	Infliximab AS	Infliximab UC	Pertuzumab	Rituximab	Trastuzumab
CL	mL/h/kg	0.15 (0.14 – 0.16)	0.12	0.16	0.24	0.13	0.15	0.13
V1	mL/kg	46.31 (45.14 – 47.48)	38.0	43.7	47.0	39.1	42.6	42.1
Q	mL/h/kg	0.27 (0.25 – 0.3)	0.35	1.02	4.25	0.33	0.39	0.29
V2	mL/kg	31.47 (28.63 – 34.31)	39.4	42.0	59.0	30.9	52.0	68.4

CL: clearance from the central compartment, V1: volume of the central compartment, Q: inter-compartment distribution clearance, V2: volume of the peripheral compartment. 95% CI: 95% confidence interval calculated from the standard error. AS: ankylosing spondylitis; UC: ulcerative colitis

Figure 4.

Median, 5th and 95th percentiles of 200 bootstrap samples with the combined species PK model using the observed dose-normalized concentration (ng/mL) vs. time (hours) data for all of the mAbs in this study. The red lines indicate the concentration vs. time profile of the 5 clinical therapeutic mAbs (bevacisumab, infliximab (for both ankylosing spondylitis and ulcerative colitis), pertuzumab, rituximab and trastuzumab).

Clearance outliers and correlation with non-specific interactions

Four of the mAbs in the human dataset had CL values 2-fold higher than the population estimate of 0.15 mL/h/kg (0.14–0.16) (Table 5). In order to investigate possible non-specific (off target) binding properties, a subset of the mAbs in the dataset (n=11) were studied using an in vitro affinity capture self-interaction nanoparticle spectroscopy (AC-SINS) assay. This assay assesses self-association, which is often coupled with other poor physicochemical characteristics, including non-specific interactions. AC-SINS uses gold nanoparticles pre-coated with anti-human Fc polyclonal antibodies to capture test mAbs. Self-interactions of immobilized mAbs lead to clustering of the gold nanoparticles, which is measured by a shift in absorbance due to changes in their optical properties. The 11 mAbs studied had AC-SINS scores ranging from 0–24 (Table 1). A plot of CL (mL/h/kg) vs. AC-SINS score for each species is shown in Fig. 5. An increase in CL was observed with an increase in AC-SINS scores in all species. This data suggests that non-specific interactions leading to off-target binding may result in faster clearance than predicted by the popPK estimates.

Comparison of observed human PK parameters for individual mAbs with ‘typical’ human popPK parameter estimates.

	CL mL/h/kg	V1 mL/kg	Q mL/h/kg	V2 mL/kg	%RMSE
Human PopPK	0.15	46.31	0.27	31.47	—
Observed Human PK
mAb5	0.26	54.50	0.30	33.29	50.6
mAb6	0.12	52.58	0.32	31.76	23.1
mAb7	0.21	57.04	0.68	16.38	31.4
mAb9	0.12	46.58	0.37	40.47	24.7
mAb10	0.07	34.62	0.25	24.39	53.8
mAb12	0.17	40.49	0.25	32.75	18.2
mAb14	0.18	38.51	0.31	15.70	32.8
mAb15	0.17	39.11	0.43	44.99	19.7
mAb16	0.10	40.23	0.36	42.42	41.6
mAb17	0.09	48.30	0.19	64.89	54.2
mAb18	0.46	65.98	0.20	167.77	361.6
mAb19	0.46	52.18	0.19	6.41	121.9
mAb20	0.14	39.54	0.39	20.48	33.5
mAb21	0.10	40.07	0.80	14.00	54.8
mAb22	0.05	41.67	0.29	42.73	111.8
mAb23	0.11	54.96	0.33	49.29	39.7
mAb24	0.07	31.37	0.21	22.75	78.4
mAb25	0.32	55.76	0.32	27.73	104.1

CL: clearance from the central compartment, V1: volume of the central compartment, Q: inter-compartment distribution clearance, V2: volume of the peripheral compartment.

%RMSE: percent root mean square error.

Figure 5.

Clearance vs. AC-SINS score for a subset of 11 mAbs in the dataset in human, cynomolgus monkey and hFcRn Tg32 mouse.

Combined analysis and estimation of allometric exponents

The Tg32 mouse and human datasets, cynomolgus monkey and human datasets and datasets from all three species were combined. 2-compartment human PK parameters were estimated in the combined datasets and allometric exponents to scale the PK parameters from preclinical species to human (Table 4). The human popPK estimates of CL, Q, V1 and V2 in the combined datasets were similar to those estimated previously (Table 2). The estimated allometric exponent for scaling Tg32 mouse CL to human was 0.9 (0.88– 0.92), cynomolgus monkey CL to human was 0.81 (0.77– 0.85) and all preclinical data CL to human was 0.89 (0.87– 0.91). Volumes of distribution from central and peripheral compartments in general scaled with an allometric exponent of approximately 1. Population estimates of the allometric exponent for Q were between 0.57 and 0.67 for Tg32 mouse, cynomolgus monkey and all species.

Table 4.

Combined species PK parameter estimates for mAbs with linear CL.

Combined Datasets		Tg32 Mouse, Cyno and Human (n=27 mAbs)		Cyno and Human (n=23 mAbs)		Tg32 Mouse and Human (n=23 mAbs)
Parameter	Unit	Value (95%-CI)	RSE (%)	Value (95%-CI)	RSE (%)	Value (95%-CI)	RSE (%)
CL	mL/h/kg	0.16 (0.15–0.16)	2.86	0.15 (0.14–0.16)	2.90	0.15 (0.14–0.16)	2.97
V1	mL/kg	45.19 (44.08–46.31)	1.26	45.89 (44.75–47.03)	1.27	46.41 (45.26–47.56)	1.26
Q	mL/h/kg	0.28 (0.25–0.31)	5.06	0.29 (0.26–0.32)	5.08	0.28 (0.25–0.31)	4.97
V2	mL/kg	30.81 (28.15–33.46)	4.40	31.14 (28.61–33.68)	4.15	32.17 (29.31–35.04)	4.55
CL: α	—	0.89 (0.87–0.91)	1.16	0.81 (0.77–0.85)	2.22	0.90 (0.88–0.92)	1.24
V1: β	—	0.98 (0.97–0.99)	0.54	1.04 (1.02–1.06)	0.87	0.97 (0.96–0.98)	0.55
Q: γ	—	0.67 (0.63–0.71)	2.87	0.57 (0.48–0.67)	8.36	0.67 (0.64–0.70)	2.19
V2: δ	—	0.95 (0.93–0.98)	1.28	1.07 (1.03–1.11)	2.02	0.93 (0.91–0.94)	0.84
IIV CL	—	0.47 (0.43–0.52)	4.61	0.45 (0.41–0.50)	4.64	0.47 (0.42–0.52)	5.00
COV CL-V1	—	0.08 (0.06–0.10)	11.25	0.08 (0.06–0.10)	10.82	0.08 (0.07–0.10)	10.83
IIV V1	—	0.11 (0.09–0.12)	6.26	0.10 (0.09–0.11)	6.44	0.10 (0.08–0.11)	6.31

CL: clearance from the central compartment, V1: volume of the central compartment, Q: inter-compartment distribution clearance, V2: volume of the peripheral compartment. %RSE: relative standard error, calculated as: standard error of estimation / estimated value x 100% 95% CI: 95% confidence interval calculated from standard error. α, β, γ, and δ: allometric exponents estimated for CL, V1, Q and V2 respectively, to scale from the preclinical species to human using the equation:

IIV: inter individual variability. COV CL-V1: covariance between clearance and volume. Residual errors per compound are shown in Supplementary Table 2. For each combined dataset the total ’n’ includes n=18 mAbs with clinical data, n=23 mAbs with cyno data and n=11 mAbs with Tg32 mouse data.

Linear mAb human PK prediction strategies

Different methods were applied to predict human PK of the mAbs in this dataset. Note, a limitation to this analysis is that there was no separate test dataset to validate the inter-species scaling.

Use of ‘typical’ PK parameters for human

Simulations of the popPK estimates of CL, Q, V1 and V2 estimated from the human dataset, compared with observed clinical PK profiles for the individual mAbs, are shown in Fig. 6. Root mean square errors (RMSEs) between observed and predicted data are shown in Table 5. Only 4 of 18 mAbs have RMSEs of >100%, indicating that the human popPK parameters can adequately predict PK for the majority of mAbs in the dataset.

Predicted human PK profiles using ‘typical’ human popPK parameter estimates ( ), compared with observed (individual) dose normalized human PK data (symbols) and profiles from fitting human data for individual mAbs ( ).

Use of allometric exponents estimated from Tg32 mouse or cynomolgus monkey

For every mAb in the dataset with both Tg32 mouse data and human data (n=8), mouse 2- compartment PK parameters were scaled to human using the allometric exponents estimated for Tg32 mouse (Table 4). This process was also completed for every mAb in the dataset with both cynomolgus monkey and human data (n=16). Simulations of the scaled PK parameters of CL, Q, V1 and V2 estimated from the Tg32 mouse dataset, cynomolgus monkey dataset or the population values estimated from the human dataset were compared with observed clinical PK for individual mAbs (Fig. 7). RMSEs between observed and predicted data are shown in Tables 6 and 7 for Tg32 mouse and cynomolgus monkey, respectively. As shown in Table 8, RMSEs were compared across different prediction methods: use of human popPK parameters or allometric exponents estimated for Tg32 mouse or cynomolgus monkey. The best prediction method (determined by the lowest RMSE) is indicated for each mAb.

Figure 7.

Predicted human PK profiles using (a) allometric exponents estimated for Tg32 mouse ( ), (b) allometric exponents estimated for cynomolgus monkey ( ) and (c) ‘typical’ human popPK parameter estimates (). Simulated data is compared with observed (individual) dose normalized human PK data (symbols) and profiles from fitting human data for individual mAbs ( ).

Table 6.

Comparison of observed human PK parameters for individual mAbs and predicted human pharmacokinetic parameters determined by allometric scaling from hFcRnTg32 mouse.

mAb Number	Human PK	CL mL/h/kg	V1 mL/kg	Q mL/h/kg	V2 mL/kg	%RMSE
mAb18	Observed	0.46	66.01	0.19	166.77	114.4
mAb18	Predicted	0.45	37.34	0.51	56.69
mAb19	Observed	0.46	52.19	0.18	5.98	47.4
mAb19	Predicted	0.28	44.71	0.47	43.87
mAb20	Observed	0.14	39.56	0.38	20.51	66.6
mAb20	Predicted	0.14	57.46	0.61	36.56
mAb21	Observed	0.1	40.57	0.74	13.33	60.9
mAb21	Predicted	0.16	48.85	0.5	33.85
mAb22	Observed	0.05	41.7	0.29	42.72	83.1
mAb22	Predicted	0.1	51.2	0.8	42.49
mAb23	Observed	0.11	54.92	0.33	49.14	137.6
mAb23	Predicted	0.06	27.42	0.14	26.41
mAb24	Observed	0.07	31.39	0.21	22.72	24.3
mAb24	Predicted	0.08	24.93	0.22	28.29
mAb25	Observed	0.32	55.83	0.32	27.63	47.7
mAb25	Predicted	0.22	59.27	0.43	33.05

CL: clearance from the central compartment, V1: volume of the central compartment, Q: inter-compartment distribution clearance, V2: volume of the peripheral compartment.

%RMSE: percent root mean square error.

Table 7.

Comparison of observed human PK parameters for individual mAbs and predicted human pharmacokinetic parameters determined by allometric scaling from cynomolgus monkey.

mAb number	Human PK	CL mL/h/kg	V1 mL/kg	Q mL/h/kg	V2 mL/kg	%RMSE
mAb5	Observed	0.26	54.51	0.3	32.97	143.5
mAb5	Predicted	0.1	27.09	0.37	21.65	.
mAb7	Observed	0.21	57.02	0.69	16.43	43.5
mAb7	Predicted	0.28	58.88	0.83	70.31	.
mAb9	Observed	0.12	46.58	0.37	40.43	20.9
mAb9	Predicted	0.09	51.18	0.15	60.19	.
mAb10	Observed	0.07	34.62	0.25	24.38	36.7
mAb10	Predicted	0.11	46.63	0.17	27.36	.
mAb12	Observed	0.17	40.49	0.25	32.75	78.4
mAb12	Predicted	0.08	36.37	0.18	40.05	.
mAb14	Observed	0.18	38.52	0.31	15.72	50
mAb14	Predicted	0.29	40.93	0.34	18.56	.
mAb16	Observed	0.1	40.23	0.36	42.43	23.7
mAb16	Predicted	0.12	35.98	0.28	38.08	.
mAb17	Observed	0.09	48.3	0.19	64.89	135.2
mAb17	Predicted	0.07	38.09	0.44	13.91	.
mAb18	Observed	0.46	65.98	0.2	167.69	40
mAb18	Predicted	0.58	72.13	0.21	134.53	.
mAb19	Observed	0.46	52.18	0.19	6.39	67.9
mAb19	Predicted	0.26	56.93	0.02	15.17	.
mab20	Observed	0.14	39.53	0.39	20.51	55
mAb20	Predicted	0.1	42.3	0.3	33.58	.
mAb21	Observed	0.1	40.07	0.8	14	61.3
mAb21	Predicted	0.13	47.35	0.51	54.07	.
mAb22	Observed	0.05	41.67	0.29	42.73	146.5
mAb22	Predicted	0.26	54.23	0.57	61.06
mAb23	Observed	0.11	54.96	0.33	49.28	46.7
mAb23	Predicted	0.08	43.34	0.21	38.39	.
mAb24	Observed	0.07	31.37	0.21	22.75	31.9
mAb24	Predicted	0.08	35.94	0.14	28.98	.
mAb25	Observed	0.32	55.76	0.32	27.73	237.2
mAb25	Predicted	0.09	38.48	1.39	13.08	.

CL: clearance from the central compartment, V1: volume of the central compartment, Q: inter-compartment distribution clearance, V2: volume of the peripheral compartment.

%RMSE: percent root mean square error.

Table 8.

Comparison of model prediction RMSE (%) across different methods: use of human PopPK parameters, allometric scaling from Tg32 mouse and allometric scaling from cynomolgus monkey. The best prediction method (determined by lowest RMSE) is indicated for each mAb.

	Prediction method/ RMSE %
mAb Number	Human PopPK	Tg32 mouse allometric scaling	Cyno allometric scaling	Best
mAb5	50.6	.	143.5	Human
mAb6	23.1	.	.	.
mAb7	31.4	.	43.5	Human
mAb8	.	.	.	.
mAb9	24.7	.	20.9	Cyno
mAb10	53.8	.	36.7	Cyno
mAb11	.	.	.	.
mAb12	18.2	.	78.4	Human
mAb13	.	.	.	.
mAb14	32.8	.	50	Human
mAb15	.	.	.	.
mAb16	41.6	.	23.7	Cyno
mAb17	54.2	.	135.2	Human
mAb18	361.6	114.4	40	Cyno
mAb19	121.9	47.4	67.9	Tg32
mAb20	33.5	66.6	55	Human
mAb21	54.8	60.9	61.3	Human
mAb22	111.8	83.1	146.5	Tg32
mAb23	39.7	137.6	46.7	Human
mAb24	78.4	24.3	31.9	Tg32
mAb25	104.1	47.7	237.2	Tg32

Discussion

In this study, we report a meta-analysis of the linear PK of mAbs across different species used in the pharmaceutical industry. Data on Pfizer mAbs were available from historical studies in human or cynomolgus monkey. In addition, for a subset of mAbs, PK data were available in transgenic mice expressing the human neonatal Fc receptor (Tg32 homozygous hFcRn mice). The hFcRn Tg32 mouse model was chosen over wild type (WT) mouse because mAb PK is often variable in WT mouse with poor predictability to human, which may be due to species differences in binding of human mAbs to mouse FcRn. , Following cellular uptake of mAbs exhibiting linear CL, FcRn functions as a salvage receptor to protect IgG from rapid intracellular catabolism. The hFcRn transgenic mouse (Jackson Laboratory) is null for the α chain mFcRn and contains 1 or 2 transgenes of hFcRn, hemizygous or homozygous, respectively. The homozygous Tg32 strain is used in this study with an hFcRn promoter. Avery et al, show that mAb CL in hFcRn homozygous mouse correlates with human PK (r2=0.83, r=0.91) better than NHP (r2=0.67, r=0,82).

The long half-life of endogenous and exogenous IgG in conjunction with small volumes of distribution results in the characteristic bi-exponential decline of mAb concentration-time profiles following intravenous (IV) administration. As such, a 2-compartmental PK model was shown to best describe the data. PopPK analysis is a useful tool that has been used for mAbs to quantify typical disposition characteristics and sources of variability within study populations. , The advantage of popPK analysis is that it can be used to simultaneously evaluate PK data from all studies and individuals available. We made use of it in this analysis to create a dataset representing a range of mAbs with linear CL, but differences in number of doses, subjects, data points and subject characteristics. A review of the literature performed by Dirks et al. showed that the popPK of different mAbs was similar despite differences in their pharmacological target and the fact that they were studied in different patient populations and disease states. PopPK analysis is often used to study the inter-subject variability of mAb PK and to explore covariates of this variability. Body weight/ surface area are the most commonly identified covariates found to influence the PK of mAbs. , , , The effects of other demographic factors, including age, sex, ethnicity, body size, genetic polymorphisms, concomitant medications, immune status and multiple other patient-specific details, have also been considered.

In this analysis, popPK was used to understand and quantify the variability in linear mAb PK. The inter-individual variability represents both inter-mAb and inter-subject variability, as we sought to estimate a single set of PK parameters across mAbs.

The popPK parameters estimated are presented in Table 2. These values represent typical 2-compartment PK estimates for mAbs with linear elimination in human, cynomolgus monkey and Tg32 mouse. Knowledge of typical parameter estimates of linear mAb PK a priori can be very useful in understanding and optimizing the PK/PD of a therapeutic mAb. They can be used at early stages to simulate the behavior of a mAb with ‘typical’ PK and to inform project teams on the benefit of extending PK half-life by altering affinity to FcRn. In animal PK studies, they can be used to inform initial parameter estimates for PK/PD modeling and for simulations to optimize in vivo study designs. The parameters can be used as part of a strategy to predict PK in the clinic (as will be discussed below). They could also be used as informative priors for a Bayesian data analysis or to construct parameter uncertainty distributions for clinical trial simulations.

In order to put the human popPK parameter estimates for our dataset into context, they were compared with population estimates for therapeutic mAbs reported to have linear PK in the scientific literature (Table 3). Five relevant mAbs were found including bevacizumab, pertuzumab and trastuzumab, which are humanized mAbs, and infliximab and rituximab, which are chimeric mAbs. In addition, different popPK estimates were considered for infliximab in patients with ankylosing spondylitis (AS) and UC.

Estimates of V1 and V2 were similar for the 5 different mAbs and consistent with the population PK parameters estimated herein (Table 3). In a review of therapeutic mAb popPK parameters reported by Dirks and Meibohm, which included the 5 mAbs in Table 3, the estimate of V1 was 3.1 (2.4– 5.5) L, which is equivalent to 44.3 (34.3– 78.6) mL/kg assuming a 70 kg body weight in human, and very similar to the popPK estimate of V1 in this report of 46.3 (45.1– 47.5) mL/kg. The value of the volume of distribution at steady state (Vss) can be calculated from the sum of V1 and V2. The population estimates for V1 (46.3 mL/kg) and V2 (31.5 mL/kg) in this analysis result in a calculated Vss of 77.8 mL/kg or, 5.4 L assuming a 70 kg individual. Data from early IgG metabolism studies in humans indicates that the mean serum IgG concentration and the total body IgG pool were 12 g/L and 1.06 g/kg, respectively, which for a 70 kg person equates to a volume of distribution of 6.2 L for IgG. In summary, the estimates of both central and peripheral volumes in the analysis reported herein appear to be consistent with published popPK of therapeutic mAbs with linear PK. In addition, calculated Vss is close to endogenous IgG.

The values of CL for the 5 mAbs from the literature varied from 0.12– 0.24 mL/h/kg and were similar to the popPK estimate for CL (0.15 (0.14– 0.16) mL/h/kg). The CL of infliximab in UC patients (0.24 mL/h/kg) was out of the popPK range and different to the CL of infliximab in AS patients (0.16 mL/h/kg). Variation in infliximab CL across patient populations could be due to a number of different factors. Elevated inflammatory status in UC could contribute to higher CL due to a higher whole body turnover rate and increased nonspecific proteolytic degradation rate. , In addition, CL of infliximab has also been shown to be affected by concomitant medication, as well as immunogenicity. In the analysis by Dirks and Meibohm, the population CL value for therapeutic mAbs ranged between 0.2–0.5 L/day (0.12– 0.3 mL/h/kg assuming a 70 kg individual). This range encompasses the popPK estimate of CL reported here (0.15 mL/h/kg). In addition, assuming a volume of distribution of 6.2 L and an elimination half-life of 21 days for IgG (which doesn't include IgG3), the CL of endogenous IgG is ∼ 0.21 L/day or 0.125 mL/h/kg.

The population estimate of the inter-compartmental CL (Q) was in general consistent with the 5 literature mAbs, again with exception of infliximab for UC. The median Q estimated in the Dirks and Meibohm analysis was 0.79 L/day (0.47 mL/h/kg). Our value was slightly lower at 0.27 (0.25–0.3) mL/h/kg, but both are consistent with slow transfer of mAbs between the central and peripheral compartments.

Closer inspection of the human analysis indicated that 12 of the 18 mAbs studied had individual CL values within 2-fold of the population estimate of CL (0.075– 0.3 mL/h/kg). Of the remaining 6 mAbs, 3 had CL < 0.075 mL/h/kg (mAb10, mAb22 and mAb24) and 3 had CL >0.3 mL/h/kg (mAb18, mAb19 and mAb25). Unexpected high CL of mAbs is particularly undesirable as it can lead to an increased clinical dose requirement and may limit clinical utility. Higher CL of the mAbs in this dataset was not associated with an alteration in FcRn binding, as all mAbs in this dataset had similar FcRn Kd values.

One potential mechanism that can contribute to faster than expected CL of mAbs is off-target binding. The mechanism of this non-specificity has not been fully elucidated, but could be in part due to hydrophobicity/ positive charged patches on mAbs. To test this, an in vitro high throughput assay measuring non-specific binding was implemented and used to identify mAbs with increased risk of having fast CL in humans. An AC-SINS assay was used to screen 11 of the mAbs in the dataset. An AC-SINS score of >11 has been associated with high self-association. The AC-SINS scores for the subset of mAbs tested from this dataset ranged from 0 to 24, with 3 mAbs exhibiting higher than typical CL having scores of 11–24. A correlation was observed with AC-SINS score and mAb CL across all species (Fig. 5). This trend has been observed previously with a larger dataset of mAbs and suggests that the AC-SINS assay is a useful screening tool to de-select mAbs that have the potential for fast clearance in humans.

Predicting human PK

Previous predictions of linear PK of mAbs in the clinic have been completed by single species scaling from cynomolgus monkey assuming allometric principles. Many examples of this type of scaling are available in the literature. The first comprehensive review of mAb CL prediction was completed by Ling et al in 2009. Their study of 14 mAbs indicated that, for mAbs with linear kinetics, CL in humans could be reasonably predicted from monkey data using simplified allometry with a fixed exponent. The optimal exponents were estimated to be 0.85 for soluble antigens and 0.9 for membrane-based antigens. In a similar analysis of 13 mAbs with linear CL, Deng et al showed that simple allometric scaling of CL in cynomolgus monkey with an exponent of 0.85 provided a good estimate of human CL. Dong et al also concluded that single species monkey PK predicted human PK of mAbs with linear CL within 2.3 fold. Oitate et al demonstrated that both human CL and Vss could be predicted reasonably well from monkey data alone using simple allometry with exponents of 0.79 on CL for soluble target antigens and 0.96 on CL for membrane target antigens. The exponent for prediction of Vss was close to 1 in each case.

In all examples cited, the allometric exponent for prediction of mAb CL is > 0.75, which is the standard exponent value used for interspecies scaling of small molecule drugs. This value was derived from the observation that basal metabolic rates and passive renal filtration could be scaled by body weight with an exponent of 0.75. This has been qualified by several groups in the interspecies scaling for prediction of small molecule drug CL. , Given the mechanisms of CL of small molecules are governed by oxidative metabolism and renal CL, an exponent of 0.75 makes sense. For mAbs, CL is driven by proteolysis a,nd therefore allometric exponents may be more dependent on proteolytic rates across species.

In this study, datasets of Tg32 mouse and human PK, cynomolgus monkey and human PK, and all three species PK were combined to estimate allometric exponents to scale preclinical data to human. The estimated allometric exponent for scaling Tg32 mouse CL to human was 0.90 (0.88– 0.92), cynomolgus monkey CL to human was 0.81 (0.77– 0.85) and all preclinical data CL to human was 0.89 (0.87– 0.91). The exponents required to predict cynomolgus monkey data to human were in agreement with the literature examples. The hFcRn Tg32 mouse data were encouraging as it predicted well to human, and provides a potential species to replace cynomolgus monkey for human PK predictions of mAbs. Volumes of distribution from central and peripheral compartments in general scaled with an allometric exponent of 1, similar to literature analyses. This was the first time that popPK methods had been used to estimate allometric exponents as a parameter within the model. This is a useful method because it facilitates separation of true parameter estimates from variability.

Another benefit of this analysis was that the use of 2-compartmental PK models enabled prediction of human PK profiles as well as PK parameters. The human predicted profiles were compared with 2 -compartmental fits to the observed human data for each mAb and RMSEs were calculated between observed and predicted parameters to assess prediction accuracy. Two approaches were taken. First, the popPK parameter estimates from the human analysis were used to simulate a ‘typical’ human profile. This was compared with profiles generated for each of the 18 individual mAbs in human. This method gave good prediction accuracy, with 14/18 mAbs with RMSEs <100%, indicating that the human popPK parameters were able to predict the human PK for the majority of mAbs in the dataset. This approach assumes that the CL mechanisms for these therapeutic mAbs are the same as each other and endogenous IgG, i.e. non-specific linear catabolic CL. Three of the 4 mAbs with RMSEs greater than 100% had high AC-SINS scores, indicating potential for rapid CL due to non-specific binding. The second approach was to predict human PK using allometric exponents determined from Tg32 mouse and cynomolgus monkey. For this approach, CL mechanisms need not be the same as each other or endogenous IgG, but must be consistent across species. Again, this method gave good prediction of human PK with RMSEs between observed and predicted data < 100% for 6/8 of the mAbs scaled from Tg32 mouse to human and 12/16 of the mAbs scaled from cynomolgus monkey to human.

In summary, good prediction accuracy was obtained using human ‘typical’ popPK parameters as an estimate of human PK or via scaling using allometric exponents from Tg32 mouse or cynomolgus monkey (Table 8). All of these methods are simple and easy to use. The use of allometric exponents from cynomolgus monkey is a common approach to linear mAb PK prediction, which is widely understood across the pharmaceutical industry. However, the use of human popPK parameters as base case scenario is more statistically informed than the allometric scaling approach. This is because the human popPK analysis is informed by rigorous analysis of rich datasets both from both this study and literature studies on the popPK of mAbs in human. Single species scaling using allometric exponents requires preclinical PK on individual mAbs, which is often only generated in low ‘n’ of 2–3 animals. Projection therefore relies on a small number of animals whose profile may be affected by many variables, e.g., assay specificities, inter-animal variability and immunogenicity. The strategies presented herein call into question the value of completing extensive in vivo preclinical PK for mAbs with linear CL and encourage refinement of PK strategies consistent with the ‘3Rs’, i.e., the reduction, refinement and replacement of animal use in research, testing and teaching This analysis provides alternatives to the use of cynomolgus monkey for PK prediction, including allometric scaling from Tg32 mouse, or use of human popPK parameters as a replacement to animal-based methods. As such, it has the potential to reduce the numbers of cynomolgus monkey PK studies completed. Use of the AC-SINS assay can also be used to screen out mAbs with high CL due to non-specific binding/self-association.

MAbs with non-linear PK due to target mediated CL mechanisms were not included in this analysis. Modeling of such data is routinely performed using a Michaelis-Menten (M-M) model with linear first-order elimination complemented by a non-linear pathway described using Vmax and Km parameters. Alternatively, a more mechanistic TMDD model that incorporates target properties can be used. Previous work has shown that the TMDD approach is more reliable for human projections because it can capture differences in target properties between species and disease populations. , , Since both the M-M and TMDD models require parameterization of the linear CL pathway, the parameters presented herein are a useful guide. Such models have a tendency for over-parameterization, and prior knowledge of typical linear PK parameters across species reduces the need for their estimation.

In conclusion, the popPK analysis completed on 27 Pfizer mAbs in human, cynomolgus monkey or Tg32 mouse showed that a single set of typical linear PK parameters could be estimated across species. These parameters will be useful to inform initial parameters for PK/PD modeling and for simulations to optimize in vivo and first in human study designs. In addition, different translational strategies were investigated for prediction of human linear PK of mAbs. Use of ‘typical’ human PK parameters gave good prediction accuracy for the majority of the mAbs in this study. Allometric exponents were estimated within the popPK model and also gave good predictions, from both Tg32 mouse or cynomolgus monkey to human. The strategies presented herein offer methods to predict linear human PK of mAbs with less reliance on cynomolgus monkey PK and use of smaller animal or in silico alternatives.

Materials and methods

mAb PK dataset

For this study, a dataset was compiled consisting of in-house historical individual concentration versus time data following IV administration of 27 Pfizer mAbs in human, cynomolgus monkey or hFcRn Tg32 transgenic mice. All procedures performed in animals were in accordance with established guidelines and were reviewed and approved by Pfizer's Institutional Animal Care and Use Committee. Research on human samples was conducted in accordance with all applicable Pfizer policies, including institutional review board / independent ethics committee approval. Data for 18/ 27 mAbs were available in healthy human volunteers or patients and consisted of single dose IV PK, at multiple dose levels with n = 3–24 individuals/ dose level. In cynomolgus monkey, single dose IV PK data were available for 23/ 27 mAbs, administered at 1–3 dose levels with n=2 monkeys/ dose. In hFcRn Tg32 transgenic mice, 11 /27 mAbs were administered as previously described at a single IV dose of 3.5 mg/kg (1 mAb) or 5 mg/kg (10 mAbs) with n=5–6 mice per mAb. For the marketed mAbs, popPK estimates were obtained from the literature.

Data inclusion and exclusion criteria: Selection of linear dose range

Non-linear data were removed from the datasets, where TMDD or factors such as immunogenicity were contributing to the overall CL. First, a visual analysis of the data was performed to check for non-linearity, then an algorithm was applied to a linear regression of dose and AUC (Fig. 1) to test for deviation of the slope from 1, and the dataset was reduced accordingly. For Tg32 mouse, all PK data were at doses ≥3.5 mg/kg and were included in the analysis. Table 1 provides details of the linear dose range, number of dose levels that were used in the popPK analysis and the non-linear dose levels removed.

PK model

mAb PK following IV administration was described by a 2-compartment disposition model with first-order elimination from the central compartment. The structural model was parameterized in terms of CL, central volume of distribution (V1), peripheral volume of distribution (V2) and inter-compartmental clearance (Q). Random effects were included as exponential terms reflecting log normal distributions of model parameters. The residual variability was implemented by proportional error model per compound. Goodness-of-fit was determined using the minimum value of the objective function defined as minus twice the log likelihood. For nested models, a decrease of 3.84 points in the objective function (MVOF; corresponding to P < 0.01 in a chi-squared distribution) by adding an additional parameter was considered significant. The goodness-of-fit was also investigated by visual inspection of the plots of individual predictions and the diagnostic plots of (weighted) residuals.

The PK model was applied to each species separately (single species PK analysis) or to a combination of two or three species together (combined species PK analysis). The combined species PK analysis used allometric scaling on all PK parameters by a scaling exponent based on bodyweight. For humans the available bodyweights were used, for cynomolgus monkey and hFcRn Tg32 transgenic mice a bodyweight of 3 kg and 0.02 kg, respectively, was assumed.

Computation

Individual concentration-time data from all subjects for all mAbs were pooled into a single dataset for pop PK analysis using nonlinear mixed effect modeling with the NONMEM software system (Version 7.3, ICON Development Solutions) using ADVAN3 and subroutine 4 and PsN (version 4.6.0). Gfortran version 4.6.0 was used as compiler. Parameters were estimated using the first-order conditional estimation method with interaction between the two levels of stochastic effects (FOCEI).

Model validation

The robustness of the final combined species PK model was evaluated using resampling techniques of a bootstrap method. The bootstrap method involves repeated random sampling of subjects in the dataset. The original dataset is replaced to produce another dataset of the same size, but with a different combination of subjects and compounds. Resampling was repeated 200 times. The obtained PK parameters from the bootstrap sets that produced successful minimization and convergence were used to simulate the concentration time profiles after a single dose of 1 mg/kg for a typical human subject of 70 kg.

In the visual predictive check (VPC), the observation versus time profile was simulated 1000 times by means of Monte Carlo simulations. In a Monte Carlo simulation, random values are drawn from the distributions of the identified random effects. Subsequently, the median and 5 and 95 percentiles of the dependent variables were calculated for each time and plotted together with the observations. A jackknife analysis was also completed with the combined species PK model to test robustness of model predictions. In this type of analysis, one mAb is omitted from the total dataset at a time to test the influence of that mAb on the resulting predictions.

AC-SINS

The affinity capture self-interaction nanoparticle spectroscopy method was implemented as previously described.21,32

PK prediction methods

To predict the human PK of mAbs in this dataset:

Using ‘typical’ PK parameters for human

Clinical data, available for 18 mAbs in the dataset, were fit individually using a 2 compartment PK model as described previously. Individual values of CL, V1, Q and V2 were determined for each mAb, and PK profiles were simulated using fitted parameters. These were compared to simulated profiles using the human popPK values of CL, V1, Q and V2 and root mean square error (RMSE) values calculated using the following equation:Where n= number of predictions, χ is predicted or observed value

Using allometric exponents estimated from Tg32 mouse or cynomolgus monkey

For every mAb in the dataset with both Tg32 mouse data and human data (n=8), or cynomolgus monkey and human data (n=16), preclinical species PK parameters were scaled to human using estimated allometric exponents (Table 4).

First, the individual Tg32 mouse data for each mAb were fit to a 2-compartment PK model (n=8) as described previously. Individual mAb values of CL, V1, Q and V2 were then scaled to human using the estimated allometric exponents for Tg32 mouse presented in Table 4. PK profiles were simulated using the scaled parameters. These were compared with the individual values of CL, V1, Q and V2 estimated for each individual mAb in human, described in section 1 above, and RMSEs calculated as before.

This process was also completed for every mAb in the dataset with both cynomolgus monkey and human data (n=16).