A predictive model of therapeutic monoclonal antibody dynamics and regulation by the neonatal Fc receptor (FcRn)

Comparison of monoclonal antibody disposition predictions using different physiologically based pharmacokinetic modelling platforms

Physiologically based pharmacokinetic (PBPK) models can be used to leverage physiological and in vitro data to predict monoclonal antibody (mAb) concentrations in serum and tissues. However, it is currently not known how consistent predictions of mAb disposition are across PBPK modelling platforms. In this work PBPK simulations of IgG, adalimumab and infliximab were compared between three platforms (Simcyp, PK-Sim, and GastroPlus). Accuracy of predicted serum and tissue concentrations was assessed using observed data collected from the literature. Physiological and mAb related input parameters were also compared and sensitivity analyses were carried out to evaluate model behavior when input values were altered. Differences in serum kinetics of IgG between platforms were minimal for a dose of 1 mg/kg, but became more noticeable at higher dosages (> 100 mg/kg) and when reference (healthy) physiological input values were altered. Predicted serum concentrations of both adalimumab and infliximab were comparable across platforms, but were noticeably higher than observed values. Tissue concentrations differed remarkably between the platforms, both for total- and interstitial fluid (ISF) concentrations. The accuracy of total tissue concentrations was within a three-fold of observed values for all tissues, except for brain tissue concentrations, which were overpredicted. Predictions of tissue ISF concentrations were less accurate and were best captured by GastroPlus. Overall, these simulations show that the different PBPK platforms generally predict similar mAb serum concentrations, but variable tissue concentrations. Caution is therefore warranted when PBPK models are used to simulate effect site tissue concentrations of mAbs without data to verify the predictions.

Emerging trends and therapeutic applications of monoclonal antibodies

Monoclonal antibodies (mAbs) are being used to prevent, detect, and treat a broad spectrum of malignancies and infectious and autoimmune diseases. Over the past few years, the market for mAbs has grown exponentially. They have become a significant part of many pharmaceutical product lines, and more than 250 therapeutic mAbs are undergoing clinical trials. Ever since the advent of hybridoma technology, antibody-based therapeutics were realized using murine antibodies which further progressed into humanized and fully human antibodies, reducing the risk of immunogenicity. Some of the benefits of using mAbs over conventional drugs include a drastic reduction in the chances of adverse reactions, interactions between drugs, and targeting specific proteins. While antibodies are very efficient, their higher production costs impede the process of commercialization. However, their cost factor has been improved by developing biosimilar antibodies, which are affordable versions of therapeutic antibodies. Along with biosimilars, innovations in antibody engineering have helped to design bio-better antibodies with improved efficacy than the conventional ones. These novel mAb-based therapeutics are set to revolutionize existing drug therapies targeting a wide spectrum of diseases, thereby meeting several unmet medical needs. In the future, mAbs generated by applying next-generation sequencing (NGS) are expected to become a powerful tool in clinical therapeutics. This article describes the methods of mAb production, pre-clinical and clinical development of mAbs, approved indications targeted by mAbs, and novel developments in the field of mAb research.

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.

A minimal physiologically based pharmacokinetic model to study the combined effect of antibody size, charge, and binding affinity to FcRn/antigen on antibody pharmacokinetics

