Next Article in Journal
Functional Activity of Cytokine-Induced Killer Cells Enhanced by CAR-CD19 Modification or by Soluble Bispecific Antibody Blinatumomab
Previous Article in Journal
Dynamics of IgM and IgG Antibody Response Profile against Linear B-Cell Epitopes from Exoerythrocytic (CelTOS and TRAP) and Erythrocytic (CyRPA) Phases of Plasmodium vivax: Follow-Up Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Physiologically Based Pharmacokinetic Model Relates the Subcutaneous Bioavailability of Monoclonal Antibodies to the Saturation of FcRn-Mediated Recycling in Injection-Site-Draining Lymph Nodes

Simcyp Division, Certara UK Ltd., Level 2 Acero, 1 Concourse Way, Sheffield S1 2BJ, UK
*
Author to whom correspondence should be addressed.
Submission received: 9 July 2024 / Revised: 8 August 2024 / Accepted: 12 August 2024 / Published: 15 August 2024

Highlights

What are the main findings?
  • Subcutaneous dosing of monoclonal antibodies (mAbs) leads to high local drug concentration at the injection site, in the dosing site draining lymph nodes, and in antigen-presenting cells (APC), resident in peripheral lymph nodes.
  • The elevated local IgG concentration is predicted to result in a transient saturation of the FcRn recycling pathway in APCs and consequently to an increased degradation of macropinocytosed IgG, which is predominantly composed of the dosed mAb.
  • Elevated degradation of macropinocytosed IgG during the first pass of the subcutaneously dosed mAb was predicted to manifest in an about 70% bioavailability of the drug.
What is the implication of the main finding?
  • The bioavailability of mAbs can mechanistically be predicted from in vitro FcRn binding data and local concentration in the draining lymph node of the subcutaneous dosing site.
  • The model can be used to understand the role of non-specific uptake of IgG in the dosing site draining lymph nodes and its impact on the bioavailability of monoclonal antibodies.

Abstract

:
The bioavailability of a monoclonal antibody (mAb) or another therapeutic protein after subcutaneous (SC) dosing is challenging to predict from first principles, even if the impact of injection site physiology and drug properties on mAb bioavailability is generally understood. We used a physiologically based pharmacokinetic model to predict pre-systemic clearance after SC administration mechanistically by incorporating the FcRn salvage pathway in antigen-presenting cells (APCs) in peripheral lymph nodes, draining the injection site. Clinically observed data of the removal rate of IgG from the arm as well as its plasma concentration after SC dosing were mostly predicted within the 95% confidence interval. The bioavailability of IgG was predicted to be 70%, which mechanistically relates to macropinocytosis in the draining lymph nodes and transient local dose-dependent partial saturation of the FcRn receptor in the APCs, resulting in higher catabolism and consequently less drug reaching the systemic circulation. The predicted free FcRn concentration was reduced to 40–45%, reaching the minimum 1–2 days after the SC administration of IgG, and returned to baseline after 8–12 days, depending on the site of injection. The model predicted the uptake into APCs, the binding affinity to FcRn, and the dose to be important factors impacting the bioavailability of a mAb.

Graphical Abstract

1. Introduction

Subcutaneous (SC) administration is increasingly used for monoclonal antibodies (mAbs) and other therapeutic proteins (TPs) to improve patient compliance and reduce treatment costs [1]. One challenge in the study of mAbs is understanding, in mechanistic and quantitative terms, the factors that affect their bioavailability, which reflects the loss of a fraction of the injected drug before reaching the systemic circulation. There is no sufficiently reliable in vitro assay or preclinical species that is predictive of mAb bioavailability in humans. It is known that bioavailability correlates with the systemic clearance of a mAb [2], which is a measure of intrinsic stability, hydrophobicity, and any charge effect [3,4]; however, this relationship is empirical and might not be predictive for novel mAbs. The lack of mechanistic modeling of bioavailability was recently identified as one of the key knowledge gaps by the SC consortium [5].
Macrophages and other APCs of monocyte lineage have been demonstrated to be important for bioavailability in rodents [6,7]. They are mainly located in lymph nodes, including those draining the SC site; although, macrophages are also present in the dermis of the skin and can migrate to the injection site as a reaction to the mechanical pressure of the injection needle [8]. Macrophages contain the neonatal Fc receptor (FcRn), which protects immunoglobulin G (IgG) and other Fc-containing TPs from lysosomal degradation [6,7]. APCs endocytose large volumes of surrounding fluid [9] whilst scanning the environment for foreign antigens as part of the innate immune system [10,11].
Physiologically based pharmacokinetic (PBPK) models can be expected to predict the bioavailability of novel mAbs or altered bioavailability in disease populations once the respective mechanisms are formulated and incorporated into these models. However, at present, mechanistic a priori predictions of bioavailability are challenging [5], because only descriptive pre-systemic clearance pathways have been fitted to observed data and these are not predictive for novel mAbs. It has been found that lymphatic clearance is correlated to bioavailability rather than clearance at the injection site [12,13,14], which is in accordance with preclinical data [6,7].
Another challenge exists in defining the SC injection site volume in a PBPK model, because its dimensions are, to an extent, arbitrary, as there is no hard boundary between exposed and unexposed regions of the skin. An arbitrary value of 5 mL has often been used in the past [15], which assumes a three- to five-fold dilution of the typical SC dose of 1–1.5 mL. The volume of the injection site is determined by the distribution space between the site of injection and the blood vessels in case of peptides, small proteins, and oligonucleotides; or between the site of injection and the initial lymphatics in case of mAbs and other larger TPs [16]. The pre-systemic clearance could depend on the administration site because of different physiological characteristics such as lymph flow, the number of draining lymph nodes, or the number of macrophages in different parts of the skin [17].
The primary aim of this study was to develop a mechanistic PBPK model for estimating mAb bioavailability that includes the local lymphatic system with APCs in SC-injection-site-draining lymph nodes. In addition, the properties of four different injection sites (the arm, abdomen, back, and thigh) were considered to simulate site-specific differences in mAb pharmacokinetics and bioavailability.

