Predicting monoclonal antibody pharmacokinetics following subcutaneous administration via whole-body physiologically-based modeling

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

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.

Inter-Antibody Variability in the Clinical Pharmacokinetics of Monoclonal Antibodies Characterized Using Population Physiologically Based Pharmacokinetic Modeling

The objective of this work was to develop a population physiologically based pharmacokinetic (popPBPK) model to characterize the variability in the clinical PK of monoclonal antibodies (mAbs) following intravenous (IV) and subcutaneous (SC) administration. An extensive literature search was conducted and clinical PK data for FDA-approved as well as non-approved mAbs were collected. Training and validation datasets of 44 and 9 mAbs exhibiting linear pharmacokinetics were used for model development. The variability in antibody PK was captured by accounting for different rate constants of pinocytosis (CLup) and intracellular degradation (kdeg) for different mAbs. Typical values for CLup and kdeg and their respective inter-antibody variabilities (ωClup, ωKdeg) were estimated to be 0.32 L/h/L and 26.1 h-1 (73% and 46%). Varied absorption profiles following SC dosing were characterized by incorporating inter-antibody variability in local degradation (kSC) and rate of lymphatic uptake (S_Lu) of mAbs. Estimates for typical kSC and S_Lu values, and ωKsc,ωS_Lu, were found to be 0.0015 h-1 and 0.54 (193%, and 49%). FDA-approved mAbs showed less local degradation (0.0014 h-1 vs. 0.0038 h-1) compared with other clinically tested mAbs, whereas no substantial differences in physiological processes involved in disposition were observed. To evaluate the generalizability of estimated PK parameters and model validation, the final popPBPK model was used to simulate the range of expected PK for mAbs following SC administration of nine different mAbs that were not used for model-building purposes. The predicted PK of all nine mAbs was within the expected range specified a priori. Thus, the popPBPK model presented here may serve as a tool to predict the clinical PK of mAbs with linear disposition before administering them to humans. The model may also support preclinical-to-clinical translation and 'first-in-human' dose determination for mAbs.

Translational two-pore PBPK model to characterize whole-body disposition of different-size endogenous and exogenous proteins

Two-pore physiologically based pharmacokinetic (PBPK) modeling has demonstrated its potential in describing the pharmacokinetics (PK) of different-size proteins. However, all existing two-pore models lack either diverse proteins for validation or interspecies extrapolation. To fill the gap, here we have developed and optimized a translational two-pore PBPK model that can characterize plasma and tissue disposition of different-size proteins in mice, rats, monkeys, and humans. Datasets used for model development include more than 15 types of proteins: IgG (150 kDa), F(ab)2 (100 kDa), minibody (80 kDa), Fc-containing proteins (205, 200, 110, 105, 92, 84, 81, 65, or 60 kDa), albumin conjugate (85.7 kDa), albumin (67 kDa), Fab (50 kDa), diabody (50 kDa), scFv (27 kDa), dAb2 (23.5 kDa), proteins with an albumin-binding domain (26, 23.5, 22, 16, 14, or 13 kDa), nanobody (13 kDa), and other proteins (110, 65, or 60 kDa). The PBPK model incorporates: (i) molecular weight (MW)-dependent extravasation through large and small pores via diffusion and filtration, (ii) MW-dependent renal filtration, (iii) endosomal FcRn-mediated protection from catabolism for IgG and albumin-related modalities, and (iv) competition for FcRn binding from endogenous IgG and albumin. The finalized model can well characterize PK of most of these proteins, with area under the curve predicted within two-fold error. The model also provides insights into contribution of renal filtration and lysosomal degradation towards total elimination of proteins, and contribution of paracellular convection/diffusion and transcytosis towards extravasation. The PBPK model presented here represents a cross-modality, cross-species platform that can be used for development of novel biologics.

Reducing Immunogenicity by Design: Approaches to Minimize Immunogenicity of Monoclonal Antibodies

Monoclonal antibodies (mAbs) have transformed therapeutic strategies for various diseases. Their high specificity to target antigens makes them ideal therapeutic agents for certain diseases. However, a challenge to their application in clinical practice is their potential risk to induce unwanted immune response, termed immunogenicity. This challenge drives the continued efforts to deimmunize these protein therapeutics while maintaining their pharmacokinetic properties and therapeutic efficacy. Because mAbs hold a central position in therapeutic strategies against an array of diseases, the importance of conducting comprehensive immunogenicity risk assessment during the drug development process cannot be overstated. Such assessment necessitates the employment of in silico, in vitro, and in vivo strategies to evaluate the immunogenicity risk of mAbs. Understanding the intricacies of the mechanisms that drive mAb immunogenicity is crucial to improving their therapeutic efficacy and safety and developing the most effective strategies to determine and mitigate their immunogenic risk. This review highlights recent advances in immunogenicity prediction strategies, with a focus on protein engineering strategies used throughout development to reduce immunogenicity.