Protein therapeutics have revolutionized the treatment of a wide range of diseases. While they have distinct physicochemical characteristics that influence their absorption, distribution, metabolism, and excretion (ADME) properties, the relationship between the physicochemical properties and PK is still largely unknown. In this work we present a minimal physiologically-based pharmacokinetic (mPBPK) model that incorporates a multivariate quantitative relation between a therapeutic's physicochemical parameters and its corresponding ADME properties. The model's compound-specific input includes molecular weight, molecular size (Stoke's radius), molecular charge, binding affinity to FcRn, and specific antigen affinity. Through derived and fitted empirical relationships, the model demonstrates the effect of these compound-specific properties on antibody disposition in both plasma and peripheral tissues using observed PK data in mice and humans. The mPBPK model applies the two-pore hypothesis to predict size-based clearance and exposure of full-length antibodies (150 kDa) and antibody fragments (50-100 kDa) within a onefold error. We quantitatively relate antibody charge and PK parameters like uptake rate, non-specific binding affinity, and volume of distribution to capture the relatively faster clearance of positively charged mAb as compared to negatively charged mAb. The model predicts the terminal plasma clearance of slightly positively and negatively charged antibody in humans within a onefold error. The mPBPK model presented in this work can be used to predict the target-mediated disposition of a drug when compound-specific and target-specific properties are known. To our knowledge, a combined effect of antibody weight, size, charge, FcRn, and antigen has not been incorporated and studied in a single mPBPK model previously. By conclusively incorporating and relating a multitude of protein's physicochemical properties to observed PK, our mPBPK model aims to contribute as a platform approach in the early stages of drug development where many of these properties can be optimized to improve a molecule's PK and ultimately its efficacy.

Kinetic modeling of the plasma pharmacokinetic profiles of ADAMTS13 fragment and its Fc-fusion counterpart in mice

Introduction: Fusion of the fragment crystallizable (Fc) to protein therapeutics is commonly used to extend the circulation time by enhancing neonatal Fc-receptor (FcRn)-mediated endosomal recycling and slowing renal clearance. This study applied kinetic modeling to gain insights into the cellular processing contributing to the observed pharmacokinetic (PK) differences between the novel recombinant ADAMTS13 fragment (MDTCS) and its Fc-fusion protein (MDTCS-Fc). Methods: For MDTCS and MDTCS-Fc, their plasma PK profiles were obtained at two dose levels following intravenous administration of the respective proteins to mice. The plasma PK profiles of MDTCS were fitted to a kinetic model with three unknown protein-dependent parameters representing the fraction recycled (FR) and the rate constants for endocytosis (kup, for the uptake into the endosomes) and for the transfer from the plasma to the interstitial fluid (kpi). For MDTCS-Fc, the model was modified to include an additional parameter for binding to FcRn. Parameter optimization was done using the Cluster Gauss-Newton Method (CGNM), an algorithm that identifies multiple sets of approximate solutions ("accepted" parameter sets) to nonlinear least-squares problems. Results: As expected, the kinetic modeling results yielded the FR of MDTCS-Fc to be 2.8-fold greater than that of MDTCS (0.8497 and 0.3061, respectively). In addition, MDTCS-Fc was predicted to undergo endocytosis (the uptake into the endosomes) at a slower rate than MDTCS. Sensitivity analyses identified the association rate constant (kon) between MDTCS-Fc and FcRn as a potentially important factor influencing the plasma half-life in vivo. Discussion: Our analyses suggested that Fc fusion to MDTCS leads to changes in not only the FR but also the uptake into the endosomes, impacting the systemic plasma PK profiles. These findings may be used to develop recombinant protein therapeutics with extended circulation time.

Non-specific irreversible 89Zr-mAb uptake in tumours: evidence from biopsy-proven target-negative tumours using 89Zr-immuno-PET

BACKGROUND

Distribution of mAbs into tumour tissue may occur via different processes contributing differently to the 89Zr-mAb uptake on PET. Target-specific binding in tumours is of main interest; however, non-specific irreversible uptake may also be present, which influences quantification. The aim was to investigate the presence of non-specific irreversible uptake in tumour tissue using Patlak linearization on 89Zr-immuno-PET data of biopsy-proven target-negative tumours. Data of two studies, including target status obtained from biopsies, were retrospectively analysed, and Patlak linearization provided the net rate of irreversible uptake (Ki).

RESULTS

Two tumours were classified as CD20-negative and two as CD20-positive. Four tumours were classified as CEA-negative and nine as CEA-positive. Ki values of CD20-negative (0.43 µL/g/h and 0.92 µL/g/h) and CEA-negative tumours (mdn = 1.97 µL/g/h, interquartile range (IQR) = 1.50-2.39) were higher than zero. Median Ki values of target-negative tumours were lower than CD20-positive (1.87 µL/g/h and 1.90 µL/g/h) and CEA-positive tumours (mdn = 2.77 µL/g/h, IQR = 2.11-3.65).

CONCLUSION

Biopsy-proven target-negative tumours showed irreversible uptake of 89Zr-mAbs measured in vivo using 89Zr-immuno-PET data, which suggests the presence of non-specific irreversible uptake in tumours. Consequently, for 89Zr-immuno-PET, even if the target is absent, a tumour-to-plasma ratio always increases over time.

Clinical validation of translational antibody PBPK model using tissue distribution data generated with 89Zr-immuno-PET imaging

The main objective of this manuscript was to validate the ability of the monoclonal antibody physiologically-based pharmacokinetic (PBPK) model to predict tissue concentrations of antibodies in the human. To accomplish this goal, preclinical and clinical tissue distribution and positron emission tomography imaging data generated using zirconium-89 (89Zr) labeled antibodies were obtained from the literature. First, our previously published translational PBPK model for antibodies was expanded to describe the whole-body biodistribution of 89Zr labeled antibody and the free 89Zr, as well as residualization of free 89Zr. Subsequently, the model was optimized using mouse biodistribution data, where it was observed that free 89Zr mainly residualizes in the bone and the extent of antibody distribution in certain tissues (e.g., liver and spleen) may be altered by labeling with 89Zr. The mouse PBPK model was scaled to rat, monkey, and human by simply changing the physiological parameters, and a priori simulations performed by the model were compared with the observed PK data. It was found that model predicted antibody PK in majority of the tissues in all the species superimposed over the observed data, and the model was also able to predict the PK of antibody in human tissues reasonably well. As such, the work presented here provides unprecedented evaluation of the antibody PPBK model for its ability to predict tissue PK of antibodies in the clinic. This model can be used for preclinical-to-clinical translation of antibodies and for prediction of antibody concentrations at the site-of-action in the clinic.

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.

Mechanisms of Obesity-Induced Changes in Pharmacokinetics of IgG in Rats

PURPOSE

To evaluate how obesity affects the pharmacokinetics of human IgG following subcutaneous (SC) and intravenous (IV) administration to rats and the homeostasis of endogenous rat IgG.

METHODS

Differences in body weight and size, body composition, and serum concentration of endogenous rat IgG in male Zucker obese (ZUC-FA/FA) and control (ZUC-LEAN) rats were measured from the age of 5 weeks up to 30 weeks. At the age of 23-24 weeks animals received a single IV or SC dose of human IgG (1 g/kg of total body weight), and serum pharmacokinetics was followed for 7 weeks. A mechanistic model linking obesity-related changes in pharmacokinetics with animal growth and changes in body composition was developed.

RESULTS

Significant differences were observed in both endogenous and exogenous IgG pharmacokinetics between obese and control groups. The AUC for human IgG was lower in obese groups (57.6% of control after IV and 48.1% after SC dosing), and clearance was 1.75-fold higher in obese animals. The mechanistic population model successfully captured the data and included several major components: endogenous rat IgG homeostasis with age-dependent synthesis rate; competition of human IgG and endogenous rat IgG for FcRn binding and its effect on endogenous rat IgG concentrations following injection of a high dose of human IgG; and the effect of body size and composition (changing over time and dependent on the obesity status) on pharmacokinetic parameters.

CONCLUSIONS

We identified important obesity-induced changes in the pharmacokinetics of IgG. Results can potentially facilitate optimization of the dosing of IgG-based therapeutics in the obese population.

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.

Modeling Blood–Brain Barrier Permeability to Solutes and Drugs In Vivo

Our understanding of the pharmacokinetic principles governing the uptake of endogenous substances, xenobiotics, and biologicals across the blood–brain barrier (BBB) has advanced significantly over the past few decades. There is now a spectrum of experimental techniques available in experimental animals and humans which, together with pharmacokinetic models of low to high complexity, can be applied to describe the transport processes at the BBB of low molecular weight agents and macromolecules. This review provides an overview of the models in current use, from initial rate uptake studies over compartmental models to physiologically based models and points out the advantages and shortcomings associated with the different methods. A comprehensive pharmacokinetic profile of a compound with respect to brain exposure requires the knowledge of BBB uptake clearance, intra-brain distribution, and extent of equilibration across the BBB. The application of proper pharmacokinetic analysis and suitable models is a requirement not only in the drug development process, but in all of the studies where the brain uptake of drugs or markers is used to make statements about the function or integrity of the BBB.

ADME and DMPK considerations for the discovery and development of antibody drug conjugates (ADCs)

 The therapeutic concept of antibody drug conjugates (ADCs) is to selectively target tumour cells with small molecule cytotoxic drugs to maximise cell kill benefit and minimise healthy tissue toxicity.An ADC generally consists of an antibody that targets a protein on the surface of tumour cells chemically linked to a warhead small molecule cytotoxic drug.To deliver the warhead to the tumour cell, the antibody must bind to the target protein and in general be internalised into the cell. Following internalisation, the cytotoxic agent can be released in the endosomal or lysosomal compartment (via different mechanisms). Diffusion or transport out of the endosome or lysosome allows the cytotoxic drug to express its cell-killing pharmacology. Alternatively, some ADCs (e.g. EDB-ADCs) rely on extracellular cleavage releasing membrane permeable warheads.One potentially important aspect of the ADC mechanism is the 'bystander effect' whereby the cytotoxic drug released in the targeted cell can diffuse out of that cell and into other (non-target expressing) tumour cells to exert its cytotoxic effect. This is important as solid tumours tend to be heterogeneous and not all cells in a tumour will express the targeted protein.The combination of large and small molecule aspects in an ADC poses significant challenges to the disposition scientist in describing the ADME properties of the entire molecule.This article will review the ADC landscape and the ADME properties of successful ADCs, with the aim of outlining best practice and providing a perspective of how the field can further facilitate the discovery and development of these important therapeutic modalities.

Leveraging a physiologically-based quantitative translational modeling platform for designing B cell maturation antigen-targeting bispecific T cell engagers for treatment of multiple myeloma

Bispecific T cell engagers (TCEs) are an emerging anti-cancer modality that redirects cytotoxic T cells to tumor cells expressing tumor-associated antigens (TAAs), thereby forming immune synapses to exert anti-tumor effects. Designing pharmacokinetically acceptable TCEs and optimizing their size presents a considerable protein engineering challenge, particularly given the complexity of intercellular bridging between T cells and tumor cells. Therefore, a physiologically-relevant and clinically-verified computational modeling framework is of crucial importance to understand the protein engineering trade-offs. In this study, we developed a quantitative, physiologically-based computational framework to predict immune synapse formation for a variety of molecular formats of TCEs in tumor tissues. Our model incorporates a molecular size-dependent biodistribution using the two-pore theory, extravasation of T cells and hematologic cancer cells, mechanistic bispecific intercellular binding of TCEs, and competitive inhibitory interactions by shed targets. The biodistribution of TCEs was verified by positron emission tomography imaging of [89Zr]AMG211 (a carcinoembryonic antigen-targeting TCE) in patients. Parameter sensitivity analyses indicated that immune synapse formation was highly sensitive to TAA expression, degree of target shedding, and binding selectivity to tumor cell surface TAAs over shed targets. Notably, the model suggested a “sweet spot” for TCEs’ CD3 binding affinity, which balanced the trapping of TCEs in T-cell-rich organs. The final model simulations indicated that the number of immune synapses is similar (~55/tumor cell) between two distinct clinical stage B cell maturation antigen (BCMA)-targeting TCEs, PF-06863135 in an IgG format and AMG420 in a BiTE format, at their respective efficacious doses in multiple myeloma patients. This result demonstrates the applicability of the developed computational modeling framework to molecular design optimization and clinical benchmarking for TCEs, thus suggesting that this framework can be applied to other targets to provide a quantitative means to facilitate model-informed best-in-class TCE discovery and development.

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.

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.

Pharmacokinetic-pharmacodynamic modelling of the anti-FcRn monoclonal antibody rozanolixizumab: Translation from preclinical stages to the clinic

Rozanolixizumab is a fully humanized high‐affinity anti‐human neonatal Fc receptor (FcRn) monoclonal antibody (mAb) that accelerates the removal of circulating immunoglobulin G (IgG), including pathogenic IgG autoantibodies, via the natural lysosomal degradation pathway. The aim of this study was to develop a pharmacokinetic/pharmacodynamic (PK/PD) model characterizing the effect of rozanolixizumab on IgG levels in cynomolgus monkeys, translate it into humans to support the first‐in‐human (FIH) rozanolixizumab clinical trial study design, and, ultimately, develop a PK/PD model in humans. Simulations from the preclinical model were performed to predict IgG responses in humans and select clinically relevant doses in the FIH study. Good alignment was observed between predicted and observed reductions in IgG, which increased with increasing dose in the FIH study. The model successfully described the PK of the 4 and 7 mg/kg intravenous (i.v.) dose groups, although the PKs were underpredicted for the 1 mg/kg i.v. dose group. Updating the model with subsequent human data identified parameters that deviated from preclinical assumptions. The updated PK/PD model was able to effectively characterize the PK FcRn‐IgG nonlinear system in response to rozanolixizumab in the FIH data.

Mechanistic Modeling of a Novel Oncolytic Virus, V937, to Describe Viral Kinetic and Dynamic Processes Following Intratumoral and Intravenous Administration

V937 is an investigational novel oncolytic non-genetically modified Kuykendall strain of Coxsackievirus A21 which is in clinical development for the treatment of advanced solid tumor malignancies. V937 infects and lyses tumor cells expressing the intercellular adhesion molecule I (ICAM-I) receptor. We integrated in vitro and in vivo data from six different preclinical studies to build a mechanistic model that allowed a quantitative analysis of the biological processes of V937 viral kinetics and dynamics, viral distribution to tumor, and anti-tumor response elicited by V937 in human xenograft models in immunodeficient mice following intratumoral and intravenous administration. Estimates of viral infection and replication which were calculated from in vitro experiments were successfully used to describe the tumor response in vivo under various experimental conditions. Despite the predicted high clearance rate of V937 in systemic circulation (t_1/2 = 4.3 min), high viral replication was observed in immunodeficient mice which resulted in tumor shrinkage with both intratumoral and intravenous administration. The described framework represents a step towards the quantitative characterization of viral distribution, replication, and oncolytic effect of a novel oncolytic virus following intratumoral and intravenous administrations in the absence of an immune response. This model may further be expanded to integrate the role of the immune system on viral and tumor dynamics to support the clinical development of oncolytic viruses.

Two-pore physiologically based pharmacokinetic model validation using whole-body biodistribution of trastuzumab and different-size fragments in mice

In the past, our lab proposed a two-pore PBPK model for different-size protein therapeutics using de novo derived parameters and the model was validated using plasma PK data of different-size antibody fragments digitized from the literature (Li Z, Shah DK, J Pharmacokinet Pharmacodynam 46(3):305-318, 2009). To further validate the model using tissue distribution data, whole-body biodistribution study of 6 different-size proteins in mice were conducted. Studied molecules covered a wide MW range (13-150 kDa). Plasma PK and tissue distribution profiles is 9 tissues were measured, including heart, lung, liver, spleen, kidney, skin, muscle, small intestine, large intestine. Tumor exposure of different-size proteins were also evaluated. The PBPK model was validated by comparing percentage predictive errors (%PE) between observed and model predicted results for each type of molecule in each tissue. Model validation showed that the two-pore PBPK model was able to predict plasma, tissues and tumor PK of all studied molecules relatively well. This model could serve as a platform for developing a generic PBPK model for protein therapeutics in the future.

Current and prospective computational approaches and challenges for developing COVID-19 vaccines

SARS-CoV-2, which causes COVID-19, was first identified in humans in late 2019 and is a coronavirus which is zoonotic in origin. As it spread around the world there has been an unprecedented effort in developing effective vaccines. Computational methods can be used to speed up the long and costly process of vaccine development. Antigen selection, epitope prediction, and toxicity and allergenicity prediction are areas in which computational tools have already been applied as part of reverse vaccinology for SARS-CoV-2 vaccine development. However, there is potential for computational methods to assist further. We review approaches which have been used and highlight additional bioinformatic approaches and PK modelling as in silico methods which may be useful for SARS-CoV-2 vaccine design but remain currently unexplored. As more novel viruses with pandemic potential are expected to arise in future, these techniques are not limited to application to SARS-CoV-2 but also useful to rapidly respond to novel emerging viruses.

Development of a physiologically-based pharmacokinetic model for ocular disposition of monoclonal antibodies in rabbits

Development of protein therapeutics for ocular disorders, particularly age-related macular degeneration (AMD), is a highly competitive and expanding therapeutic area. However, the application of a predictive and translatable ocular PK model to better understand ocular disposition of protein therapeutics, such as a physiologically-based pharmacokinetic (PBPK) model, is missing from the literature. Here, we present an expansion of an antibody platform PBPK model towards rabbit and incorporate a novel anatomical and physiologically relevant ocular component. Parameters describing all tissues, flows, and binding events were obtained from existing literature and fixed a priori. First, translation of the platform PBPK model to rabbit was confirmed by evaluating the model's ability to predict plasma PK of a systemically administered exogenous antibody. Then, the PBPK model with the new ocular component was validated by estimation of serum and ocular (i.e. aqueous humor, retina, and vitreous humor) PK of two intravitreally administered monoclonal antibodies. We show that the proposed PBPK model is capable of accurately (i.e. within twofold) predicting ocular exposure of antibody-based drugs. The proposed PBPK model can be used for preclinical-to-clinical translation of antibodies developed for ocular disorders, and assessment of ocular toxicity for systemically administered antibody-based therapeutics.

A Physiologically-Based Pharmacokinetic Model for the Prediction of "Half-Life Extension" and "Catch and Release" Monoclonal Antibody Pharmacokinetics

Monoclonal antibodies (mAbs) can be engineered to have “extended half‐life” and “catch and release” properties to improve target coverage. We have developed a mAb physiologically‐based pharmacokinetic model that describes intracellular trafficking, neonatal Fc receptor (FcRn) recycling, and nonspecific clearance of mAbs. We extended this model to capture target binding as a function of target affinity, expression, and turnover. For mAbs engineered to have an extended half‐life, the model was able to accurately predict the terminal half‐life (82% within 2‐fold error of the observed value) in the human FcRn transgenic (Tg32) homozygous mouse and human. The model also accurately captures the trend in pharmacokinetic and target coverage data for a set of mAbs with differing catch and release properties in the Tg32 mouse. The mechanistic nature of this model allows us to explore different engineering techniques early in drug discovery, potentially expanding the number of “druggable” targets.

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

Use of the subcutaneous (SC) route for administering monoclonal antibodies (mAbs) to treat chronic conditions has been hindered because of an incomplete understanding of fundamental mechanisms controlling mAb absorption from the SC site, and due to the limited translatability of preclinical studies. In this paper, we report on the development and evaluation of a whole-body physiologically-based model to predict mAb pharmacokinetics following SC administration. The circulatory model is based on the physiological processes governing mAb transport and includes two mAb-specific parameters representing differences in pinocytosis rate and the diffusive/convective transport rates among mAbs. At the SC administration site, two additional parameters are used to represent mAb differences in lymphatic capillary uptake and in pre-systemic clearance. Model development employed clinical intravenous (IV) plasma PK data from 20 mAbs and SC plasma PK data from 12 of these mAbs, as obtained from the literature. The resulting model reliably described both the IV and SC measured plasma concentration data. In addition, a metric based on the positive charge across the mAb's complementarity determining region vicinity was found to positively correlate with the model-based estimates of the mAb-specific parameter governing organ/tissue pinocytosis transport and with estimates of the mAb's SC lymphatic capillary clearance. These two relationships were incorporated into the model and accurately predicted the SC PK profiles of three out of four separate mAbs not included in model development. The whole-body physiologically-based model reported herein, provides a platform to characterize and predict the plasma disposition of monoclonal antibodies following SC administration in humans.

The biodistribution of therapeutic proteins: Mechanism, implications for pharmacokinetics, and methods of evaluation

Therapeutic proteins (TPs) are a diverse drug class that include monoclonal antibodies (mAbs), recombinantly expressed enzymes, hormones and growth factors, cytokines (e.g. chemokines, interleukins, interferons), as well as a wide range of engineered fusion scaffolds containing IgG1 Fc domain for half-life extension. As the pharmaceutical industry advances more potent and selective protein-based medicines through discovery and into the clinical stages of development, it has become widely appreciated that a comprehensive understanding of the mechanisms of TP biodistribution can aid this endeavor. This review aims to highlight the literature that has advanced our understanding of the determinants of TP biodistribution. A particular emphasis is placed on the multi-faceted role of the neonatal Fc receptor (FcRn) in mAb and Fc-fusion protein disposition. In addition, characterization of the TP-target interaction at the cell-level is discussed as an essential strategy to establish pharmacokinetic-pharmacodynamic (PK/PD) relationships that may lead to more informed human dose projections during clinical development. Methods for incorporation of tissue and cell-level parameters defining these characteristics into higher-order mechanistic and semi-mechanistic PK models will also be presented.

Understanding the Monoclonal Antibody Disposition after Subcutaneous Administration using a Minimal Physiologically based Pharmacokinetic Model

PURPOSE

Monoclonal antibodies (mAbs) are commonly administered by subcutaneous (SC) route. However, bioavailability is often reduced after SC administration. In addition, the sequential transfer of mAbs through the SC tissue and lymphatic system is not completely understood. Therefore, major objectives of this study were a) To understand absorption of mAbs via the lymphatic system after SC administration using physiologically based pharmacokinetic (PBPK) modeling, and b) to demonstrate application of the model for prediction of SC pharmacokinetics (PK) of mAbs.