2. Material and Methods

Four steps were taken in the construction of this SC PBPK model. First, a literature review was conducted to understand the physiological processes occurring following SC dosing. Second, the model was developed and parameterized with published physiological data. Third, a case study was conducted with IgG to test the model performance without optimization of any physiological or drug parameters. And fourth, we explored drug-related factors that have an impact on bioavailability.

2.1. Model Structure

The structure of a previously published SC injection site model was enhanced by the inclusion of local lymph capillaries and a peripheral draining lymph node compartment next to the injection site (Figure 1; [15]). APCs were incorporated as an additional compartment in the peripheral lymph node. The processes in APCs included fluid-phase endocytosis, FcRn recycling, and catabolism of unbound mAbs. Endogenous IgG was included during the simulation [18]. The SC site model was coupled to a whole-body PBPK model, implemented in the Simcyp Simulator V22 (Certara UK Ltd., Sheffield, UK). The SC site had a permeability-limited structure, similar to other compartments of the PBPK model. It contained the vascular space, the endothelial endosome, and the interstitial space. Dose administration was included in the interstitial space and a bolus injection was always simulated. The two-pore filtration–diffusion description was used to describe the distribution of dosed TPs between the vascular and interstitial spaces of the injection site [19]. The FcRn pathway was incorporated into the endothelial endosomal compartment [18]. The endothelial cell volume was parameterized as 5% of the vascular volume and endosomes covered 25% of the cell volume [20,21]. The ordinary differential equations to describe the disposition of the mAb at the SC site can be found in the Supplementary Materials.

2.2. Model Parameterization

First, the SC site volume was defined by using a study that investigated the travel distance of labeled albumin between the arm/back and sentinel lymph nodes in humans [22]. Four individuals (two women) aged 45–70 years had a mean (min–max) travel distance of 19.6 (11.5–32.5) cm and a mean (min–max) travel time of 2.0 (1.0–2.3) minutes after the injection of labeled albumin into the arm. Measurements in the back were conducted in six individuals (four women) aged 31–72 years. The mean (min–max) travel distance was 23.9 (11.5–40.0) cm, and the mean travel time (min–max) was 4.4 (1.0–10.5) minutes. The high variability observed in that study could not be correlated to age or sex. Given that the distance was measured to sentinel lymph nodes, these distance values will likely include the initial lymphatics. By assuming a circular distribution radius and by taking a subcutaneous thickness of 1.91 mm [23], the SC site volume could be estimated (Table 1). An average value of arm-and-back injection site volume was assumed for the abdomen and thigh in the absence of measured values for these two sites. The fraction of the injection site volume accounted for by the vascular, endosomal, and interstitial compartments was assumed to be the same as that accounted for by the skin.
Second, the blood flows were parameterized. Blood flows were defined as 9.2 mL/min/100 g, 2.9 mL/min/100 g, 7.4 mL/min/100 g, and 6.0 mL/min/100 g for the arm, abdomen, back, and thigh, respectively [23]. These values were incorporated as a % of cardiac output by using the derived SC site volume, a skin density of 1.1 g/mL [24], and the cardiac output of a population representative of the Simcyp Simulator V22 (397.5 L/h).
Third, the local lymphatic system was incorporated into the PBPK model. An important parameter is the number of peripheral lymph nodes that a mAb has to distribute through before reaching the central lymphatic system and consequently the blood circulation. This number of lymph nodes was found to be in the ranges of 25–36, 18–26, 20–29, and 18–26 for the arm, abdomen, back, and thigh, respectively [25]. The study that was used to measure the travel distance and time of labeled albumin to sentinel lymph nodes was used to determine the afferent lymph flow. The mean flow in the arm and back was 4.80 × 10−5 L/h and 4.87 × 10−5 L/h, respectively [22]. Each mean flow value was converted to the % of total lymph flow by multiplying it with the number of lymph nodes and using the value of total lymph flow from the Simcyp Simulator V22 (0.391 L/h). The average arm-and-back lymph flow was used for the abdomen and the thigh in the absence of measured values for these sites. The efferent lymph flow was estimated, using a mathematical model, to be 93% of the afferent lymph flow in mice [26]. The remaining 7% of the afferent lymph flow travels to the systemic circulation at the peripheral lymph nodes. Convection through aquaporins was incorporated to maintain fluid mass balance between afferent and efferent lymph flow in the model. Aquaporins have a radius of 0.2 nm [27], which means that mAbs are excluded from distribution. These calculations are similar to the two-pore hypothesis. The lymph capillary volume was parameterized by using the aforementioned study with labeled albumin, which reported a relationship between lymph flow and the lymph capillaries; the study measured the travel distance to sentinel lymph nodes and was thus corrected by the density of lymph vessels in the skin [28]. The average volume of a cervical lymph node is 0.284 mL, with high interindividual variability [29]. This average volume was multiplied with the number of lymph nodes to obtain the volume of the peripheral SC-site-draining lymph node compartment. The peripheral lymph nodes collect the fluid from interstitial-space watershed areas to filter it and APCs remove any foreign particles or proteins [30].
Table 1. Physiological input parameters for the four different SC sites.
Table 1. Physiological input parameters for the four different SC sites.
ParameterArmAbdomenBackThighReferences
SC site volume (L/kg)0.000770.000990.001150.00099Calculated based on [22]
Vascular volume (% of the SC site)5555[15]
Endosomal volume (% of the SC site)0.060.060.060.06[20,21]
Interstitial volume (% of the SC site)59.559.559.559.5[15]
Blood flow (% of the cardiac output)0.0880.0350.1050.073[23]
Afferent lymph flow (% of the total lymph flow)0.1910.2050.2620.272[22]
Efferent lymph flow (% of the afferent lymph flow)93.093.093.093.0[26]
Lymph capillary volume (L/kg)0.0003520.0003920.0004150.000387[22]
Peripheral lymph node volume (L/kg)0.0001168.38 × 10−59.52 × 10−58.38 × 10−5[25]
It is rather difficult to measure the number of macrophages in a lymph node. In mice, 42,000 subcapsular sinus macrophages have been measured in a single lymph node [31], which represent 29.7% of all macrophages [32]. The volume of a representative macrophage and lymph node in mice has been measured to be 1.66 × 10−9 [33] and 0.00095 mL [34], respectively. This leads to a macrophage volume of 24.5% of the peripheral lymph node volume in mice, and we assumed that the macrophage volume in humans would be similar to this value (Table S1). The fluid-phase endocytosis rate into macrophages has been reported to be 25% of their volume per hour [9]. Macrophages have been reported to contribute to site-specific bioavailability differences; however, we assumed the macrophage parameters to be similar between injection sites given the limited experimental data. Other cells present in lymph nodes also contribute to macropinocytosis and bioavailability, like dendritic cells and B-cells. In that sense, the “macrophage” population in the PBPK model includes all cells which are involved in high-rate macropinocytosis and mAb recycling.