A compartment model for subcutaneous injection of monoclonal antibodies

Despite the growing popularity of subcutaneous (SC) administration for monoclonal antibodies (mAbs), there remains a limited understanding of the significance of mAb transport rate constants within the interstitial space and the lymphatic system on their pharmacokinetics. To bridge this knowledge gap, we introduce a compartmental model for subcutaneously administered mAbs. Our model differentiates FcRn-expressing cells across various sites, and the model predictions agree with experimental data from both human and rat studies. Our findings indicate that the time to reach the maximum mAb concentration in the plasma, denoted by Tmax, displays a weak positive correlation with mAb half-life and a negligible correlation with bioavailability. In contrast, the half-life of mAbs exhibits a strong positive correlation with bioavailability. Moreover, the rate of mAb transport from lymph to plasma significantly affects the mAb half-life. Increasing the transport rates of mAbs from the injection site to the lymph or from lymph to plasma enhances bioavailability. These insights, combined with our compartmental model, contribute to a deeper understanding of the pharmacokinetics of subcutaneously administered mAbs.

Predicting the clinical subcutaneous absorption rate constant of monoclonal antibodies using only the primary sequence: a machine learning approach

Subcutaneous injections are an increasingly prevalent route of administration for delivering biological therapies including monoclonal antibodies (mAbs). Compared with intravenous delivery, subcutaneous injections reduce administration costs, shorten the administration time, and are strongly preferred from a patient experience point of view. An understanding of the absorption process of a mAb from the injection site to the systemic circulation is critical to the process of subcutaneous mAb formulation development. In this study, we built a model to predict the absorption rate constant (ka), which denotes how fast a mAb is absorbed from the site of administration. Once trained, our model (enabled by the XGBoost algorithm in machine learning) can predict the ka of a mAb following a subcutaneous injection using in silico molecular properties alone (generated from the primary sequence). Our model does not need clinically observed plasma concentration-time data; this is a novel capability not previously achieved in predictive pharmacokinetic models. The model also showed improved performance when benchmarked against a recently reported mechanistic model that relied on clinical data to predict subcutaneous absorption of mAbs. We further interpreted the model to understand which molecular properties affect the absorption rate and showed that our findings are consistent with previous studies evaluating subcutaneous absorption through direct experimentation. Taken altogether, this study reports the development, validation, benchmarking, and interpretation of a model that can predict the clinical ka of a mAb using its primary sequence as the only input.

Physiologically Based Pharmacokinetic Modeling to Characterize the Effect of Molecular Charge on Whole-Body Disposition of Monoclonal Antibodies

Motivated by a series of work demonstrating the effect of molecular charge on antibody pharmacokinetics (PK), physiological-based pharmacokinetic (PBPK) models are emerging that relate in silico calculated charge or in vitro measures of polyspecificity to antibody PK parameters. However, only plasma data has been used for model development in these studies, leading to unvalidated assumptions. Here, we present an extended platform PBPK model for antibodies that incorporate charge-dependent endothelial cell pinocytosis rate and nonspecific off-target binding in the interstitial space and on circulating blood cells, to simultaneously characterize whole-body disposition of three antibody charge variants. Predictive potential of various charge metrics was also explored, and the difference between positive charge patches and negative charge patches (i.e., PPC-PNC) was used as the charge parameter to establish quantitative relationships with nonspecific binding affinities and endothelial cell uptake rate. Whole-body disposition of these charge variants was captured well by the model, with less than 2-fold predictive error in area under the curve of most plasma and tissue PK data. The model also predicted that with greater positive charge, nonspecific binding was more substantial, and pinocytosis rate increased especially in brain, heart, kidney, liver, lung, and spleen, but remained unchanged in adipose, bone, muscle, and skin. The presented PBPK model contributes to our understanding of the mechanisms governing the disposition of charged antibodies and can be used as a platform to guide charge engineering based on desired plasma and tissue exposures.