METHODS

A minimal PBPK model was constructed using various physiological parameters related to the SC injection site and lymphatic system. The remainder of the body organs were represented using a 2-compartment model (central and peripheral compartments), with parameters derived from available intravenous (IV) PK data. The IV and SC clinical PK data of a total of 10 mAbs were obtained from literature. The SC PK data were used to estimate the lymphatic trunk-lymph node (LN) clearance.

RESULTS

The mean estimated lymphatic trunk-LN clearance obtained from 37 SC PK profiles of mAbs was 0.00213 L/h (0.001332 to 0.002928, 95% confidence intervals). The estimated lymphatic trunk-LN clearance was greater for the mAbs with higher isoelectric point (pI). In addition, the estimated clearance increased with decrease in the bioavailability.

CONCLUSION

The minimal PBPK model identified SC injection site lymph flow, afferent and efferent lymph flows, and volumes associated with the SC injection site, lymphatic capillaries and lymphatic trunk-LN as important physiological parameters governing the absorption of mAbs after SC administration. The model may be used to predict PK of mAbs using the relationship of lymphatic trunk-LN clearance and the pI. In addition, the model can be used as a bottom platform to incorporate SC and lymphatic in vitro clearance data for mAb PK prediction in the future.

A Physiologically-Based Pharmacokinetic Model for the Prediction of Monoclonal Antibody Pharmacokinetics From In Vitro Data