2.3. Bioavailability Calculation

Bioavailability is defined as the fraction of a dosed drug that reaches the systemic circulation. We therefore split the concentration at the injection site into one coming from the dose and one returning from the systemic circulation. The concentration coming from the dose considered only outflows from the injection site and was used to calculate the bioavailability (see Supplementary Materials). The accumulated eliminated amount of the different clearance pathways was integrated, namely catabolic clearance at the SC site in macrophages and in the central lymph node compartment.

2.4. Model Validation

The model was evaluated using clinically observed data for subcutaneously dosed mAbs [35] as well as the clearance of labeled IgG from the SC injection site [15]. Visual predictive checking was used to analyze the model’s performance for the observed vs. predicted data. The published peak concentration (Cmax) and time to Cmax (tmax) were compared against the model’s predictions and predictions were considered successful if any differences they showed were within a 1.5-fold alteration of the observed data.
Cmax and tmax for IgG were predicted at the four different SC sites and the results were compared against clinically observed data for 47 different IgG-based mAbs for which site-specific pharmacokinetic data had been published [17].
We evaluated the impact of the FcRn recycling pathway in the vascular endothelium at the SC site and in APCs by simulating four virtual patient cohorts along the lines of the preclinical study by Richter et al. [6], who studied in Tg32 transgenic mice the relative fractions of human IgG catabolized in bone-marrow-derived hematopoietic cells versus endothelial and parenchymal cells in the body. The relative changes in bioavailability were compared against the findings of Richter et al. [6].

2.5. Sensitivity Analyses

Pre-systemic first-pass catabolism in draining lymph nodes by resident APCs appears to be a critical determinant of the SC bioavailability of mAbs. Sensitivity analysis was performed on the pinocytotic uptake rate (range: 0–2.5 1/h) and on the binding affinity to FcRn (range: 0.128–2.0 μM). Additionally, sensitivity analysis was performed on the dose, which ranged from 0.01 to 1000 mg, to investigate the general trend, knowing that the used dose boundaries are not clinically relevant. All sensitivity analyses were performed using the arm as a representative SC injection site.

3. Results

Three steps were performed to evaluate the model.
First, the model performance was validated against the IgG plasma concentration after SC injection [35]. The predicted Cmax and tmax of 33.5% of the injected dose (predicted/observed ratio: 0.96) and after 5.0 days (predicted/observed ratio: 0.98) agreed with the published clinical data (Figure 2A).
Second, the removal of IgG from the SC dosing site in the arm was simulated and compared against clinically observed data [15]. The predicted average removal rate agreed with the observed data; however, the variability was slightly underpredicted (Figure 2B).
Third, the predicted impact of site-specific administration was compared against clinically observed data from 47 mAbs ([17]; Table 2). The Cmax ratio for the arm/abdomen and the thigh/abdomen agreed with the observed data (predicted/observed ratio: 0.99). The arm was predicted to have a higher tmax than the thigh in accordance with the observed data; however, the arm had the highest observed tmax, but the model predicted the highest tmax in the abdomen. The trend in the AUC was captured by the model with a predicted AUC ratio for the arm/abdomen of 0.95 (predicted/observed ratio: 0.97) and a predicted AUC ratio for the thigh/abdomen of 1.06 (predicted/observed ratio: 0.94). Overall, the site-specific differences agreed with clinical data but were predicted to be marginal based on the physiological data of the respective site.
The impact of the APC compartment in the model was also investigated. The mouse study of Richter et al. was simulated in humans with subcutaneously dosed IgG to analyze whether the general trend in bioavailability observed in transgenic mice can be captured [6]. Richter et al. had four cohorts of mice: cohort 1 with wildtype mice, cohort 2 comprising mice with FcRn knocked-out in hematopoietic cells, cohort 3 comprising mice with FcRn knocked-out in endothelial cells, and cohort 4 comprising mice with a global FcRn knock-out. We simulated the bioavailability of IgG in humans using a similar study design to test the contribution of APCs to bioavailability. The predicted bioavailability in cohorts 1–4 was 65.1%, 25.6%, 66.1%, and 24.8%, respectively. Thus, cohort 1 and 3, as well as cohort 2 and 4, share roughly the same bioavailability, similar to the preclinical findings of Richter et al. The predicted/observed ratios for the bioavailability in cohorts 2, 3, and 4 relative to cohort 1 were 0.87, 1.00, and 0.60, respectively.
In the next step, we investigated the contribution of macrophages to bioavailability in the model. After the injection of 10 mg, the mAb concentration in the interstitial space of the SC site was around 1000 mg/L but became diluted when distributing into the lymph capillaries and peripheral lymph nodes, draining the SC site. Nevertheless, even the diluted mAb concentration in APCs (550–750 mg/L, depending on the SC site) remained sufficiently high to result in a transient saturation of the FcRn recycling pathway, as manifested by the decrease in the free fraction of the receptor (Figure 3). The free FcRn concentration was predicted to drop to 40–45%, depending on the injection site. The minimal free FcRn level was reached 35 h after the injections into the arm and abdomen and this was observed half a day earlier for the back and thigh. The saturation of FcRn in the macrophages was transient and the free FcRn was predicted to return to baseline after 8–12 days, depending on the SC site. The transient saturation of FcRn in the SC-dosing-site-draining lymph nodes results in elevated catabolism in macrophages and reduces the bioavailability to below 100%.
Sensitivity analyses were conducted investigating the impact of the endocytosis uptake into macrophages, the binding affinity to FcRn, and the administered dose on bioavailability (Figure 4). At a default uptake rate of 0.25 1/h, a binding affinity to FcRn of 0.728 μM, and a dose of 10 mg, the model predicted the bioavailability to be 65.1% for IgG injected into the arm. If the uptake rate was two-fold lower, the predicted bioavailability increased to 79.2%, and if the uptake rate was two-fold higher, the predicted bioavailability decreased to 47.3%. If the KD value for IgG (0.728 µM) was two-fold lower, the predicted bioavailability increased to 76.1%, and if the KD value was two-fold higher, bioavailability decreased to 53.6%. Dose had a marginal impact on bioavailability up to 2 mg. The predicted bioavailability decreased by 4.1% and 24.7% when the dose was increased to 10 mg and 100 mg, respectively. At a very high and theoretical dose of 1 g, bioavailability was predicted to be 32.9%.