Intravital Microscopy Reveals Unforeseen Biodistribution Within the Liver and Kidney Mechanistically Connected to the Clearance of a Bifunctional Antibody

Bifunctional antibody (BfAb) therapeutics offer the potential for novel functionalities beyond those of the individual monospecific entities. However, combining these entities into a single molecule can have unpredictable effects, including changes in pharmacokinetics that limit the compound's therapeutic profile. A better understanding of how molecular modifications affect in vivo tissue interactions could help inform BfAb design. The present studies were predicated on the observation that a BfAb designed to have minimal off-target interactions cleared from the circulation twice as fast as the monoclonal antibody (mAb) from which it was derived. The present study leverages the spatial and temporal resolution of intravital microscopy (IVM) to identify cellular interactions that may explain the different pharmacokinetics of the two compounds. Disposition studies of mice demonstrated that radiolabeled compounds distributed similarly over the first 24 hours, except that BfAb accumulated approximately two- to -three times more than mAb in the liver. IVM studies of mice demonstrated that both distributed to endosomes of liver endothelia but with different kinetics. Whereas mAb accumulated rapidly within the first hour of administration, BfAb accumulated only modestly during the first hour but continued to accumulate over 24 hours, ultimately reaching levels similar to those of the mAb. Although neither compound was freely filtered by the mouse or rat kidney, BfAb, but not mAb, was found to accumulate over 24 hours in endosomes of proximal tubule cells. These studies demonstrate how IVM can be used as a tool in drug design, revealing unpredicted cellular interactions that are undetectable by conventional analyses. SIGNIFICANCE STATEMENT: Bifunctional antibodies offer novel therapeutic functionalities beyond those of the individual monospecific entities. However, combining these entities into a single molecule can have unpredictable effects, including undesirable changes in pharmacokinetics. Studies of the dynamic distribution of a bifunctional antibody and its parent monoclonal antibody presented here demonstrate how intravital microscopy can expand our understanding of the in vivo disposition of therapeutics, detecting off-target interactions that could not be detected by conventional pharmacokinetics approaches or predicted by conventional physicochemical analyses.

Utility of physiologically based pharmacokinetic modeling to predict inter-antibody variability in monoclonal antibody pharmacokinetics in mice

In this investigation, we tested the hypothesis that a physiologically based pharmacokinetic (PBPK) model incorporating measured in vitro metrics of off-target binding can largely explain the inter-antibody variability in monoclonal antibody (mAb) pharmacokinetics (PK). A diverse panel of 83 mAbs was evaluated for PK in wild-type mice and subjected to 10 in vitro assays to measure major physiochemical attributes. After excluding for target-mediated elimination and immunogenicity, 56 of the remaining mAbs with an eight-fold variability in the area under the curve (A U C 0 - 672 h : 1.74 × 106 -1.38 × 107 ng∙h/mL) and 10-fold difference in clearance (2.55-26.4 mL/day/kg) formed the training set for this investigation. Using a PBPK framework, mAb-dependent coefficients F1 and F2 modulating pinocytosis rate and convective transport, respectively, were estimated for each mAb with mostly good precision (coefficient of variation (CV%) <30%). F1 was estimated to be the mean and standard deviation of 0.961 ± 0.593, and F2 was estimated to be 2.13 ± 2.62. Using principal component analysis to correlate the regressed values of F1/F2 versus the multidimensional dataset composed of our panel of in vitro assays, we found that heparin chromatography retention time emerged as the predictive covariate to the mAb-specific F1, whereas F2 variability cannot be well explained by these assays. A sigmoidal relationship between F1 and the identified covariate was incorporated within the PBPK framework. A sensitivity analysis suggested plasma concentrations to be most sensitive to F1 when F1 > 1. The predictive utility of the developed PBPK model was evaluated against a separate panel of 14 mAbs biased toward high clearance, among which area under the curve of PK data of 12 mAbs was predicted within 2.5-fold error, and the positive and negative predictive values for clearance prediction were 85% and 100%, respectively. MAb heparin chromatography assay output allowed a priori identification of mAb candidates with unfavorable PK.

Monoclonal Antibody Pharmacokinetics in Cynomolgus Monkeys Following Subcutaneous Administration: Physiologically Based Model Predictions from Physiochemical Properties