Monoclonal antibody (mAb) pharmacokinetics (PK) have largely been predicted via allometric scaling with little consideration for cross-species differences in neonatal Fc receptor (FcRn) affinity or clearance/distribution mechanisms. To address this, we developed a mAb physiologically-based PK model that describes the intracellular trafficking and FcRn recycling of mAbs in a human FcRn transgenic homozygous mouse and human. This model uses mAb-specific in vitro data together with species-specific FcRn tissue expression, tissue volume, and blood-flow physiology to predict mAb in vivo linear PK a priori. The model accurately predicts the terminal half-life of 90% of the mAbs investigated within a twofold error. The mechanistic nature of this model allows us to not only predict linear PK from in vitro data but also explore the PK and target binding of mAbs engineered to have pH-dependent binding to its target or FcRn and could aid in the selection of mAbs with optimal PK and pharmacodynamic properties.

A translational platform PBPK model for antibody disposition in the brain

In this manuscript, we have presented the development of a novel platform physiologically-based pharmacokinetic (PBPK) model to characterize brain disposition of mAbs in the mouse, rat, monkey and human. The model accounts for known anatomy and physiology of the brain, including the presence of distinct blood-brain barrier and blood-cerebrospinal fluid (CSF) barrier. CSF and interstitial fluid turnover, and FcRn mediated transport of mAbs are accounted for. The model was first used to characterize published and in-house pharmacokinetic (PK) data on the disposition of mAbs in rat brain, including the data on PK of mAb in different regions of brain determined using microdialysis. Majority of model parameters were fixed based on literature reported values, and only 3 parameters were estimated using rat data. The rat PBPK model was translated to mouse, monkey, and human, simply by changing the values of physiological parameters corresponding to each species. The translated PBPK models were validated by a priori predicting brain PK of mAbs in all three species, and comparing predicted exposures with observed data. The platform PBPK model was able to a priori predict all the validation PK profiles reasonably well (within threefold), without estimating any parameters. As such, the platform PBPK model presented here provides an unprecedented quantitative tool for prediction of mAb PK at the site-of-action in the brain, and preclinical-to-clinical translation of mAbs being developed against central nervous system (CNS) disorders. The proposed model can be further expanded to account for target engagement, disease pathophysiology, and novel mechanisms, to support discovery and development of novel CNS targeting mAbs.

Computer-assembled cross-species/cross-modalities two-pore physiologically based pharmacokinetic model for biologics in mice and rats