4. Discussion

SC injection is the preferred route of administration for mAbs, but a priori predictions of bioavailability remain challenging. In this study, we evaluated and modeled the role of APCs in peripheral lymph nodes draining the SC site and demonstrated that a high drug concentration after injection transiently saturates the FcRn recycling pathway, resulting in a higher fraction of the endosomal mAb being degraded. This results in the bioavailability for mAbs being in the range of 60–70% of the dose given.
Predicting the bioavailability of novel mAbs from first principles is the next step for PBPK models, wherein the current models predominantly rely on fitted clearances and empirical relationships. Many different approaches have been taken to explore the impact of factors on bioavailability. Li et al. developed a whole-body PBPK model that included an SC site and fitted a degradation rate to represent catabolism at the injection site and in the local lymphatic system [36], while Varkhede et al. and Hu et al. implemented an SC absorption compartment that included the local lymphatic system [12,13]. In the latter approach, Varkhede et al. implemented lymphatic degradation and not clearance at the injection site as the main factor affecting bioavailability [12], with Hu et al. taking a similar approach by linking the empirically fitted lymphatic clearance [13]. In both cases, the lymphatic clearance—and thus, bioavailability—was linked to the charge as the sole predictor for a mAb.
While protein charge can play a role in mAb pharmacokinetics, we decided to focus on the role of cellular catabolism as exercised by APCs in the peripheral draining lymph nodes of the SC dosing site, because charge interaction measured by heparin binding can only be used empirically in PBPK models [37] and other factors such as hydrophobicity might also have an impact [4]. Lymphatic circulation forms an integral part of antibody circulation, with lymph nodes collecting the fluid from interstitial-space watershed areas before passing through the lymph nodes, where the resident APCs constantly sample the surrounding medium for any foreign particles or proteins as part of the early innate immune response [30]. Any mAb or protein with a great enough molecular weight to enter the lymphatic system [16] will therefore travel through the peripheral lymph nodes and be exposed to APCs before reaching the systemic circulation. APCs macropinocytose large volumes of the surrounding medium (20–40% of cell volume per hour) [9] and contain FcRn [38] to recycle endogenous compounds such as albumin and IgG, and will therefore also recycle any TP with an Fc domain. In our model, we track how a subcutaneously dosed mAb at a very high initial concentration travels from the dosing site to the draining lymph nodes, where the transient spike in total IgG causes a temporary fall in the free fraction of FcRn in APCs. As a result, a significant fraction of the endocytosed mAb cannot bind the receptor for recycling and is catabolized in lysosomes. This loss is irreversible and manifests itself as decreased bioavailability (i.e., a decreased fraction of the dose entering systemic circulation), which is typically around 60–70% for therapeutic mAbs. In comparison with the previous PBPK models that demonstrated lymphatic clearance to be important [12,13], we derived the degradation pathway a priori based on the macropinocytosis rate into APCs and saturable FcRn-mediated recycling.
The link to the protein charge found in previous models [12,13,39] can be explained by the dependency between the fluid-phase pinocytosis into APCs and the charge distribution of an antibody [40], but this effect is rather weak and correlates most strongly with the charge carried by the complementary determining regions in the antibody paratope region, which contributes modestly to overall protein charge. Nevertheless, it is possible that positively charged surface areas at the variable site of an antibody can interact with negatively charged components of cell membranes such as phospholipids or the glycocalyx. The latter can be evaluated by measuring mAb retention time on heparin affinity chromatography columns, which correlates with systemic clearance [3]. Previous PBPK models have linked the heparin retention time to the fluid-phase endocytosis uptake rate into endothelial cells [37,41]. The same effect can be expected for pinocytosis into APCs. The more positive the charge of an antibody is (or the higher the relative heparin retention time), the greater the uptake into APCs. This would result in a greater transient saturation of the FcRn pathway and lower the bioavailability (Figure 4). Antibodies with enhanced uptake into APCs in SC-site-draining lymph nodes may be predisposed to local uptake during the first pass and continued uptake during the elimination phase, establishing the correlation between bioavailability and terminal half-life. In the future, in vitro methods such as heparin column chromatography could be used as a measure of adhesion-related cellular uptake of proteins and more readily implemented in predictive PBPK models.
Another important factor that determines bioavailability in our model is the binding affinity to FcRn (Figure 4). A 10-fold-higher affinity for FcRn, with everything else being equal, results in an enhanced bioavailability of 92.5%. However, the “bottleneck” is the uptake into APCs, because only what distributes into the cell can bind to FcRn.
The current model has some limitations that should be acknowledged. First, APCs were only incorporated as a pre-systemic clearance pathway and are not included directly in our model for the homeostasis of IgG. Second, the physiological parameters used in this model carry a degree of uncertainty. The SC site volume in our model is based on the distribution space of a mAb between the site of injection and the initial lymphatics. As smaller proteins and peptides can diffuse through the pores of the vascular endothelium, resulting in a tmax within hours, the modeled distribution space, and thus the modeled SC site volume, is expected to be different compared with that of an actual mAb. We used data based on labeled albumin to parameterize the SC site volume, and thus this value is applicable for large proteins, including mAbs [16,22]. Data for smaller TPs were not found and thus the current value should be used with caution. In addition, the distribution of labeled albumin was measured until the peripheral lymph node and thus the volume can be expected to already include the initial lymphatics [22]. The afferent lymph flow was determined based on labeled albumin from the same study [22]. Despite the uncertainty, the removal of IgG from the arm of breast cancer patients could be predicted well using our model without optimizing any parameters (Figure 2B), which gives confidence to the volume and lymph flow values we used for the arm. However, when comparing the different sites, the highest tmax was predicted in the abdomen but was observed after injection into the arm. tmax would mainly be impacted by the afferent lymph flow. Another physiological parameter with uncertainty is the number of APCs in lymph nodes, which is difficult to measure and was obtained from previous measurements taken in mice. Ideally, the model would require observed concentration data of IgG in the lymph nodes draining the injection site, but as these data are currently not available, we rely on the assumptions made in the model and the observed plasma concentration time-course values, which were predicted in accordance with the observed data for an IgG-based mAb. One further limitation is that the model verification beside IgG focused on the Cmax and tmax of 47 mAbs and not the concentration–time profiles of these mAbs. Cmax and tmax values embed bioavailability, a key characteristic of subcutaneous absorption, and determining this was the goal of our model evaluation. Once in circulation, there is no further significant impact of the dosing site on the mAb plasma concentration time course, and thus Cmax and tmax evaluation can be considered sufficient for the evaluation of subcutaneous absorption. A strength of the verification method we used is that plasma Cmax and tmax values can be predicted without empirical fitting of any absorption-related parameters.
Given the limitations, this model describes a typical IgG-based mAb and is most useful for evaluating the impact of factors like dose volume, concentration, and interval time, in addition to generic parameters like the FcRn-binding affinity. We also envisage the possibility that if the excipients play a dominant role in keeping the mAb soluble in the formulation, these would be lost from the dosing site more rapidly than the mAb, thus increasing the risk of aggregation and dosage loss.