An integrated physiologically based modeling framework is presented for predicting pharmacokinetics and bioavailability of subcutaneously administered monoclonal antibodies in cynomolgus monkeys, based on in silico structure-derived metrics characterizing antibody size, overall charge, local charge, and hydrophobicity. The model accounts for antibody-specific differences in pinocytosis, transcapillary transport, local lymphatic uptake, and pre-systemic degradation at the subcutaneous injection site and reliably predicts the pharmacokinetics of five different wild-type mAbs and their Fc variants following intravenous and subcutaneous administration. Significant associations were found between subcutaneous injection site degradation rate and the antibody's local positive charge of its complementarity-determining region (R = 0.56, p = 0.0012), antibody pinocytosis rate and its overall positive charge (R = 0.59, p = 0.00063), and antibody paracellular transport and its overall charge together with hydrophobicity (R = 0.63, p = 0.00096). Based on these results, population simulations were performed to predict the relationship between bioavailability and antibody local positive charge. In addition, model simulations were conducted to calculate the relative contribution of absorption pathways (lymphatic and blood), pre-systemic degradation pathways (interstitial and lysosomal), and the influence of injection site lymph flow on antibody bioavailability and pharmacokinetics. The proposed physiologically based modeling framework integrates fundamental mechanisms governing antibody subcutaneous absorption and disposition, with structured-based physiochemical properties, to predict antibody bioavailability and pharmacokinetics in vivo.

Application of physiologically-based pharmacokinetic models for therapeutic proteins and other novel modalities

The past two decades have seen diversification of drug development pipelines and approvals from traditional small molecule therapies to alternative modalities including monoclonal antibodies, engineered proteins, antibody drug conjugates (ADCs), oligonucleotides and gene therapies. At the same time, physiologically-based pharmacokinetic (PBPK) models for small molecules have seen increased industry and regulatory acceptance.This review focusses on the current status of the application of PBPK models to these newer modalities and give a perspective on the successes, challenges and future directions of this field.There is greatest experience in the development of PBPK models for therapeutic proteins, and PBPK models for ADCs benefit from prior experience for both therapeutic proteins and small molecules. For other modalities, the application of PBPK models is in its infancy.Challenges are discussed and a common theme is lack of availability of physiological and experimental data to characterise systems and drug parameters to enable a priori prediction of pharmacokinetics. Furthermore, sufficient clinical data are required to build confidence in developed models.The PBPK modelling approach provides a quantitative framework for integrating knowledge and data from multiple sources and can be built on as more data becomes available.

Mathematical Models to Characterize the Absorption, Distribution, Metabolism, and Excretion of Protein Therapeutics

Therapeutic proteins (TPs) have ranked among the most important and fastest-growing classes of drugs in the clinic, yet the development of successful TPs is often limited by unsatisfactory efficacy. Understanding pharmacokinetic (PK) characteristics of TPs is key to achieving sufficient and prolonged exposure at the site of action, which is a prerequisite for eliciting desired pharmacological effects. PK modeling represents a powerful tool to investigate factors governing in vivo disposition of TPs. In this mini-review, we discuss many state-of-the-art models that recapitulate critical processes in each of the absorption, distribution, metabolism/catabolism, and excretion pathways of TPs, which can be integrated into the physiologically-based pharmacokinetic framework. Additionally, we provide our perspectives on current opportunities and challenges for evolving the PK models to accelerate the discovery and development of safe and efficacious TPs. SIGNIFICANCE STATEMENT: This minireview provides an overview of mechanistic pharmacokinetic (PK) models developed to characterize absorption, distribution, metabolism, and elimination (ADME) properties of therapeutic proteins (TPs), which can support model-informed discovery and development of TPs. As the next-generation of TPs with diverse physicochemical properties and mechanism-of-action are being developed rapidly, there is an urgent need to better understand the determinants for the ADME of TPs and evolve existing platform PK models to facilitate successful bench-to-bedside translation of these promising drug molecules.

Recent Advances in Translational Pharmacokinetics and Pharmacodynamics Prediction of Therapeutic Antibodies Using Modeling and Simulation

Therapeutic monoclonal antibodies (mAbs) have been a promising therapeutic approach for several diseases and a wide variety of mAbs are being evaluated in clinical trials. To accelerate clinical development and improve the probability of success, pharmacokinetics and pharmacodynamics (PKPD) in humans must be predicted before clinical trials can begin. Traditionally, empirical-approach-based PKPD prediction has been applied for a long time. Recently, modeling and simulation (M&S) methods have also become valuable for quantitatively predicting PKPD in humans. Although several models (e.g., the compartment model, Michaelis–Menten model, target-mediated drug disposition model, and physiologically based pharmacokinetic model) have been established and used to predict the PKPD of mAbs in humans, more complex mechanistic models, such as the quantitative systemics pharmacology model, have been recently developed. This review summarizes the recent advances and future direction of M&S-based approaches to the quantitative prediction of human PKPD for mAbs.