Two-pore physiologically-based pharmacokinetic (PBPK) models can be expected to describe the tissue distribution and elimination kinetics of soluble proteins, endogenous or dosed, as function of their size. In this work, we amalgamated our previous two-pore PBPK model for an inert domain antibody (dAb) in mice with the cross-species platform PBPK model for monoclonal antibodies described in literature into a unified two-pore platform that describes protein modalities of different sizes and includes neonatal Fc receptor (FcRn) mediated recycling. This unified PBPK model was parametrized for organ-specific lymph flow rates and the endosomal recycling rate constant using an extended tissue distribution time-course dataset that included an inert dAb, albumin and IgG in rats and mice. The model was evaluated by comparing the ab initio predictions for the tissue distribution and elimination properties of albumin-binding dAbs (AlbudAbs^TM) in mice and rats with the experimental observations. Due to the large number of molecular species and reactions involved in large-scale PBPK models, we have also developed and deployed a Matlab^TM script for automating the assembly of SimBiology^TM-based two-pore biologics PBPK models which drastically cuts the time and effort required for model building.

Parameter Identification for a Model of Neonatal Fc Receptor-Mediated Recycling of Endogenous Immunoglobulin G in Humans

Salvage of endogenous immunoglobulin G (IgG) by the neonatal Fc receptor (FcRn) is implicated in many clinical areas, including therapeutic monoclonal antibody kinetics, patient monitoring in IgG multiple myeloma, and antibody-mediated transplant rejection. There is a clear clinical need for a fully parameterized model of FcRn-mediated recycling of endogenous IgG to allow for predictive modeling, with the potential for optimizing therapeutic regimens for better patient outcomes. In this paper we study a mechanism-based model incorporating nonlinear FcRn-IgG binding kinetics. The aim of this study is to determine whether parameter values can be estimated using the limited in vivo human data, available in the literature, from studies of the kinetics of radiolabeled IgG in humans. We derive mathematical descriptions of the experimental observations-timecourse data and fractional catabolic rate (FCR) data-based on the underlying physiological model. Structural identifiability analyses are performed to determine which, if any, of the parameters are unique with respect to the observations. Structurally identifiable parameters are then estimated from the data. It is found that parameter values estimated from timecourse data are not robust, suggesting that the model complexity is not supported by the available data. Based upon the structural identifiability analyses, a new expression for the FCR is derived. This expression is fitted to the FCR data to estimate unknown parameter values. Using these parameter estimates, the plasma IgG response is simulated under clinical conditions. Finally a suggestion is made for a reduced-order model based upon the newly derived expression for the FCR. The reduced-order model is used to predict the plasma IgG response, which is compared with the original four-compartment model, showing good agreement. This paper shows how techniques for compartmental model analysis-structural identifiability analysis, linearization, and reparameterization-can be used to ensure robust parameter identification.

Two-pore physiologically based pharmacokinetic model with de novo derived parameters for predicting plasma PK of different size protein therapeutics