5. Conclusions

In conclusion, the incorporation of APCs in the SC-dosing-site-draining lymph nodes allows prediction of the bioavailability of IgG. The reason for the reduced bioavailability of mAbs is described as a consequence of the transient saturation of the FcRn receptor in APCs, which leads to a higher catabolism and consequently a lower bioavailability. Important antibody-related factors that have an impact on the bioavailability are the uptake rate into APCs and the binding affinity to FcRn. Higher doses, a positive charge distribution, or a low binding affinity to FcRn can result in a lower bioavailability. The presented model describes the bioavailability of a mAb mechanistically and represents a transformative advancement in drug discovery and development. It offers a valuable tool that can be used to optimize dosing regimens, which will ultimately lead to safe and efficacious therapies for patients.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/antib13030070/s1: Section S1. Ordinary differential equations for the proposed subcutaneous absorption model. Section S2. Bioavailability calculation. Section S3. List of physiological parameters used to parameterize the antigen-presenting cell compartment. Table S1: Physiological parameters to parameterize the antigen-presenting cell compartment. Reference [42] is cited in the Table S1.

Author Contributions

Conceptualization, F.S. and A.S.; data curation, F.S., C.L. and A.D.; formal analysis, F.S., H.M. and X.P.; investigation, F.S., C.L., A.D., H.M. and X.P.; methodology, F.S., I.G. and A.S.; project administration, F.S.; supervision, I.G., M.J. and A.S.; validation, F.S., C.L., A.D., H.M. and X.P.; visualization, F.S.; writing—original draft preparation, F.S.; writing—review and editing, F.S., C.L., A.D., H.M., X.P., I.G., M.J. and A.S. All authors have read and agreed to the published version of the manuscript.

Funding

The activities of Simcyp Division, Certara UK, are supported by the Simcyp Consortium of pharmaceutical companies.

Data Availability Statement

Data can be made available upon request from the authors if not given in this publication.

Acknowledgments

The authors would like to acknowledge Abi Reader and Eleanor Savill for their assistance with the manuscript submission.

Conflicts of Interest

Felix Stader, Cong Liu, Abdallah Derbalah, Hiroshi Momiji, Xian Pan, Iain Gardner, Masoud Jamei, and Armin Sepp are employees of and may be shareholders of Certara.