Opportunities and challenges of physiologically based pharmacokinetic modeling in drug delivery

Physiologically based pharmacokinetic (PBPK) modeling is an important in silico tool to bridge drug properties and in vivo PK behaviors during drug development. Over the recent decade, the PBPK method has been largely applied to drug delivery systems (DDS), including oral, inhaled, transdermal, ophthalmic, and complex injectable products. The related therapeutic agents have included small-molecule drugs, therapeutic proteins, nucleic acids, and even cells. Simulation results have provided important insights into PK behaviors of new dosage forms, which strongly support drug regulation. In this review, we comprehensively summarize recent progress in PBPK applications in drug delivery, which shows large opportunities for facilitating drug development. In addition, we discuss the challenges of applying this methodology from a practical viewpoint.

Minimal Physiologically-based Pharmacokinetic Model to Investigate the Effect of Charge on the Pharmacokinetics of Humanized anti-HCV-E2 IgG Antibodies in Sprague-Dawley Rats

PURPOSE

To develop a minimal physiologically-based pharmacokinetic (mPBPK) model in quantifying the relationships between the charge and pharmacokinetics (PK) of therapeutic monoclonal IgG antibody (TMAb).

METHODS

PK data used in this study were native IgG and five humanized anti-HCVE2-IgG antibodies in rats. Different models that related the effect of charge on interstitial distribution, transcapillary transport, and cellular uptake for FcRn-mediated metabolism were tested. External validation was conducted to assess if the charge-parameter relationships derived from rats could be used to predict the PK of TMAbs in mice. The final mPBPK model was used to construct the relationships between the FcRn binding and charge on the PK of TMAbs.

RESULTS

Increasing the isoelectric point (pI) of IgG was associated with higher interstitial space distribution and cellular uptake. The transcapillary transport of IgG from plasma to interstitial space remains constant with pI values below 7.96 and then increased linearly with pI. The model-based simulation results suggested that improving the FcRn binding affinity can overcome the problems of low plasma/interstitial space exposures associated with TMAbs with higher pI values by reducing the FcRn-mediated metabolism and hence increasing drug exposure in the interstitial space that has close contact with many solid tumors.

CONCLUSIONS

The final mPBPK model was developed and used to construct complex quantitative relationships between the pI/FcRn binding affinity and PK of TMAbs and such relationships are useful to select the discovery of a "sweet spot" of designing future generation of TMAbs with optimal PK properties to achieve desirable plasma and tissue drug exposures.

Physiologically Based Modeling to Predict Monoclonal Antibody Pharmacokinetics in Humans from in vitro Physiochemical Properties

A model-based framework is presented to predict monoclonal antibody (mAb) pharmacokinetics (PK) in humans based on in vitro measures of antibody physiochemical properties. A physiologically based pharmacokinetic (PBPK) model is used to explore the predictive potential of 14 in vitro assays designed to measure various antibody physiochemical properties, including nonspecific cell-surface interactions, FcRn binding, thermal stability, hydrophobicity, and self-association. Based on the mean plasma PK time course data of 22 mAbs from humans reported in the literature, we found a significant positive correlation (R = 0.64, p = .0013) between the model parameter representing antibody-specific vascular to endothelial clearance and heparin relative retention time, an in vitro measure of nonspecific binding. We also found that antibody-specific differences in paracellular transport due to convection and diffusion could be partially explained by antibody heparin relative retention time (R = 0.52, p = .012). Other physiochemical properties, including antibody thermal stability, hydrophobicity, cross-interaction and self-association, in and of themselves were not predictive of model-based transport parameters. In contrast to other studies that have reported empirically derived expressions relating in vitro measures of antibody physiochemical properties directly to antibody clearance, the proposed PBPK model-based approach for predicting mAb PK incorporates fundamental mechanisms governing antibody transport and processing, informed by in vitro measures of antibody physiochemical properties, and can be expanded to include more descriptive representations of each of the antibody processing subsystems, as well as other antibody-specific information.

Non-human primates in the PKPD evaluation of biologics: Needs and options to reduce, refine, and replace. A BioSafe White Paper