Two-pore PBPK models have been used for characterizing the PK of protein therapeutics since 1990s. However, widespread utilization of these models is hampered by the lack of a priori parameter values, which are typically estimated using the observed data. To overcome this hurdle, here we have presented the development of a two-pore PBPK model using de novo derived parameters. The PBPK model was validated using plasma PK data for different size proteins in mice. Using the "two pore theory" we were able to establish the relationship between protein size and key model parameters, such as: permeability-surface area product (PS), vascular reflection coefficient (σ), peclet number (Pe), and glomerular sieving coefficient (θ). The model accounted for size dependent changes in tissue extravasation and glomerular filtration. The model was able to a priori predict the PK of 8 different proteins: IgG (150 kDa), scFv-Fc (105 kDa), F(ab)₂ (100 kDa, minibody (80 kDa), scFv₂ (55 kDa), Fab (50 kDa), diabody (50 kDa), scFv (27 kDa), and nanobody (13 kDa). In addition, the model was able to provide unprecedented quantitative insight into the relative contribution of convective and diffusive pathway towards trans-capillary mass transportation of different size proteins. The two-pore PBPK model was also able to predict systemic clearance (CL) versus Molecular Weight relationship for different size proteins reasonably well. As such, the PBPK model proposed here represents a bottom-up systems PK model for protein therapeutics, which can serve as a generalized platform for the development of truly translational PBPK model for protein therapeutics.

Physiologically-based modeling of monoclonal antibody pharmacokinetics in drug discovery and development

Over the past few decades, monoclonal antibodies (mAbs) have become one of the most important and fastest growing classes of therapeutic molecules, with applications in a wide variety of disease areas. As such, understanding of the determinants of mAb pharmacokinetic (PK) processes (absorption, distribution, metabolism, and elimination) is crucial in developing safe and efficacious therapeutics. In the present review, we discuss the use of physiologically-based pharmacokinetic (PBPK) models as an approach to characterize the in vivo behavior of mAbs, in the context of the key PK processes that should be considered in these models. Additionally, we discuss current and potential future applications of PBPK in the drug discovery and development timeline for mAbs, spanning from identification of potential target molecules to prediction of potential drug-drug interactions. Finally, we conclude with a discussion of currently available PBPK models for mAbs that could be implemented in the drug development process.

Development and Evaluation of a Physiologically Based Pharmacokinetic Model for Predicting the Effects of Anti-FcRn Therapy on the Disposition of Endogenous IgG in Humans

This work scaled up a previously developed physiologically based pharmacokinetic model to predict the effects of anti-FcRn agents on the disposition of endogenous IgG in human subjects. Simulations were performed with the scaled model to predict the effects of single- and multiple-dose administration of anti-FcRn monoclonal antibodies (1-256 mg/kg) and high-dose intravenous immune globulin (0.4-2 g/kg). The model was evaluated for prediction accuracy through comparison to the effects of rozanolixizumab, an anti-FcRn monoclonal antibodies under current clinical evaluation, on the disposition of endogenous IgG in healthy human subjects. The model provided reasonably accurate predictions of the effects of rozanolixizumab. Prediction errors for the maximum reduction in endogenous IgG concentrations were -8.50% (90% model prediction interval: -14.0% to 1.44%), 3.33% (90% model prediction interval: -13.9% to 21.2%), and 6.85% (90% model prediction interval: -35.2% to 10.5%) for rozanolixizumab doses of 1, 4, and 7 mg/kg, respectively. Model simulations predict that anti-FcRn therapies will exhibit greater dose potency in healthy volunteers than in patients with elevated IgG production rates (e.g., as typically found in autoimmune disease). The model appears to have potential for use in assessing and predicting novel dosing strategies for anti-FcRn therapies.

Application of Physiologically Based Pharmacokinetic Modeling to Predict the Effects of FcRn Inhibitors in Mice, Rats, and Monkeys

There is a growing interest in developing inhibitors of the neonatal Fc-receptor, FcRn, for use in the treatment for humoral autoimmune conditions. We have developed a new physiologically based pharmacokinetic model that is capable of characterizing the pharmacokinetics and pharmacodynamics of anti-FcRn monoclonal antibodies (mAb) in mice, rats, and monkeys. The model includes incorporation of FcRn recycling of immune gamma globulin (IgG) in hematopoietic cells in addition to FcRn recycling of IgG in vascular endothelial cells and considers FcRn turnover and intracellular cycling. The model captured antibody disposition in wild-type and FcRn-knockout mice and rats, and also predicted the effects of intravenous immune globulin and anti-FcRn mAb on IgG disposition. Simulations predicted the change in IgG clearance in response to intravenous immune globulin with good accuracy in rats (mean prediction error of 7.15% ± 7.67%). In monkeys, prediction windows for simulated IgG concentration versus time data, as generated through Monte Carlo simulation, were able to capture the effects of anti-FcRn mAb on endogenous IgG. The model may have utility in guiding preclinical evaluations of anti-FcRn therapies in development, potentially assisting in the identification of optimal dosing strategies for this emerging class of immunosuppressive drugs.

Targeting therapeutics to endothelium: are we there yet?

Vascular endothelial cells represent an important therapeutic target in many pathologies, including inflammation, oxidative stress, and thrombosis; however, delivery of drugs to this site is often limited by the lack of specific affinity of therapeutics for these cells. Selective delivery of both small molecule drugs and therapeutic proteins to the endothelium has been achieved through the use of targeting ligands, such as monoclonal antibodies, directed against endothelial cell surface markers, particularly cell adhesion molecules (CAMs). Careful selection of target molecules and targeting agents allows for precise delivery to sites of inflammation, thereby maximizing therapeutic drug concentrations at the site of injury. A good understanding of the physiological and pathological determinants of drug and drug carrier pharmacokinetics and biodistribution may allow for a priori identification of optimal properties of drug carrier and targeting agent. Targeted delivery of therapeutics such as antioxidants and antithrombotic agents to the injured endothelium has shown efficacy in preclinical models, suggesting the potential for translation into clinical practice. As with all therapeutics, demonstration of both efficacy and safety are required for successful clinical implementation, which must be considered not only for the individual components (drug, targeting agent, etc.) but also for the sum of the parts (e.g., the drug delivery system), as unexpected toxicities may arise with complex delivery systems. While the use of endothelial targeting has not been translated into the clinic to date, the preclinical results summarized here suggest that there is hope for successful implementation of these agents in the years to come.

Physiologically Based Modeling of the Pharmacokinetics of "Catch-and-Release" Anti-Carcinoembryonic Antigen Monoclonal Antibodies in Colorectal Cancer Xenograft Mouse Models

Engineered monoclonal antibodies (mAbs) with pH-sensitive target release, or "catch-and-release" (CAR) binding, have shown promise in decreasing the extent of target-mediated mAb elimination, increasing mAb exposure, and increasing dose potency. This study developed a mechanistic physiologically based pharmacokinetic (PBPK) model to evaluate the effects of pH-sensitive CAR target binding on the disposition of anti-carcinoembryonic antigen (CEA) mAbs in mouse models of colorectal cancer. The PBPK model was qualified by comparing model-predicted plasma concentration-time data with data observed in tumor-bearing mice following the administration of T84.66, a "standard" anti-CEA mAb that demonstrates strong binding at pH 7.4 and 5.5. Further simulations evaluated the effects CAR pH-dependent binding, with decreasing CEA affinity with decreasing pH, on anti-CEA mAb plasma pharmacokinetics. Simulated data were compared with data observed for a novel CAR mAb, 10H6. The PBPK model provided precise parameter estimates, and excellent data characterization (median prediction error 18.4%) following fitting to T84.66 data. Simulations well predicted 10H6 data (median prediction error 21.4%). Sensitivity analyses demonstrated that key determinants of the disposition of CAR mAbs include the following: antigen binding affinity, the rate constant of mAb-CEA dissociation in acidified endosomes, antigen concentration, and the tumor vasculature reflection coefficient.

Tissue Physiology of Cynomolgus Monkeys: Cross-Species Comparison and Implications for Translational Pharmacology

We previously performed a comparative assessment of tissue-level vascular physiological parameters in mice and rats, two of the most commonly utilized species in translational drug development. The present work extends this effort to non-human primates by measuring tissue- and organ-level vascular volumes (V_v), interstitial volumes (V_i), and blood flow rates (Q) in cynomolgus monkeys. These measurements were accomplished by red blood cell labeling, extracellular marker infusion, and rubidium chloride bolus distribution, respectively, the same methods used in previous rodent measurements. In addition, whole-body blood volumes (BV) were determined across species. The results demonstrate that V_v, V_i, and Q, measured using our methods scale approximately by body weight across mouse, rat, and monkey in the tissues considered here, where allometric analysis allowed extrapolation to human parameters. Significant differences were observed between the values determined in this study and those reported in the literature, including V_v in muscle, brain, and skin and Q in muscle, adipose, heart, thymus, and spleen. The impact of these differences for selected tissues was evaluated via sensitivity analysis using a physiologically based pharmacokinetic model. The blood-brain barrier in monkeys was shown to be more impervious to an infused radioactive tracer, indium-111-pentetate, than in mice or rats. The body weight-normalized total BV measured in monkey agreed well with previously measured value in rats but was lower than that in mice. These findings have important implications for the common practice of scaling physiological parameters from rodents to primates in translational pharmacology.

Effects of the FcRn developmental pharmacology on the pharmacokinetics of therapeutic monoclonal IgG antibody in pediatric subjects using minimal physiologically-based pharmacokinetic modelling

The aim of this study was to investigate neonatal Fc receptor (FcRn) concentration developmental pharmacology in adult and pediatric subjects using minimal physiologically-based pharmacokinetic (mPBPK) modelling. Three types of pharmacokinetic (PK) data for three agents (endogenous/exogenous native IgG, bevacizumab and palivizumab) were used. The adult group contained six subjects with weights from 50 to 100 kg. For pediatric subjects, seven age groups were assumed, with five subjects each having the weight of 95%, 75%, 50%, 25% and 5% percentile of the population. A first evidence-based rating system to evaluate the quality of the source data used to derive pediatric-specific mPBPK model parameter was proposed. A stepwise approach was used to examine the best combination of age/weight effect on the parameters of the mPBPK model in adult and pediatric subjects. IgG synthesis rate (K_syn), extravasation rate (ER) and FcRn were fitted simultaneously to the PK of bevacizumab and native-IgG in both adult and pediatric. All fitting showed good fits based on the graphs and the coefficient of variation of the fitted parameters (< 50%). Estimated weight-normalized K_syn increased while weight-normalized FcRn and ER decreased with increasing age. The age and weight effect on FcRn were successfully estimated from the data. The final mPBPK model developed with native IgG and bevacizumab was able to predict the PK of palivizumab in pediatric subjects. Implementation of the mPBPK model enables us to analyze the relationships of age, weight, FcRn, ER and K_syn in both adult and pediatric subject. This information may benefit the understanding of complex interaction between the FcRn developmental pharmacology and PK parameters, and improve the prediction of the antibody disposition in pediatric subjects.

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.

Hematopoietic cells as site of first-pass catabolism after subcutaneous dosing and contributors to systemic clearance of a monoclonal antibody in mice

The neonatal Fc receptor (FcRn) has been demonstrated to contribute to a high bioavailability of monoclonal antibodies (mAbs). In this study, we explored the cellular sites of FcRn-mediated protection after subcutaneous (SC) and intravenous (IV) administration. SC absorption and IV disposition kinetics of a mAb were studied in hFcRn transgenic (Tg) bone marrow chimeric mice in which hFcRn was restricted to radioresistant cells or hematopoietic cells. SC bioavailabilities close to 90% were observed in hFcRn Tg mice and chimeric mice with hFcRn expression in hematopoietic cells, whereas SC bioavailabilities were markedly lower when FcRn was missing in hematopoietic cells. Our study demonstrates: 1) FcRn in radiosensitive hematopoietic cells is required for high SC bioavailability, indicating first-pass catabolism after SC administration by hematopoietic cells; 2) FcRn-mediated transcytosis or recycling by radioresistent cells is not required for high SC bioavailability; and 3) after IV administration hematopoietic and radioresistent cells contribute about equally to clearance of the mAb. A pharmacokinetic model was devised to describe a mixed elimination via radioresistent and hematopoietic cells from vascular and extravascular compartments, respectively. Overall, the study indicates a relevant role of hematopoietic cells for first-pass clearance of mAbs after SC administration and confirms their role in the overall clearance of mAbs.

Development and validation of a physiology-based model for the prediction of pharmacokinetics/toxicokinetics in rabbits

The environmental fates of pharmaceuticals and the effects of crop protection products on non-target species are subjects that are undergoing intense review. Since measuring the concentrations and effects of xenobiotics on all affected species under all conceivable scenarios is not feasible, standard laboratory animals such as rabbits are tested, and the observed adverse effects are translated to focal species for environmental risk assessments. In that respect, mathematical modelling is becoming increasingly important for evaluating the consequences of pesticides in untested scenarios. In particular, physiologically based pharmacokinetic/toxicokinetic (PBPK/TK) modelling is a well-established methodology used to predict tissue concentrations based on the absorption, distribution, metabolism and excretion of drugs and toxicants. In the present work, a rabbit PBPK/TK model is developed and evaluated with data available from the literature. The model predictions include scenarios of both intravenous (i.v.) and oral (p.o.) administration of small and large compounds. The presented rabbit PBPK/TK model predicts the pharmacokinetics (Cmax, AUC) of the tested compounds with an average 1.7-fold error. This result indicates a good predictive capacity of the model, which enables its use for risk assessment modelling and simulations.

A Minimal Physiologically Based Pharmacokinetic Model with a Nested Endosome Compartment for Novel Engineered Antibodies

We proposed here a minimal physiologically based pharmacokinetic (mPBPK) model for a group of novel engineered antibodies in mice and humans. These antibodies are designed with altered binding properties of their Fc domain with neonatal Fc receptor (FcRn) or the Fab domain with their cognate targets (recycling antibodies) in acidic endosomes. To enable simulations of such binding features in the change of antibody pharmacokinetics and its target suppression, we nested an endothelial endosome compartment in parallel with plasma compartment based on our previously established mPBPK model. The fluid-phase pinocytosis rate from plasma to endothelial endosomes was reflected by the clearance of antibodies in FcRn dysfunctional humans or FcRn-knockout mice. The endosomal recycling rate of FcRn-bound antibodies was calculated based on the reported endosomal transit time. The nonspecific catabolism in endosomes was fitted using pharmacokinetic data of a human wild-type IgG₁ adalimumab in humans and B21M in human FcRn (hFcRn) transgenic mice. The developed model adequately predicted the pharmacokinetics of infliximab, motavizumab, and an Fc variant of motavizumab in humans and the pharmacokinetics of bevacizumab, an Fc variant of bevacizumab, and a recycling antibody PH-IgG₁ and its non-pH dependent counterpart NPH-IgG₁ in hFcRn transgenic mice. Our proposed model provides a platform for evaluation of the pharmacokinetics and disposition behaviors of Fc-engineered antibodies and recycling antibodies.

Two-Pore Minimum Physiologically-based Pharmacokinetic Model to Describe the Disposition of Therapeutic Monoclonal IgG Antibody in Humans

PURPOSE

The aim of this study was to develop a two-pore minimum physiologically-based pharmacokinetic (mPBPK) model in describing the pharmacokinetic (PK) of therapeutic monoclonal antibody (TMAb) in human subjects.

METHODS

PK data used in this study were endogenous/exogenous native IgG and two TMAbs (palivizumab and Motavizumab-YTE) in normal volunteer or familial hypercatabolic hypoproteinemia (FIHH) patient. Several important components were implemented to overcome the limitations of the early mPBPK model, e.g. two-pore model to describe the transcapillary transport of IgG from vascular to interstitial space. Six mPBPK models with different osmotic reflection coefficient (OFC) of transcapillary transport, endocytosis rates (ETR) and plasma clearance for the TMAbs/IgG were tested and the best model was selected using AICc values.

RESULTS

The final model consisted of different OFC and ETR values for native IgG and TMAbs, supporting the hypothesis that the dynamics in the endosomal space had an important role in the compliant FcRn salvage mechanism to determine the clearance of TMAbs. The estimated FcRn concentration of FIHH subjects was 2.72 μmol/l. The final two-pore mPBPK model has a better performance for native IgG than previously developed mPBPK model.

CONCLUSIONS

The final two-pore mPBPK model not only overcome the limitations of the early mPBPK model but also has a better performance to describe the disposition of the IgG antibody in human subjects.

A generic whole body physiologically based pharmacokinetic model for therapeutic proteins in PK-Sim

Proteins are an increasingly important class of drugs used as therapeutic as well as diagnostic agents. A generic physiologically based pharmacokinetic (PBPK) model was developed in order to represent at whole body level the fundamental mechanisms driving the distribution and clearance of large molecules like therapeutic proteins. The model was built as an extension of the PK-Sim model for small molecules incorporating (i) the two-pore formalism for drug extravasation from blood plasma to interstitial space, (ii) lymph flow, (iii) endosomal clearance and (iv) protection from endosomal clearance by neonatal Fc receptor (FcRn) mediated recycling as especially relevant for antibodies. For model development and evaluation, PK data was used for compounds with a wide range of solute radii. The model supports the integration of knowledge gained during all development phases of therapeutic proteins, enables translation from pre-clinical species to human and allows predictions of tissue concentration profiles which are of relevance for the analysis of on-target pharmacodynamic effects as well as off-target toxicity. The current implementation of the model replaces the generic protein PBPK model available in PK-Sim since version 4.2 and becomes part of the Open Systems Pharmacology Suite.

Interstitial IgG antibody pharmacokinetics assessed by combined in vivo- and physiologically-based pharmacokinetic modelling approaches

KEY POINTS

For therapeutic antibodies, total tissue concentrations are frequently reported as a lump sum measure of the antibody in residual plasma, interstitial fluid and cells. In terms of correlating antibody exposure to a therapeutic effect, however, interstitial pharmacokinetics might be more relevant. In the present study, we collected total tissue and interstitial antibody biodistribution data in mice and assessed the composition of tissue samples aiming to correct total tissue measurements for plasma and cellular content. All data and parameters were integrated into a refined physiologically-based pharmacokinetic model for monoclonal antibodies to enable the tissue-specific description of antibody pharmacokinetics in the interstitial space. We found that antibody interstitial concentrations are highly tissue-specific and dependent on the underlying capillary structure but, in several tissues, they reach relatively high interstitial concentrations, contradicting the still-prevailing view that both the distribution to tissues and the interstitial concentrations for antibodies are generally low.

ABSTRACT

For most therapeutic antibodies, the interstitium is the target space. Although experimental methods for measuring antibody pharmacokinetics (PK) in this space are not well established, thus making quantitative assessment difficult, the interstitial antibody concentration is assumed to be low. In the present study, we combined direct quantification of antibodies in the interstitial fluid with a physiologically-based PK (PBPK) modelling approach, with the aim of better describing the PK of monoclonal antibodies in the interstitial space of different tissues. We isolated interstitial fluid by tissue centrifugation and conducted an antibody biodistribution study in mice, measuring total tissue and interstitial concentrations in selected tissues. Residual plasma, interstitial volumes and lymph flows, which are important PBPK model parameters, were assessed in vivo. We could thereby refine the PBPK modelling of monoclonal antibodies, better interpret antibody biodistribution data and more accurately predict their PK in the different tissue spaces. Our results indicate that, in tissues with discontinuous capillaries (liver and spleen), interstitial concentrations are reflected by the plasma concentration. In tissues with continuous capillaries (e.g. skin and muscle), ∼50-60% of the plasma concentration is found in the interstitial space. In the brain and kidney, on the other hand, antibodies are restricted to the vascular space. Our data may significantly impact the interpretation of biodistribution data of monoclonal antibodies and might be important when relating measured concentrations to a therapeutic effect. By contrast to the view that the antibody distribution to the interstitial space is limited, using direct measurements and model-based data interpretation, we show that high antibody interstitial concentrations are reached in most tissues.

Potential Sources of Inter-Subject Variability in Monoclonal Antibody Pharmacokinetics

Understanding inter-subject variability in drug pharmacokinetics and pharmacodynamics is important to ensure that all patients attain suitable drug exposure to achieve efficacy and avoid toxicity. Inter-subject variability in the pharmacokinetics of therapeutic monoclonal antibodies (mAbs) is generally moderate to high; however, the factors responsible for the high inter-subject variability have not been comprehensively reviewed. In this review, the extent of inter-subject variability for mAb pharmacokinetics is presented and potential factors contributing to this variability are explored and summarised. Disease status, age, sex, ethnicity, body size, genetic polymorphisms, concomitant medication, co-morbidities, immune status and multiple other patient-specific details have been considered. The inter-subject variability for mAb pharmacokinetics most likely depends on the complex interplay of multiple factors. However, studies aimed at investigating the reasons for the inter-subject variability are sparse. Population pharmacokinetic models and physiologically based pharmacokinetic models are useful tools to identify important covariates, aiding in the understanding of factors contributing to inter-subject variability. Further understanding of inter-subject variability in pharmacokinetics should aid in development of dosing regimens that are more appropriate.

Quantification of IgG monoclonal antibody clearance in tissues

Monoclonal antibodies are an important therapeutic entity, and knowledge of antibody pharmacokinetics has steadily increased over the years. Despite this effort, little is known about the extent of IgG antibody degradation in different tissues of the body. While studies have been published identifying sites of degradation with the use of residualizing and non-residualizing radiolabels, quantitative tissue clearances have not yet been derived. Here, we show that in physiologically-based pharmacokinetic (PBPK) models we can combine mouse data of Indium-111 and Iodine-125 labeled antibodies with prior physiologic knowledge to determine tissue-specific intrinsic clearances. Unspecific total tissue clearance (mL/day) in the mouse was estimated to be: liver = 4.75; brain = 0.02; gut = 0.40; heart = 0.07; kidney = 0.97; lung = 0.20; muscle = 3.02; skin = 3.89; spleen = 0.45; rest of body = 2.16. The highest catabolic activity (per g tissue) was in spleen for an FcRn wild-type antibody, but shifts to the liver for an antibody with reduced FcRn affinity. In the model developed, this shift can be explained by the liver having a greater FcRn-mediated protection capacity than the spleen. The quantification of tissue intrinsic clearances and FcRn salvage capacity increases our understanding of quantitative processes that drive the therapeutic responses of antibodies. This knowledge is critical, for instance to estimate the non-specific cellular uptake and degradation of antibodies used for targeted delivery of payloads.

Emerging Role of the Spleen in the Pharmacokinetics of Monoclonal Antibodies, Nanoparticles and Exosomes

After being absorbed, drugs distribute in the body in part to reach target tissues, in part to be disposed in tissues where they do not exert clinically-relevant effects. Therapeutically-relevant effects are usually terminated by drug metabolism and/or elimination. The role that has been traditionally ascribed to the spleen in these fundamental pharmacokinetic processes was definitely marginal. However, due to its high blood flow and to the characteristics of its microcirculation, this organ would be expected to be significantly exposed to large, new generation drugs that can hardly penetrate in other tissues with tight endothelial barriers. In the present review, we examine the involvement of the spleen in the disposition of monoclonal antibodies, nanoparticles and exosomes and the possible implications for their therapeutic efficacy and toxicity. The data that we will review lead to the conclusion that a new role is emerging for the spleen in the pharmacokinetics of new generation drugs, hence suggesting that this small, neglected organ will certainly deserve stronger attention by pharmacologists in the future.

Impact of altered endogenous IgG on unspecific mAb clearance

Immunodeficient mice are crucial models to evaluate the efficacy of monoclonal antibodies (mAbs). When studying mAb pharmacokinetics (PK), protection from elimination by binding to the neonatal Fc receptor (FcRn) is known to be a major process influencing the unspecific clearance of endogenous and therapeutic IgG. The concentration of endogenous IgG in immunodeficient mice, however is reduced, and this effect on the FcRn protection mechanism and subsequently on unspecific mAb clearance is unknown, yet of great importance for the interpretation of mAb PK data. We used a PBPK modelling approach to elucidate the influence of altered endogenous IgG concentrations on unspecific mAb clearance. To this end, we used PK data in immunodeficient mice, i.e. nude and severe combined immunodeficiency mice. To avoid impact of target-mediated clearance processes, we focussed on mAbs without affinity to a target antigen in these mice. In addition, intravenous immunoglobulin (IVIG) data of immunocompetent mice was used to study the impact of increased total IgG concentrations on unspecific therapeutic antibody clearance. The unspecific clearance is linear, whenever therapeutic IgG concentrations, i.e. mAb and IVIG concentrations are lower than FcRn; it can be non-linear if therapeutic IgG concentrations are larger than FcRn and endogenous IgG concentrations (e.g., under IVIG therapy). Unspecific mAb clearance of immunodeficient mice is effectively linear (under mAb doses as typically used in human). Studying the impact of reduced endogenous IgG concentrations on unspecific mAb clearance is of great relevance for the extrapolation to clinical species, e.g., when predicting mAb PK in immunosuppressed cancer patients.

Antigenic burden and serum IgG concentrations influence rituximab pharmacokinetics in rheumatoid arthritis patients

AIMS

Rituximab is a monoclonal antibody directed against CD20, which is approved in rheumatoid arthritis (RA). This study aimed at assessing the influence of CD19+ cell counts as target-antigen amount, and of immunoglobulin G (IgG) serum concentrations on rituximab pharmacokinetics in RA patients.

METHODS

In a cohort of 64 RA patients who had received repetitive courses of rituximab, the influence of CD19+ cell count, IgG serum concentration, body surface area, sex and disease activity score in 28 joints on rituximab pharmacokinetic parameters was assessed using a population pharmacokinetic analysis.

RESULTS

A two-compartment model, with first-order distribution and elimination best described the data. The volume of distribution of central compartment and clearance of rituximab were estimated at 4.7 l and 0.56 l day^-1 , respectively. Distribution and elimination half-lives were 0.9 days and 17.3 days, respectively. As expected, the central volume of distribution increased with body surface area (P = 0.012) and was higher in male than in female (P = 0.004). We found that the elimination rate constant (k₁₀ ) increased with CD19+ count (P = 0.00022) and IgG concentration (P = 7.4 × 10^-8 ), and that k₁₀ decreased with time (P = 0.00015), partly explained by a change in target-antigen amount.

CONCLUSIONS

The association between CD19+ count and k₁₀ may be explained by target-mediated drug disposition, while the association between IgG serum concentration and k₁₀ may be explained by a saturation of the neonatal Fc receptor at high IgG concentrations, resulting in decreased recycling of rituximab.