References

  1. Viola, M.; Sequeira, J.; Seiça, R.; Veiga, F.; Serra, J.; Santos, A.C.; Ribeiro, A.J. Subcutaneous delivery of monoclonal antibodies: How do we get there? J. Control. Release 2018, 286, 301–314. [Google Scholar] [CrossRef]
  2. Zheng, Y.; Tesar, D.B.; Benincosa, L.; Birnböck, H.; Boswell, C.A.; Bumbaca, D.; Cowan, K.J.; Danilenko, D.M.; Daugherty, A.L.; Fielder, P.J. Minipig as a potential translatable model for monoclonal antibody pharmacokinetics after intravenous and subcutaneous administration. mAbs 2012, 4, 243–255. [Google Scholar] [CrossRef]
  3. Kraft, T.E.; Richter, W.F.; Emrich, T.; Knaupp, A.; Schuster, M.; Wolfert, A.; Kettenberger, H. Heparin chromatography as an in vitro predictor for antibody clearance rate through pinocytosis. mAbs 2020, 12, 1683432. [Google Scholar] [CrossRef]
  4. Grinshpun, B.; Thorsteinson, N.; Pereira, J.N.; Rippmann, F.; Nannemann, D.; Sood, V.D.; Fomekong Nanfack, Y. Identifying biophysical assays and in silico properties that enrich for slow clearance in clinical-stage therapeutic antibodies. mAbs 2021, 13, 1932230. [Google Scholar] [CrossRef] [PubMed]
  5. Sánchez-Félix, M.; Burke, M.; Chen, H.H.; Patterson, C.; Mittal, S. Predicting bioavailability of monoclonal antibodies after subcutaneous administration: Open innovation challenge. Adv. Drug Deliv. Rev. 2020, 167, 66–77. [Google Scholar] [CrossRef]
  6. Richter, W.F.; Christianson, G.J.; Frances, N.; Grimm, H.P.; Proetzel, G.; Roopenian, D.C. Hematopoietic cells as site of first-pass catabolism after subcutaneous dosing and contributors to systemic clearance of a monoclonal antibody in mice. mAbs 2018, 10, 803–813. [Google Scholar] [CrossRef]
  7. Challa, D.K.; Wang, X.; Montoyo, H.P.; Velmurugan, R.; Ober, R.J.; Ward, E.S. Neonatal Fc receptor expression in macrophages is indispensable for IgG homeostasis. mAbs 2019, 11, 848–860. [Google Scholar] [CrossRef] [PubMed]
  8. Fathallah, A.M.; Bankert, R.B.; Balu-Iyer, S.V. Immunogenicity of subcutaneously administered therapeutic proteins—A mechanistic perspective. AAPS J. 2013, 15, 897–900. [Google Scholar] [CrossRef] [PubMed]
  9. Steinman, R.M.; Brodie, S.E.; Cohn, Z.A. Membrane flow during pinocytosis. A stereologic analysis. J. Cell Biol. 1976, 68, 665–687. [Google Scholar] [CrossRef]
  10. Sallusto, F.; Cella, M.; Danieli, C.; Lanzavecchia, A. Dendritic cells use macropinocytosis and the mannose receptor to concentrate macromolecules in the major histocompatibility complex class II compartment: Downregulation by cytokines and bacterial products. J. Exp. Med. 1995, 182, 389–400. [Google Scholar] [CrossRef]
  11. Blank, F.; Rothen-Rutishauser, B.; Gehr, P. Dendritic cells and macrophages form a transepithelial network against foreign particulate antigens. Am. J. Respir. Cell Mol. Biol. 2007, 36, 669–677. [Google Scholar] [CrossRef]
  12. Varkhede, N.; Forrest, M.L. Understanding the monoclonal antibody disposition after subcutaneous administration using a minimal physiologically based pharmacokinetic model. J. Pharm. Pharm. Sci. 2018, 21, 130s. [Google Scholar] [CrossRef] [PubMed]
  13. Hu, S.; D’Argenio, D.Z. Predicting monoclonal antibody pharmacokinetics following subcutaneous administration via whole-body physiologically-based modeling. J. Pharmacokinet. Pharmacodyn. 2020, 47, 385–409. [Google Scholar] [CrossRef] [PubMed]
  14. Zheng, F.; Hou, P.; Corpstein, C.D.; Park, K.; Li, T. Multiscale pharmacokinetic modeling of systemic exposure of subcutaneously injected biotherapeutics. J. Control. Release 2021, 337, 407–416. [Google Scholar] [CrossRef] [PubMed]
  15. Gill, K.L.; Gardner, I.; Li, L.; Jamei, M. A bottom-up whole-body physiologically based pharmacokinetic model to mechanistically predict tissue distribution and the rate of subcutaneous absorption of therapeutic proteins. AAPS J. 2016, 18, 156–170. [Google Scholar] [CrossRef] [PubMed]
  16. Supersaxo, A.; Hein, W.R.; Steffen, H. Effect of molecular weight on the lymphatic absorption of water-soluble compounds following subcutaneous administration. Pharm. Res. 1990, 7, 167–169. [Google Scholar] [CrossRef]
  17. Zou, P.; Wang, F.; Wang, J.; Lu, Y.; Tran, D.; Seo, S.K. Impact of injection sites on clinical pharmacokinetics of subcutaneously administered peptides and proteins. J. Control. Release 2021, 336, 310–321. [Google Scholar] [CrossRef] [PubMed]
  18. Li, L.; Gardner, I.; Dostalek, M.; Jamei, M. Simulation of monoclonal antibody pharmacokinetics in humans using a minimal physiologically based model. AAPS J. 2014, 16, 1097–1109. [Google Scholar] [CrossRef] [PubMed]
  19. Rippe, B.; Haraldsson, B. Transport of macromolecules across microvascular walls: The two-pore theory. Physiol. Rev. 1994, 74, 163–219. [Google Scholar] [CrossRef] [PubMed]
  20. Crapo, J.D.; Barry, B.E.; Gehr, P.; Bachofen, M.; Weibel, E.R. Cell number and cell characteristics of the normal human lung. Am. Rev. Respir. Dis. 1982, 126, 332–337. [Google Scholar]
  21. Herring, N.; Paterson, D.J. Levick’s Introduction to Cardiovascular Physiology; CRC Press: Boca Raton, FL, USA, 2018. [Google Scholar]
  22. Nathanson, S.D.; Nelson, L.; Karvelis, K.C. Rates of flow of technetium 99m-labeled human serum albumin from peripheral injection sites to sentinel lymph nodes. Ann. Surg. Oncol. 1996, 3, 329–335. [Google Scholar] [CrossRef] [PubMed]
  23. Patel, N.; Clarke, J.F.; Salem, F.; Abdulla, T.; Martins, F.; Arora, S.; Tsakalozou, E.; Hodgkinson, A.; Arjmandi-Tash, O.; Cristea, S. Multi-phase multi-layer mechanistic dermal absorption (MPML MechDermA) model to predict local and systemic exposure of drug products applied on skin. CPT Pharmacomet. Syst. Pharmacol. 2022, 11, 1060–1084. [Google Scholar] [CrossRef]
  24. Valentin, J. Guide for the Practical Application of the ICRP Human Respiratory Tract Model: ICRP Supporting Guidance 3Approved by ICRP Committee 2 in October 2000. Ann. ICRP 2002, 32, 13. [Google Scholar] [CrossRef]
  25. University of Arkansas for Medical Science. Lymphatic Tables. 2023. Available online: https://medicine.uams.edu/neurobiology/education/medical-school-courses/human-structure-module/anatomy-tables/lymphatic-tables/. (accessed on 22 June 2023).
  26. Jafarnejad, M.; Woodruff, M.C.; Zawieja, D.C.; Carroll, M.C.; Moore, J., Jr. Modeling lymph flow and fluid exchange with blood vessels in lymph nodes. Lymphat. Res. Biol. 2015, 13, 234–247. [Google Scholar] [CrossRef] [PubMed]
  27. Savage, D.F.; Egea, P.F.; Robles-Colmenares, Y.; O’Connell, J.D., III; Stroud, R.M. Architecture and selectivity in aquaporins: 2.5 Å X-ray structure of aquaporin Z. PLoS Biol. 2003, 1, e72. [Google Scholar] [CrossRef] [PubMed]
  28. Mellor, R.; Stanton, A.; Azarbod, P.; Sherman, M.; Levick, J.; Mortimer, P. Enhanced cutaneous lymphatic network in the forearms of women with postmastectomy oedema. J. Vasc. Res. 2000, 37, 501–512. [Google Scholar] [CrossRef]
  29. Ying, M.; Pang, B. Three-dimensional ultrasound measurement of cervical lymph node volume. Br. J. Radiol. 2009, 82, 617–625. [Google Scholar] [CrossRef]
  30. Zuther, J.; Norton, S. Lymphedema Management; Thieme: Stuttgart, Germany, 2013; Volume 3. [Google Scholar]
  31. Camara, A.; Lavanant, A.C.; Abe, J.; Desforges, H.L.; Alexandre, Y.O.; Girardi, E.; Igamberdieva, Z.; Asano, K.; Tanaka, M.; Hehlgans, T. CD169+ macrophages in lymph node and spleen critically depend on dual RANK and LTbetaR signaling. Proc. Natl. Acad. Sci. USA 2022, 119, e2108540119. [Google Scholar] [CrossRef] [PubMed]
  32. Komohara, Y.; Ohnishi, K.; Takeya, M. Possible functions of CD 169-positive sinus macrophages in lymph nodes in anti-tumor immune responses. Cancer Sci. 2017, 108, 290–295. [Google Scholar] [CrossRef]
  33. Mellman, I.S.; Steinman, R.M.; Unkeless, J.C.; Cohn, Z.A. Selective iodination and polypeptide composition of pinocytic vesicles. J. Cell Biol. 1980, 86, 712–722. [Google Scholar] [CrossRef] [PubMed]
  34. Walk, E.L.; McLaughlin, S.; Coad, J.; Weed, S.A. Use of high frequency ultrasound to monitor cervical lymph node alterations in mice. PLoS ONE 2014, 9, e100185. [Google Scholar] [CrossRef] [PubMed]
  35. Smith, G.; Griffiths, B.; Mollison, D.; Mollison, P. Uptake of IgG after intramuscular and subcutaneous injection. Lancet 1972, 299, 1208–1212. [Google Scholar] [CrossRef]
  36. Li, Z.; Yu, X.; Li, Y.; Verma, A.; Chang, H.P.; Shah, D.K. A two-pore physiologically based pharmacokinetic model to predict subcutaneously administered different-size antibody/antibody fragments. AAPS J. 2021, 23, 62. [Google Scholar] [CrossRef] [PubMed]
  37. Hu, S.; Datta-Mannan, A.; D’Argenio, D.Z. Physiologically based modeling to predict monoclonal antibody pharmacokinetics in humans from in vitro physiochemical properties. mAbs 2022, 14, 2056944. [Google Scholar] [CrossRef] [PubMed]
  38. Fan, Y.-Y.; Farrokhi, V.; Caiazzo, T.; Wang, M.; O’Hara, D.M.; Neubert, H. Human FcRn tissue expression profile and half-life in PBMCs. Biomolecules 2019, 9, 373. [Google Scholar] [CrossRef] [PubMed]
  39. Zou, P. Predicting human bioavailability of subcutaneously administered fusion proteins and monoclonal antibodies using human intravenous clearance or antibody isoelectric point. AAPS J. 2023, 25, 31. [Google Scholar] [CrossRef] [PubMed]
  40. Schoch, A.; Kettenberger, H.; Mundigl, O.; Winter, G.; Engert, J.; Heinrich, J.; Emrich, T. Charge-mediated influence of the antibody variable domain on FcRn-dependent pharmacokinetics. Proc. Natl. Acad. Sci. USA 2015, 112, 5997–6002. [Google Scholar] [CrossRef] [PubMed]
  41. Liu, S.; Shah, D.K. Physiologically based pharmacokinetic modeling to characterize the effect of molecular charge on whole-body disposition of monoclonal antibodies. AAPS J. 2023, 25, 48–64. [Google Scholar] [CrossRef]
  42. Mellman, I.S.; Plutner, H.; Steinman, R.M.; Unkeless, J.C.; Cohn, Z.A. Internalization and degradation of macrophage Fc receptors during receptor-mediated phagocytosis. J. Cell Biol. 1983, 96, 887–895. [Google Scholar] [CrossRef]