Monoclonal antibodies (mAbs) deliver great benefits to patients with chronic and/or severe diseases thanks to their strong specificity to the therapeutic target. As a result of this specificity, non-human primates (NHP) are often the only preclinical species in which therapeutic antibodies cross-react with the target. Here, we highlight the value and limitations that NHP studies bring to the design of safe and efficient early clinical trials. Indeed, data generated in NHPs are integrated with in vitro information to predict the concentration/effect relationship in human, and therefore the doses to be tested in first-in-human trials. The similarities and differences in the systems defining the pharmacokinetics and pharmacodynamics (PKPD) of mAbs in NHP and human define the nature and the potential of the preclinical investigations performed in NHPs. Examples have been collated where the use of NHP was either pivotal to the design of the first-in-human trial or, inversely, led to the termination of a project prior to clinical development. The potential impact of immunogenicity on the results generated in NHPs is discussed. Strategies to optimize the use of NHPs for PKPD purposes include the addition of PD endpoints in safety assessment studies and the potential re-use of NHPs after non-terminal studies or cassette dosing several therapeutic agents of interest. Efforts are also made to reduce the use of NHPs in the industry through the use of in vitro systems, alternative in vivo models, and in silico approaches. In the case of prediction of ocular PK, the body of evidence gathered over the last two decades renders the use of NHPs obsolete. Expert perspectives, advantages, and pitfalls with these alternative approaches are shared in this review.

Effect of variable domain charge on in vitro and in vivo disposition of monoclonal antibodies

A growing body of evidence supports the important role of molecular charge on antibody pharmacokinetics (PK), yet a quantitative description of the effect of charge on systemic and tissue disposition of antibodies is still lacking. Consequently, we have systematically engineered complementarity-determining regions (CDRs) of trastuzumab to create a series of variants with an isoelectric point (pI) range of 6.3–8.9 and a variable region (Fv) charge range of −8.9 to +10.9 (at pH 5.5), and have investigated in vitro and in vivo disposition of these molecules. These monoclonal antibodies (mAbs) exhibited incrementally enhanced binding to cell surfaces and cellular uptake with increased positive charge in antigen-negative cells. After single intravenous dosing in mice, a bell-shaped relationship between systemic exposure and Fv charge was observed, with both extended negative and positive charge patches leading to more rapid nonspecific clearance. Whole-body PK experiments revealed that, although overall exposures of most variants in the tissues were very similar, positive charge of mAbs led to significantly enhanced tissue:plasma concentration ratios for most tissues. In well-perfused organs such as liver, spleen, and kidney, the positive charge variants show superior accumulation. In tissues with continuous capillaries such as fat, muscle, skin, and bone, plasma concentrations governed tissue exposures. The in vitro and in vivo disposition data presented here facilitate better understanding of the impact of charge modifications on antibody PK, and suggest that alteration in the charge may help to improve tissue:plasma concentration ratios for mAbs in certain tissues. The data presented here also paves the way for the development of physiologically based pharmacokinetic models of mAbs that incorporate charge variations.

A Two-Pore Physiologically Based Pharmacokinetic Model to Predict Subcutaneously Administered Different-Size Antibody/Antibody Fragments

Quantitative modeling of the subcutaneous absorption processes of protein therapeutics is challenging. Here we have proposed a "two-pore" PBPK model that is able to simultaneously characterize plasma PK of different-size protein therapeutics in mice. The skin compartment is evolved to mechanistically account for the absorption pathways through lymph and blood capillaries, as well as local degradation at the SC injection site. The model is developed using in-house plasma PK data generated following subcutaneous administration of 6 different-size protein therapeutics (13-150 kDa) in mice. The model was able to capture plasma PK of all molecules following intravenous and subcutaneous administration relatively well. From the observed plasma PK profiles, as well as from the model simulation result, several important PK descriptors were found to be dependent on protein size for FcRn nonbinding molecules. A positive correlation was found between T_max and protein size. A "U" shape relationship was found between C_max and protein size. Negative correlations were observed between bioavailability (F) and local degradation rate (k_deg,SC), and F and protein size. Pathway analysis of the model was conducted for the subcutaneous absorption process, and continuous relationships were established between the percentage of absorption through lymphatic and vascular pathways and protein size. This PBPK model could serve as a platform for the development of different-size protein therapeutics and will be scaled up to humans for translational studies in the future.