Figure 1. Structure of the subcutaneous absorption model, coupled to a whole-body PBPK model in the Simcyp Simulator V22. Arrows represent the disposition of the simulated TP. APCs = antigen-presenting cells, CLcat = catabolic clearance, fcat = fraction of catabolic clearance, FR = fraction recycled, Kcat = catabolic clearance, KD = binding affinity to FcRn at pH 6.0, Krc = recycling rate, Kint = internalization rate, Kup = fluid-phase endocytosis uptake rate, Laff = afferent lymph flow, Leff = efferent lymph flow, PS = permeability surface area product of small and large pores, Qsc = blood flow to the SC site, σl = lymphatic reflection coefficient, σv = vascular reflection coefficient through small and large pores.
Figure 1. Structure of the subcutaneous absorption model, coupled to a whole-body PBPK model in the Simcyp Simulator V22. Arrows represent the disposition of the simulated TP. APCs = antigen-presenting cells, CLcat = catabolic clearance, fcat = fraction of catabolic clearance, FR = fraction recycled, Kcat = catabolic clearance, KD = binding affinity to FcRn at pH 6.0, Krc = recycling rate, Kint = internalization rate, Kup = fluid-phase endocytosis uptake rate, Laff = afferent lymph flow, Leff = efferent lymph flow, PS = permeability surface area product of small and large pores, Qsc = blood flow to the SC site, σl = lymphatic reflection coefficient, σv = vascular reflection coefficient through small and large pores.
Antibodies 13 00070 g001
Figure 2. Predicted vs. observed IgG plasma concentration (A) and removal rate from the arm (B). The red markers represent the observed data [15,35]. The solid lines, the dashed line, and the shaded area show the mean of each virtual trial and the mean and the confidence interval of the entire virtual population.
Figure 2. Predicted vs. observed IgG plasma concentration (A) and removal rate from the arm (B). The red markers represent the observed data [15,35]. The solid lines, the dashed line, and the shaded area show the mean of each virtual trial and the mean and the confidence interval of the entire virtual population.
Antibodies 13 00070 g002
Figure 3. Percentage of free FcRn in antigen-presenting cells (APCs) in peripheral lymph nodes draining the SC site. Dash dotted line = arm, dashed line = abdomen, solid line = back, dotted line = thigh.
Figure 3. Percentage of free FcRn in antigen-presenting cells (APCs) in peripheral lymph nodes draining the SC site. Dash dotted line = arm, dashed line = abdomen, solid line = back, dotted line = thigh.
Antibodies 13 00070 g003
Figure 4. Population sensitivity analysis on the fluid-phase endocytosis uptake (Kup) into APCs (A), the binding affinity to FcRn (B), and the dose (C) after injection of IgG into the arm. The black line shows the mean and the gray shaded area denotes the 95% confidence interval of a virtual population containing 100 individuals (50% female) aged 20 to 50 years.
Figure 4. Population sensitivity analysis on the fluid-phase endocytosis uptake (Kup) into APCs (A), the binding affinity to FcRn (B), and the dose (C) after injection of IgG into the arm. The black line shows the mean and the gray shaded area denotes the 95% confidence interval of a virtual population containing 100 individuals (50% female) aged 20 to 50 years.
Antibodies 13 00070 g004
Table 2. Site-specific IgG predictions against observed data from 47 IgG-based antibodies [17].
Table 2. Site-specific IgG predictions against observed data from 47 IgG-based antibodies [17].
Cmax Ratiotmax RatioAUC Ratio
ObservedPredictedObservedPredictedObservedPredicted
Arm/Abdomen0.920.961.260.980.980.95
Thigh/Abdomen1.141.130.970.771.131.06
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Stader, F.; Liu, C.; Derbalah, A.; Momiji, H.; Pan, X.; Gardner, I.; Jamei, M.; Sepp, A. A Physiologically Based Pharmacokinetic Model Relates the Subcutaneous Bioavailability of Monoclonal Antibodies to the Saturation of FcRn-Mediated Recycling in Injection-Site-Draining Lymph Nodes. Antibodies 2024, 13, 70. https://doi.org/10.3390/antib13030070

AMA Style

Stader F, Liu C, Derbalah A, Momiji H, Pan X, Gardner I, Jamei M, Sepp A. A Physiologically Based Pharmacokinetic Model Relates the Subcutaneous Bioavailability of Monoclonal Antibodies to the Saturation of FcRn-Mediated Recycling in Injection-Site-Draining Lymph Nodes. Antibodies. 2024; 13(3):70. https://doi.org/10.3390/antib13030070

Chicago/Turabian Style

Stader, Felix, Cong Liu, Abdallah Derbalah, Hiroshi Momiji, Xian Pan, Iain Gardner, Masoud Jamei, and Armin Sepp. 2024. "A Physiologically Based Pharmacokinetic Model Relates the Subcutaneous Bioavailability of Monoclonal Antibodies to the Saturation of FcRn-Mediated Recycling in Injection-Site-Draining Lymph Nodes" Antibodies 13, no. 3: 70. https://doi.org/10.3390/antib13030070

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop