Towards a platform PBPK model to characterize the plasma and tissue disposition of monoclonal antibodies in preclinical species and human

A quantitative systems pharmacology (QSP) platform for preclinical to clinical translation of in-vivo CRISPR-Cas therapy

Background: In-vivo CRISPR Cas genome editing is a complex therapy involving lipid nanoparticle (LNP), messenger RNA (mRNA), and single guide RNA (sgRNA). This novel modality requires prior modeling to predict dose-exposure-response relationships due to limited information on sgRNA and mRNA biodistribution. This work presents a QSP model to characterize, predict, and translate the Pharmacokinetics/Pharmacodynamics (PK/PD) of CRISPR therapies from preclinical species (mouse, non-human primate (NHP)) to humans using two case studies: transthyretin amyloidosis and LDL-cholesterol reduction. Methods: PK/PD data were sourced from literature. The QSP model incorporates mechanisms post-IV injection: 1) LNP binding to opsonins in liver vasculature; 2) Phagocytosis into the Mononuclear Phagocytotic System (MPS); 3) LNP internalization via endocytosis and LDL receptor-mediated endocytosis in the liver; 4) Cellular internalization and transgene product release; 5) mRNA and sgRNA disposition via exocytosis and clathrin-mediated endocytosis; 6) Renal elimination of LNP and sgRNA; 7) Exonuclease degradation of sgRNA and mRNA; 8) mRNA translation into Cas9 and RNP complex formation for gene editing. Monte-Carlo simulations were performed for 1000 subjects and showed a reduction in serum TTR. Results: The rate of internalization in interstitial layer was 0.039 1/h in NHP and 0.007 1/h in humans. The rate of exocytosis was 6.84 1/h in mouse, 2690 1/h in NHP, and 775 1/h in humans. Pharmacodynamics were modeled using an indirect response model, estimating first-order degradation rate (0.493 1/d) and TTR reduction parameters in NHP. Discussion: The QSP model effectively characterized biodistribution and dose-exposure relationships, aiding the development of these novel therapies. The utility of platform QSP model can be paramount in facilitating the discovery and development of these novel agents.

Exploring How Antibody Format Drives Clearance from the Brain

Monoclonal antibodies (mAbs) have high binding specificity and affinity, making them attractive for treating brain diseases. However, their effectiveness is limited by poor blood-brain barrier (BBB) penetration and rapid central nervous system (CNS) clearance. Our group identified blood-brain barrier modulator (BBBM) peptides that improved mAb penetration across the BBB into the brain. In this study, we investigated the pharmacokinetics of a mAb delivered to the brain using BBBMs after intravenous (IV) administration and explored the impact of antibody format (size, neonatal Fc receptor (FcRn) binding, hyaluronic acid binding) on brain clearance following direct injection into the central nervous system (CNS) via intracerebroventricular (ICV) injection. IRDye800CW-labeled antibodies were administered into C57BL/6 mice via ICV or IV injection, and organ concentrations were measured after various time points. When a mAb was coadministered with a BBBM peptide, the permeation of mAb across the BBB was increased compared to mAb alone at early time points; however, the mAb was cleared within 2 h from the brain. ICV experiments revealed that an antibody Fab fragment had a higher brain exposure than a mAb, and that a Fab fused to a hyaluronic acid binding domain (Fab-VG1) showed remarkable improvement in brain exposure. These findings suggest that BBBMs and antibody format optimization may be promising strategies for enhancing brain retention of therapeutic antibodies.

Viral Kinetics Model of SARS-CoV-2 Infection Informs Drug Discovery, Clinical Dose, and Regimen Selection

Quantitative systems pharmacology (QSP) has been an important tool to project safety and efficacy of novel or repurposed therapies for the SARS-CoV-2 virus. Here, we present a QSP modeling framework to predict response to antiviral therapeutics with three mechanisms of action (MoA): cell entry inhibitors, anti-replicatives, and neutralizing biologics. We parameterized three distinct model structures describing virus-host interaction by fitting to published viral kinetics data of untreated COVID-19 patients. The models were used to test theoretical behaviors and map therapeutic design criteria of the different MoAs, identifying the most rapid and robust antiviral activity from neutralizing biologic and anti-replicative MoAs. We found good agreement between model predictions and clinical viral load reduction observed with anti-replicative nirmatrelvir/ritonavir (Paxlovid®) and neutralizing biologics bamlanivimab and casirivimab/imdevimab (REGEN-COV®), building confidence in the modeling framework to inform a dose selection. Finally, the model was applied to predict antiviral response with ensovibep, a novel DARPin therapeutic designed as a neutralizing biologic. We developed a new in silico measure of antiviral activity, area under the curve (AUC) of free spike protein concentration, as a metric with larger dynamic range than viral load reduction. By benchmarking to bamlanivimab predictions, we justified dose levels of 75, 225, and 600 mg ensovibep to be administered intravenously in a Phase 2 clinical investigation. Upon trial completion, we found model predictions to be in good agreement with the observed patient data. These results demonstrate the utility of this modeling framework to guide the development of novel antiviral therapeutics.

PBPK model of pegylated liposomal doxorubicin to simultaneously predict the concentration-time profile of encapsulated and free doxorubicin in tissues

The objective of this study was to develop a physiologically based pharmacokinetic (PBPK) model to predict the concentrations of encapsulated and free doxorubicin in plasma and tissues in mice after intravenous injection of PEGylated liposomes (Doxil®). The PBPK model used in this study contains liposomes and free doxorubicin disposition components. The free doxorubicin disposition component was used to simulate the disposition of free doxorubicin produced by mononuclear phagocyte system (MPS)-degrading liposomes. The liver, spleen, kidneys, and lungs contain an additional MPS subcompartment. These compartments are interconnected through blood and lymphatic circulation. The model was validated strictly by four doses of external observed plasma and tissue concentration-time profiles. The fold error (FE) values were almost all within threefold. The sensitivity analysis revealed that the MPS-related parameters greatly influenced the model. The predicted in vivo distribution characteristics of the doxorubicin liposomes and doxorubicin solution were consistent with the observed values. The PBPK model was established based on the physiological mechanism and parameters of practical significance that can be measured in vitro. Thus, it can be used to study the pharmacokinetic properties of liposomes. This study also provides a reference for the establishment of liposome PBPK model.

Physiologically Based Pharmacokinetic Model of OATP1B Substrates with a Nonlinear Mixed Effect Approach: Estimating Empirical In Vitro-to-In Vivo Scaling Factors

BACKGROUND AND OBJECTIVE

Physiologically based pharmacokinetic (PBPK) models are valuable for translating in vitro absorption, distribution, metabolism, and excretion (ADME) data to predict clinical pharmacokinetics, and can enable discovery and early clinical stages of pharmaceutical research. However, in predicting pharmacokinetics of organic anion transporting polypeptide (OATP) 1B substrates based on in vitro transport and metabolism data, PBPK models typically require additional empirical in vitro-to-in vivo scaling factors (ESFs) in order to accurately recapitulate observed clinical profiles. As model simulation is very sensitive to ESFs, a critical evaluation of ESF estimation is prudent. Previously studies have applied classic 'two-stage' and 'naïve pooled data' approaches in identifying a set of compound independent ESFs. However, the 'two-stage' approach has the parameter identification issue in separately fitting data for individual compounds, while the 'naïve pooled data' approach ignores interstudy variability, leading to potentially biased ESF estimates.

METHODS

In this study, we have applied a nonlinear mixed-effect approach in estimating ESF of the PBPK model and incorporated additional data from 86 runs of in vitro uptake assay and 49 clinical studies of 12 training compounds in model development to further enhance the translation of in vitro data to predict the pharmacokinetics of OATP1B substrate drugs. To test predication accuracy of the model, a 'leave-one-out' analysis has been performed.

RESULTS

The established model can reasonably describe the clinical observations, with both mean values and interstudy variabilities quantified for ESF and volume of distribution parameters. The mean estimates are largely consistent with values in the previous reports. The interstudy variabilities of these parameters are estimated to be at least 50% (as coefficient of variation). Most compounds can be reasonably predicted in the 'leave-one-out' analysis.

CONCLUSION

This study improves the confidence in predicting the pharmacokinetics of OATP1B substrates in individual studies of small sample sizes, and quantifies the variability associated with the prediction.

Determining tissue distribution of the oral antileishmanial agent miltefosine: a physiologically-based pharmacokinetic modeling approach

Miltefosine (MTS) is the only approved oral drug for treating leishmaniasis caused by intracellular Leishmania parasites that localize in macrophages of the liver, spleen, skin, bone marrow, and lymph nodes. MTS is extensively distributed in tissues and has prolonged elimination half-lives due to its high plasma protein binding, slow metabolic clearance, and minimal urinary excretion. Thus, understanding and predicting the tissue distribution of MTS help assess therapeutic and toxicologic outcomes of MTS, especially in special populations, e.g., pediatrics. In this study, a whole-body physiologically-based pharmacokinetic (PBPK) model of MTS was built on mice and extrapolated to rats and humans. MTS plasma and tissue concentration data obtained by intravenous and oral administration to mice were fitted simultaneously to estimate model parameters. The resulting high tissue-to-plasma partition coefficient values corroborate extensive distribution in all major organs except the bone marrow. Sensitivity analysis suggests that plasma exposure is most susceptible to changes in fraction unbound in plasma. The murine oral-PBPK model was further validated by assessing overlay of simulations with plasma and tissue profiles obtained from an independent study. Subsequently, the murine PBPK model was extrapolated to rats and humans based on species-specific physiological and drug-related parameters, as well as allometrically scaled parameters. Fold errors for pharmacokinetic parameters were within acceptable range in both extrapolated models, except for a slight underprediction in the human plasma exposure. These animal and human PBPK models are expected to provide reliable estimates of MTS tissue distribution and assist dose regimen optimization in special populations.

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.

An asymptotic description of a basic FcRn-regulated clearance mechanism and its implications for PBPK modelling of large antibodies

A basic FcRn-regulated clearance mechanism is investigated using the method of matched asymptotic expansions. The broader aim of the work is to obtain further insight on the mechanism, thereby providing theoretical support for future pharmacologically-based pharmacokinetic modelling efforts. The corresponding governing equations are first non-dimensionalised and the order of magnitudes of the model parameters are assessed based on their values reported in the literature. Under the assumption of high FcRn-binding affinity, analytical approximations are derived that are valid over the characteristic phases of the problem. Additionally, relatively simple equations relating clearance and AUC to physiological model parameters are derived, which are valid over the longest characteristic time scale of the problem. For lower to moderate doses clearance is effectively linear, whereas for higher doses it is nonlinear. It is shown that for all doses sufficiently high the leading-order approximation for the IgG concentration in plasma, over the longest characteristic time scale, is independent of the initial dose. This is because IgG that is in 'excess' of FcRn is eliminated over a time scale much shorter than that of the terminal phase. In conclusion, analytical approximations of the basic FcRn mechanism have been derived using matched asymptotic expansions, leading to a simple equation relating clearance to FcRn binding affinity, the ratio of degradation and FcRn concentration, and the volumes of the system.

Physiologically Based Pharmacokinetic Modeling: The Reversible Metabolism and Tissue-Specific Partitioning of Methylprednisolone and Methylprednisone in Rats

The pharmacokinetics (PK) of methylprednisolone (MPL) exhibited tissue-specific saturable binding and reversible conversion with its metabolite, methylprednisone (MPN). Blood and 11 tissues were collected in male rats after intravenous (i.v.) bolus doses of 50 mg/kg MPL and 20 mg/kg MPN and upon i.v. infusion of MPL and MPN at 0.3, 3, and 10 mg/h per kg. The concentrations of MPL and MPN were simultaneously measured. A comprehensive physiologically based pharmacokinetic (PBPK) model was applied to describe the plasma and tissue profiles and estimate PK parameters of the MPL/MPN interconversion system. Both dosed and formed MPL and MPN were in rapid equilibrium or achieved steady-state rapidly in plasma and tissues. MPL tissue partitioning was nonlinear, with highest capacity in liver (322.9 ng/ml) followed by kidney, heart, intestine, skin, spleen, bone, brain, muscle, and lowest in adipose (2.74 ng/ml) and displaying high penetration in lung. The tissue partition coefficient of MPN was linear but widely variable (0.15∼5.38) across most tissues, with nonlinear binding in liver and kidney. The conversion of MPL to MPN occurred in kidney, lung, and intestine with total clearance of 429 ml/h, and the back conversion occurred in liver and kidney at 1342 ml/h. The irreversible elimination clearance of MPL was 789 ml/h from liver and that of MPN was 2758 ml/h with liver accounting for 44%, lung 35%, and kidney 21%. The reversible metabolism elevated MPL exposure in rats by 13%. This highly complex PBPK model provided unique and comprehensive insights into the disposition of a major corticosteroid. SIGNIFICANCE STATEMENT: Our dual physiologically based pharmacokinetic (PBPK) study and model of methylprednisolone/methylprednisone (MPL/MPN) with multiple complexities reasonably characterized and parameterized their disposition, and provided greater insights into the interpretation of their pharmacodynamics in rats. Drug knowledge gained in this study may be translatable to higher-order species to appreciate the clinical utility of MPL. The complex model itself is instructive for advanced PBPK analysis of drugs with reversible metabolism and/or nonlinear tissue partitioning features.

Comparison of dose selection based on target engagement versus inhibition of receptor-ligand interaction for checkpoint inhibitors

Immune checkpoint inhibitors block the interaction between a receptor on one cell and its ligand on another cell, thus preventing the transduction of an immunosuppressive signal. While inhibition of the receptor-ligand interaction is key to the pharmacological activity of these drugs, it can be technically challenging to measure these intercellular interactions directly. Instead, target engagement (or receptor occupancy) is commonly measured, but may not always be an accurate predictor of receptor-ligand inhibition, and can be misleading when used to inform clinical dose projections for this class of drugs. In this study, a mathematical model explicitly representing the intercellular receptor-ligand interaction is used to compare dose prediction based on target engagement or receptor-ligand inhibition for two checkpoint inhibitors, atezolizumab and magrolimab. For atezolizumab, there is little difference between target engagement and receptor-ligand inhibition, but for magrolimab, the model predicts that receptor-ligand inhibition is significantly less than target engagement. The key variables explaining the difference between these two drugs are the relative concentrations of the target receptors and their ligands. Drug-target affinity and receptor-ligand affinity can also have divergent effects on target engagement and inhibition. These results suggest that it is important to consider ligand-receptor inhibition in addition to target engagement and demonstrate the impact of using modeling for efficacious dose estimation.

A quantitative systems pharmacology model of plasma kallikrein-kinin system dysregulation in hereditary angioedema

Hereditary angioedema (HAE) due to C1-inhibitor deficiency is a rare, debilitating, genetic disorder characterized by recurrent, unpredictable, attacks of edema. The clinical symptoms of HAE arise from excess bradykinin generation due to dysregulation of the plasma kallikrein-kinin system (KKS). A quantitative systems pharmacology (QSP) model that mechanistically describes the KKS and its role in HAE pathophysiology was developed based on HAE attacks being triggered by autoactivation of factor XII (FXII) to activated FXII (FXIIa), resulting in kallikrein production from prekallikrein. A base pharmacodynamic model was constructed and parameterized from literature data and ex vivo assays measuring inhibition of kallikrein activity in plasma of HAE patients or healthy volunteers who received lanadelumab. HAE attacks were simulated using a virtual patient population, with attacks recorded when systemic bradykinin levels exceeded 20 pM. The model was validated by comparing the simulations to observations from lanadelumab and plasma-derived C1-inhibitor clinical trials. The model was then applied to analyze the impact of nonadherence to a daily oral preventive therapy; simulations showed a correlation between the number of missed doses per month and reduced drug effectiveness. The impact of reducing lanadelumab dosing frequency from 300 mg every 2 weeks (Q2W) to every 4 weeks (Q4W) was also examined and showed that while attack rates with Q4W dosing were substantially reduced, the extent of reduction was greater with Q2W dosing. Overall, the QSP model showed good agreement with clinical data and could be used for hypothesis testing and outcome predictions.

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.

Integrative model-based comparison of target site-specific antimicrobial effects: A case study with ceftaroline and lefamulin

OBJECTIVE

Predictions of antimicrobial effects typically rely on plasma-based pharmacokinetic-pharmacodynamic (PK-PD) targets, ignoring target-site concentrations and potential differences in tissue penetration between antibiotics. In this study, we applied PK-PD modelling to compare target site-specific effects of antibiotics by integrating clinical microdialysis data, in vitro time-kill curves, and antimicrobial susceptibility distributions. As a case study, we compared the effect of lefamulin and ceftaroline against methicillin-resistant Staphylococcus aureus (MRSA) at soft-tissue concentrations.

METHODS

A population PK model describing lefamulin concentrations in plasma, subcutaneous adipose and muscle tissue was developed. For ceftaroline, a similar previously reported PK model was adopted. In vitro time-kill experiments were performed with six MRSA isolates and a PD model was developed to describe bacterial growth and antimicrobial effects. The clinical PK and in vitro PD models were linked to compare antimicrobial effects of ceftaroline and lefamulin at the different target sites.

RESULTS

Considering minimum inhibitory concentration (MIC) distributions and standard dosages, ceftaroline showed superior anti-MRSA effects compared to lefamulin both at plasma and soft-tissue concentrations. Looking at the individual antibiotics, lefamulin effects were highest at soft-tissue concentrations, while ceftaroline effects were highest at plasma concentrations, emphasising the importance of considering target-site PK-PD in antibiotic treatment optimisation.

CONCLUSION

Given standard dosing regimens, ceftaroline appeared more effective than lefamulin against MRSA at soft-tissue concentrations. The PK-PD model-based approach applied in this study could be used to compare or explore the potential of antibiotics for specific indications or in populations with unique target-site PK.

Mathematical modeling of endogenous and exogenously administered T cell recirculation in mouse and its application to pharmacokinetic studies of cell therapies

Introduction

In vivo T cell migration has been of interest to scientists for the past 60 years. T cell kinetics are important in the understanding of the immune response to infectious agents. More recently, adoptive T cell therapies have proven to be a most promising approach to treating a wide range of diseases, including autoimmune and cancer diseases, whereby the characterization of cellular kinetics represents an important step towards the prediction of therapeutic efficacy.

Methods

Here, we developed a physiologically-based pharmacokinetic (PBPK) model that describes endogenous T cell homeostasis and the kinetics of exogenously administered T cells in mouse. Parameter calibration was performed using a nonlinear fixed-effects modeling approach based on published data on T cell kinetics and steady-state levels in different tissues of mice. The Partial Rank Correlation Coefficient (PRCC) method was used to perform a global sensitivity assessment. To estimate the impact of kinetic parameters on exogenously administered T cell dynamics, a local sensitivity analysis was conducted.

Results

We simulated the model to analyze cellular kinetics following various T cell doses and frequencies of CCR7+ T cells in the population of infused lymphocytes. The model predicted the effects of T cell numbers and of population composition of infused T cells on the resultant concentration of T cells in various organs. For example, a higher percentage of CCR7+ T cells among exogenously administered T lymphocytes led to an augmented accumulation of T cells in the spleen. The model predicted a linear dependence of T cell dynamics on the dose of adoptively transferred T cells.

Discussion

The mathematical model of T cell migration presented here can be integrated into a multi-scale model of the immune system and be used in a preclinical setting for predicting the distribution of genetically modified T lymphocytes in various organs, following adoptive T cell therapies.

A PBPK model for PRRT with [177Lu]Lu-DOTA-TATE: Comparison of model implementations in SAAM II and MATLAB/SimBiology

Physiologically based pharmacokinetic (PBPK) models offer the ability to simulate and predict the biodistribution of radiopharmaceuticals and have the potential to enable individualised treatment planning in molecular radiotherapy. The objective of this study was to develop and implement a whole-body compartmental PBPK model for peptide receptor radionuclide therapy (PRRT) with [177Lu]Lu-DOTA-TATE in SimBiology to allow for more complex analyses. The correctness of the model implementation was ensured by comparing its outputs, such as the time-integrated activity (TIA), with those of a PBPK model implemented in SAAM II software.

METHODS

A combined PBPK model for [68Ga]Ga-DOTA-TATE and [177Lu]Lu-DOTA-TATE was developed and implemented in both SAAM II and SimBiology. A retrospective analysis of 12 patients with metastatic neuroendocrine tumours (NETs) was conducted. First, time-activity curves (TACs) and TIAs from the two software were calculated and compared for identical parameter values. Second, pharmacokinetic parameters were fitted to activity concentrations, analysed and compared.

RESULTS

The PBPK model implemented in SimBiology produced TIA results comparable to those generated by the model implemented in SAAM II, with a relative deviation of less than 0.5% when using the same input parameters. The relative deviation of the fitted TIAs was less than 5% when model parameter values were fitted to the measured activity concentrations.

CONCLUSION

The proposed PBPK model implemented in SimBiology can be used for dosimetry in radioligand therapy and TIA prediction. Its outputs are similar to those generated by the PBPK model implemented in SAAM II, confirming the correctness of the model implementation in SimBiology.

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.

Advances and challenges in immunoPET methodology

Immuno-positron emission tomography (immunoPET) enables imaging of specific targets that play a role in targeted therapy and immunotherapy, such as antigens on cell membranes, targets in the disease microenvironment, or immune cells. The most common immunoPET applications use a monoclonal antibody labeled with a relatively long-lived positron emitter such as 89Zr (T 1/2 = 78.4 h), but smaller antibody-based constructs labeled with various other positron emitting radionuclides are also being investigated. This molecular imaging technique can thus guide the development of new drugs and may have a pivotal role in selecting patients for a particular therapy. In early phase immunoPET trials, multiple imaging time points are used to examine the time-dependent biodistribution and to determine the optimal imaging time point, which may be several days after tracer injection due to the slow kinetics of larger molecules. Once this has been established, usually only one static scan is performed and semi-quantitative values are reported. However, total PET uptake of a tracer is the sum of specific and nonspecific uptake. In addition, uptake may be affected by other factors such as perfusion, pre-/co-administration of the unlabeled molecule, and the treatment schedule. This article reviews imaging methodologies used in immunoPET studies and is divided into two parts. The first part summarizes the vast majority of clinical immunoPET studies applying semi-quantitative methodologies. The second part focuses on a handful of studies applying pharmacokinetic models and includes preclinical and simulation studies. Finally, the potential and challenges of immunoPET quantification methodologies are discussed within the context of the recent technological advancements provided by long axial field of view PET/CT scanners.

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.

"Digital twins elucidate critical role of Tscm in clinical persistence of TCR-engineered cell therapy"

Despite recent progress in adoptive T cell therapy for cancer, understanding and predicting the kinetics of infused T cells remains a challenge. Multiple factors can impact the distribution, expansion, and decay or persistence of infused T cells in patients. We have developed a novel quantitative systems pharmacology (QSP) model of TCR-transgenic T cell therapy in patients with solid tumors to describe the kinetics of endogenous T cells and multiple memory subsets of engineered T cells after infusion. These T cells undergo lymphodepletion, proliferation, trafficking, differentiation, and apoptosis in blood, lymph nodes, tumor site, and other peripheral tissues. Using the model, we generated patient-matched digital twins that recapitulate the circulating T cell kinetics reported from a clinical trial of TCR-engineered T cells targeting E7 in patients with metastatic HPV-associated epithelial cancers. Analyses of key parameters influencing cell kinetics and differences among digital twins identify stem cell-like memory T cells (Tscm) cells as an important determinant of both expansion and persistence and suggest that Tscm-related differences contribute significantly to the observed variability in cellular kinetics among patients. We simulated in silico clinical trials using digital twins and predict that Tscm enrichment in the infused product improves persistence of the engineered T cells and could enable administration of a lower dose. Finally, we verified the broader relevance of the QSP model, the digital twins, and findings on the importance of Tscm enrichment by predicting kinetics for two patients with pancreatic cancer treated with KRAS G12D targeting T cell therapy. This work offers insight into the key role of Tscm biology on T cell kinetics and provides a quantitative framework to evaluate cellular kinetics for future efforts in the development and clinical application of TCR-engineered T cell therapies.

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.

Integrated PK/PD Modeling Relates Smoothened Inhibitor Biomarkers to The Heterogeneous Intratumor Disposition of Cetuximab in Pancreatic Cancer Tumor Models

Therapeutic antibodies have shown little efficacy in the treatment of pancreatic ductal adenocarcinomas (PDAC). Tumor desmoplasia, hypovascularity, and poor perfusion result in insufficient tumor cell exposure, contributing to treatment failure. Smoothened inhibitors of hedgehog signaling (sHHi) increase PDAC tumor permeability, perfusion, and drug delivery, and provide a tool to develop a quantitative, mechanistic understanding as to how the temporal dynamics of tumor priming can impact intratumor distribution of monoclonal antibodies (mAb). A linked pharmacokinetic (PK)/pharmacodynamic (PD) model was developed to integrate the plasma and tumor PK of a sHHi priming agent with its effects upon downstream stromal biomarkers Gli1, hyaluronic acid, and interstitial fluid pressure in PDAC patient-derived xenograft (PDX) tumors. In parallel, in situ tumor concentrations of cetuximab (CTX: anti-epidermal growth factor receptor; EGFR) were quantified as a marker for tumor delivery of mAb or antibody-drug conjugates. A minimal, physiologically-based pharmacokinetic (mPBPK) model was constructed to link sHHi effects upon mechanistic effectors of tumor barrier compromise with the intratumor distribution of CTX, and CTX occupancy of EGFR in tumors. Integration of the mPBPK model of mAb deposition and intratumor distribution with the PK/PD model of tumor responses to priming not only identified physiological parameters that are critical for tumor antibody distribution, but also provides insight into dosing regimens that could achieve maximal tumor disposition of therapeutic antibodies under conditions of transient PDAC tumor permeability barrier compromise that mechanistically-diverse tumor priming strategies may achieve.

Whole-Body Disposition and Physiologically Based Pharmacokinetic Modeling of Adeno-Associated Viruses and the Transgene Product

To facilitate model-informed drug development (MIDD) of adeno-associated virus (AAV) therapy, here we have developed a physiologically based pharmacokinetic (PBPK) model for AAVs following preclinical investigation in mice. After 2E11 Vg/mouse dose of AAV8 and AAV9 encoding a monoclonal antibody (mAb) gene, whole-body disposition of both the vector and the transgene mAb was evaluated over 3 weeks. At steady-state, the following tissue-to-blood (T/B) concentration ratios were found for AAV8/9: ∼50 for liver; ∼10 for heart and muscle; ∼2 for brain, lung, kidney, adipose, and spleen; ≤1 for bone, skin, and pancreas. T/B values for mAb were compared with the antibody biodistribution coefficients, and five different clusters of organs were identified based on their transgene expression profile. All the biodistribution data were used to develop a novel AAV PBPK model that incorporates: (i) whole-body distribution of the vector; (ii) binding, internalization, and intracellular processing of the vector; (iii) transgene expression and secretion; and (iv) whole-body disposition of the secreted transgene product. The model was able to capture systemic and tissue PK of the vector and the transgene-produced mAb reasonably well. Pathway analysis of the PBPK model suggested that liver, muscle, and heart are the main contributors for the secreted transgene mAb. Unprecedented PK data and the novel PBPK model developed here provide the foundation for quantitative systems pharmacology (QSP) investigations of AAV-mediated gene therapies. The PBPK model can also serve as a quantitative tool for preclinical study design and preclinical-to-clinical translation of AAV-based gene therapies.

Mechanistic Pharmacokinetics and Pharmacodynamics of GalNAc-siRNA: Translational Model Involving Competitive Receptor-Mediated Disposition and RISC-Dependent Gene Silencing Applied to Givosiran

Triantennary N-acetyl-D galactosamine (GalNAc)3-conjugated small interfering RNA (siRNA) have majorly advanced the development of RNA-based therapeutics. Chemically stabilized GalNAc-siRNAs exhibit extensive albeit capacity-limited (nonlinear) distribution into hepatocytes with additional complexities in intracellular liver disposition and pharmacology. A mechanism-based pharmacokinetic-pharmacodynamic (PK-PD) model of GalNAc-siRNA was developed to i) quantitate ASGPR-mediated disposition and downstream RNA-induced silencing complex (RISC)-dependent pharmacology following intravenous (IV) and subcutaneous (SC) dosing, ii) assess the kinetics of formed active metabolite, iii) leverage, as an example, published experimental data for givosiran, and iv) demonstrate PK translation across two preclinical species (rat and monkey) with subsequent prediction of human plasma PK. The structural model is based on competition between parent and formed active metabolite for occupancy and uptake via ASGPR into hepatocytes, intracellular sequestration and degradation, and downstream engagement of RNA-induced silencing complex (RISC) governing target mRNA degradation. The model jointly and accurately captured available concentration-time profiles of givosiran and/or AS(N-1)3' givosiran in rat and/or monkey plasma, liver, and/or kidney following givosiran administered both IV and SC. RISC-dependent gene silencing of ALAS1 mRNA was well-characterized. The model estimated an in vivo affinity (KD) value of 27.7 nM for GalNAc-ASGPR and weight-based allometric exponents of -0.27 and -0.24 for SC absorption and intracellular (endolysosomal) degradation rate constants. The model well-predicted reported givosiran plasma PK profiles in humans. PK simulations revealed net-shifts in liver-to-kidney distribution ratios with increasing IV and SC dose. Importantly, decreases in the relative liver uptake efficiency were demonstrated following IV and, to a lesser extent, following SC dosing explained by differential ASGPR occupancy profiles over time.

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.

Time-Activity data fitting in molecular Radiotherapy: Methodology and pitfalls

Absorbed radiation doses are essential in assessing the effects, e.g. safety and efficacy, of radiopharmaceutical therapy (RPT). Patient-specific absorbed dose calculations in the target or the organ at risk require multiple inputs. These include the number of disintegrations in the organ, i.e. the time-integrated activities (TIAs) of the organs, as well as other parameters describing the process of radiation energy deposition in the target tissue (i.e. mean energy per disintegration, radiation dose constants, etc). TIAs are then estimated by incorporating the area under the radiopharmaceutical's time-activity curve (TAC), which can be obtained by quantitative measurements of the biokinetics in the patient (typically based on imaging data such as planar scintigraphy, SPECT/CT, PET/CT, or blood and urine samples). The process of TAC determination/calculation for RPT generally depends on the user, e.g., the chosen number and schedule of measured time points, the selection of the fit function, the error model for the data and the fit algorithm. These decisions can strongly affect the final TIA values and thus the accuracy of calculated absorbed doses. Despite the high clinical importance of the TIA values, there is currently no consensus on processing time-activity data or even a clear understanding of the influence of uncertainties and variations in personalised RPT dosimetry related to user-dependent TAC calculation. As a first step towards minimising site-dependent variability in RPT dosimetry, this work provides an overview of quality assurance and uncertainty management considerations of the TIA estimation.

Toxicokinetics of a humanized anti-cocaine monoclonal antibody in male and female rats and lack of cross-reactivity

A humanized monoclonal antibody h2E2 designed to bind cocaine with high affinity, specificity, and a long half-life (~7 d in rats) is being developed as a treatment for cocaine use disorder. We report here a pharmacokinetic (PK) study of h2E2 using male and female rats conducted under a Good Laboratory Practice (GLP) protocol over a dose range of 40 to 1200 mg/kg. The maximum concentration measured in rat plasma (Cmax) varied proportionately to the dose administered in both male and female rats. The terminal elimination half-lives (t1/2β) were not significantly different in male and female rats at all doses tested. Importantly, this study reports pharmacokinetics for a humanized monoclonal antibody at a dose never tested before. h2E2 has a high affinity for cocaine, whereas low or no affinity was demonstrated for cocaine metabolites (all except cocaethylene), endogenous monoamines, and methamphetamine. This demonstrates its specificity and a potential lack of interactions with physiological and endocrine systems. A review of the clinical signs in single-dose toxicity studies in rats revealed no effects on the central nervous, respiratory, or cardiovascular systems following single intravenous doses of 40 to 1200 mg/kg. This study predicts that this monoclonal antibody may be safe and effective in humans.

Multiphysics-Informed Pharmacokinetic Modeling of Systemic Exposure of Intramuscularly Injected LNPs

Lipid nanoparticle (LNP) constructs have been widely developed for gene therapy delivery. Understanding local absorption and presystemic clearance kinetics of LNPs, however, remains limited. This subsequently restrains the prediction and assessment of the systemic exposure of locally injected LNPs. As such, a multiscale computational approach was developed by integrating multiphysics simulation of intramuscular absorption kinetics of LNPs with whole-body pharmacokinetics modeling, bridged by a presystemic lymphatic kinetic model. The overall framework was enabled by utilizing physiological parameters obtained from the literature and drug-related parameters derived from experiments. The multiscale modeling and simulation approach predicted the systemic exposure of LNPs administered intramuscularly, with a high degree of agreement between the predicted and the experimental data. Sensitivity analyses revealed that the local absorption rate, pinocytosis presystemic clearance rate, and lymph flow rate of the presystemic lymphatic compartment had the most significant impacts on Cmax. The study yielded refreshing perspectives on estimating systemic exposures of locally injected LNPs and their safety and effectiveness.

Physiologically based pharmacokinetic (PBPK) model that describes enhanced FcRn-dependent distribution of monoclonal antibodies (mAbs) by pI-engineering in mice

We previously reported that monoclonal antibodies (mAbs) with a high isoelectric point (pI) value tended to exhibit fast plasma clearance (CL) and large steady-state volume of distribution (Vdss) in mice. However, the positive correlation between pI, CL, and Vdss cannot be described by the reported physiologically based pharmacokinetic (PBPK) models, in which FcRn-mediated transcytosis of mAbs is set to be minimal compared to convection-mediated transport. To address this issue, physiological parameters (lymph flow rate, reflection coefficient, endothelial uptake clearance, and FcRn concentration) were optimized based on the pharmacokinetic profiles of mAbs with various pI values in wild type and FcRn-deficient (beta-2-microglobulin knockout [KO]) mice. Simulations using the PBPK model developed in this study showed a positive correlation between pI, CL and Vdss observed in wild-type mice. Therefore, this model successfully characterized our hypothetical mechanism that an electrostatic positive interaction between mAbs and the endothelial membrane enhances FcRn-mediated transcytosis of mAbs, resulting in large Vdss. We sought to determine the right contribution of the two pathways of antibody distribution to the interstitial space and established a new model that could effectively capture the effect of pI on FcRn-mediated distribution of mAbs in the body.

Balancing safety and efficacy: tuning the biodistribution and pharmacokinetics of cytokine immunotherapies

Modulating the immune system shows promise in treating various conditions such as autoimmune, malignant, inflammatory, and infectious diseases. While immunotherapies can provide significant clinical benefits, they can also trigger debilitating immune-related toxicities. Achieving a balance between safety and efficacy of immunotherapy remains a significant engineering challenge. A complex immune response can be simplified into a sequence of coordinated signals with precise spatial and temporal arrangements. Mimicking or inhibiting these signals with protein immunotherapies relies on engineering them with specific biodistribution and pharmacokinetic properties. This review summarizes principles governing the movement of therapeutic proteins (i.e. biologics), focusing on cytokine immunotherapies injected intravenously or locally, and highlights approaches and considerations to balance their efficacy and safety.

Systemically Administered Anti-uPAR Antibody Plasma and Lung ELF Pharmacokinetics Characterized by Minimal Lung PBPK Model

Antibody-based therapeutics have recently gained keen attention for the treatment of pulmonary indications. However, systemically administered antibody exposure in the lungs needs to be better understood and remains a topic of interest. In this study, we evaluated the exposure of two different uPAR (urokinase-type plasminogen activator receptor) targeting full-length monoclonal IgGs in plasma and lung epithelial lining fluid (ELF) of mice after IP and IV administration. Antibody AK17 exhibited linear pharmacokinetics (PK) in plasma and ELF at 3 and 30 mg/kg single IV dose. The average plasma and ELF half-lives for AK17 and AK21 ranged between ~321-411 h and ~230-345 h, respectively, indicating sustained systemic and lung exposure of antibodies. The average ELF to the plasma concentration ratio of antibodies was ~0.01 and ~0.03 with IP and IV dosing, respectively, over 2 weeks post single dose. We simultaneously characterized plasma and ELF PK of antibody in mice by developing a minimal lung PBPK model for antibody. This model reasonably captured the plasma and ELF PK data while estimating three parameters. The model accounts for the convective transport of antibody into the tissues via blood and lymph flow. FcRn-mediated transcytosis was incorporated into the model for antibody distribution across the lung epithelial barrier. This model serves as a platform to predict the pulmonary PK of systemically administered antibodies and to support optimal dose selection for desired exposure in the lungs as the site of action.

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.

Expansion of platform physiologically-based pharmacokinetic model for monoclonal antibodies towards different preclinical species: cats, sheep, and dogs

Monoclonal antibodies (mAbs) are becoming an important therapeutic option in veterinary medicine, and understanding the pharmacokinetic (PK) of mAbs in higher-order animal species is also important for human drug development. To better understand the PK of mAbs in these animals, here we have expanded a platform physiological-based pharmacokinetic (PBPK) model to characterize the disposition of mAbs in three different preclinical species: cats, sheep, and dogs. We obtained PK data for mAbs and physiological parameters for the three different species from the literature. We were able to describe the PK of mAbs following intravenous (IV) or subcutaneous administration in cats, IV administration in sheep, and IV administration dogs reasonably well by fixing the physiological parameters and just estimating the parameters related to the binding of mAbs to the neonatal Fc receptor. The platform PBPK model presented here provides a quantitative tool to predict the plasma PK of mAbs in dogs, cats, and sheep. The model can also predict mAb PK in different tissues where the site of action might be located. As such, the mAb PBPK model presented here can facilitate the discovery, development, and preclinical-to-clinical translation of mAbs for veterinary and human medicine. The model can also be modified in the future to account for more detailed compartments for certain organs, different pathophysiology in the animals, and target-mediated drug disposition.

Mechanistically modeling peripheral cytokine dynamics following bispecific dosing in solid tumors

Bispecific antibodies exhibit proven clinical benefit, and many bispecifics are currently in clinical development for oncology. Cytokine release syndrome (CRS) is a common clinical adverse effect observed following CD3-based bispecific dosing. However, the pathophysiology of CRS is not fully understood, and no computational model mechanistically describing clinical cytokine dynamics following bispecific dosing in solid tumors exists. Here, a quantitative systems pharmacology (QSP) model describing peripheral clinical cytokine dynamics following bispecific dosing in solid tumors is presented. Using tebentafusp as a case study, a CD3-bispecific approved for uveal melanoma, the model successfully captures the dynamics of five cytokines. The QSP model was shown to predict observed phenomena, such as cytokine maximum concentration suppression using step-up dosing regimens and the importance of on-target off-tumor binding toward CRS and toxicity. Furthermore, the QSP model provides rationale for these biological phenomena based on dynamics of immune cell activation and desensitization in tumors and healthy tissues. Overall, the QSP model structure presented here serves as a basis to infer cytokine dynamics for other CD3-based bispecifics or tumor types by altering model parameters to capture the scenario of interest, supporting applications including dose selection, candidate nomination, and disease area selection.

Development and Evaluation of Amorphous Solid Dispersion of Riluzole with PBPK Model to Simulate the Pharmacokinetic Profile

In the current work, screening of polymers viz. polyacrylic acid (PAA), polyvinyl pyrrolidone vinyl acetate (PVP VA), and hydroxypropyl methyl cellulose acetate succinate (HPMC AS) based on drug-polymer interaction and wetting property was done for the production of a stable amorphous solid dispersion (ASD) of a poorly water-soluble drug Riluzole (RLZ). PAA showed maximum interaction and wetting property hence, was selected for further studies. Solid state characterization studies confirmed the formation of ASD with PAA. Saturation solubility, dissolution profile, and in vivo pharmacokinetic data of the ASD formulation were generated in rats against its marketed tablet Rilutor. The RLZ:PAA ASD showed exponential enhancement in the dissolution of RLZ. Predicted and observed pharmacokinetic data in rats showed enhanced area under curve (AUC) and Cmax in plasma and brain with respect to Rilutor. Furthermore, a physiologically based pharmacokinetic (PBPK) model of rats for Rilutor and RLZ ASD was developed and then extrapolated to humans where physiological parameters were changed along with a biochemical parameter. The partition coefficient was kept similar in both species. The model was used to predict different exposure scenarios, and the simulated data was compared with observed data points. The PBPK model simulated Cmax and AUC was within two times the experimental data for plasma and brain. The Cmax and AUC in the brain increased with ASD compared to Rilutor for humans showing its potential in improving its biopharmaceutical performance and hence enhanced therapeutic efficacy. The model can predict the RLZ concentration in multiple compartments including plasma and liver.

Towards a platform quantitative systems pharmacology (QSP) model for preclinical to clinical translation of antibody drug conjugates (ADCs)

A next generation multiscale quantitative systems pharmacology (QSP) model for antibody drug conjugates (ADCs) is presented, for preclinical to clinical translation of ADC efficacy. Two HER2 ADCs (trastuzumab-DM1 and trastuzumab-DXd) were used for model development, calibration, and validation. The model integrates drug specific experimental data including in vitro cellular disposition data, pharmacokinetic (PK) and tumor growth inhibition (TGI) data for T-DM1 and T-DXd, as well as system specific data such as properties of HER2, tumor growth rates, and volumes. The model incorporates mechanistic detail at the intracellular level, to account for different mechanisms of ADC processing and payload release. It describes the disposition of the ADC, antibody, and payload inside and outside of the tumor, including binding to off-tumor, on-target sinks. The resulting multiscale PK model predicts plasma and tumor concentrations of ADC and payload. Tumor payload concentrations predicted by the model were linked to a TGI model and used to describe responses following ADC administration to xenograft mice. The model was translated to humans and virtual clinical trial simulations were performed that successfully predicted progression free survival response for T-DM1 and T-DXd for the treatment of HER2+ metastatic breast cancer, including differential efficacy based upon HER2 expression status. In conclusion, the presented model is a step toward a platform QSP model and strategy for ADCs, integrating multiple types of data and knowledge to predict ADC efficacy. The model has potential application to facilitate ADC design, lead candidate selection, and clinical dosing schedule optimization.

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.

A physiologically based pharmacokinetic model for V937 oncolytic virus in mice

Introduction: Oncolytic viruses (OVs) represent a novel therapeutic strategy in oncology due to their capability to selectively infect and replicate in cancer cells, triggering a direct and/or immune-induced tumor lysis. However, the mechanisms governing OV pharmacokinetics are still poorly understood. This work aims to develop a physiologically based pharmacokinetic model of the novel OV, V937, in non-tumor-bearing mice to get a quantitative understanding of its elimination and tissue uptake processes. Materials and methods: Model development was performed using data obtained from 60 mice. Viral levels were quantified from eight tissues after a single intravenous V937 dose. An external dataset was used for model validation. This test set included multiple-dose experiments with different routes of administration. V937 distribution in each organ was described using a physiological structure based on mouse-specific organ blood flows and volumes. Analyses were performed using the non-linear mixed-effects approach with NONMEM 7.4. Results: Viral levels showed a drop from 108 to 105 copies/µg RNA at day 1 in blood, reflected in a high estimate of total clearance (18.2 mL/h). A well-stirred model provided an adequate description for all organs except the muscle and heart, where a saturable uptake process improved data description. The highest numbers of viral copies were observed in the brain, lymph node, kidney, liver, lung, and spleen on the first day after injection. On the other hand, the maximum amount of viral copies in the heart, muscle, and pancreas occurred 3 days after administration. Conclusion: To the best of our knowledge, this is the first physiologically based pharmacokinetic model developed to characterize OV biodistribution, representing a relevant source of quantitative knowledge regarding the in vivo behavior of OVs. This model can be further expanded by adding a tumor compartment, where OVs could replicate.

Model-Informed Approaches and Innovative Clinical Trial Design for Adeno-Associated Viral Vector-Based Gene Therapy Product Development: A White Paper

The promise of viral vector-based gene therapy (GT) as a transformative paradigm for treating severely debilitating and life-threatening diseases is slowly coming to fruition with the recent approval of several drug products. However, they have a unique mechanism of action often necessitating a tortuous clinical development plan. Expertise in such complex therapeutic modality is still fairly limited in this emerging class of adeno-associated virus (AAV) vector-based gene therapies. Because of the irreversible mode of action and incomplete understanding of genotype-phenotype relationship and disease progression in rare diseases careful considerations should be given to GT product's benefit-risk profile. In particular, special attention needs to be paid to safe dose selection, reliable dose exposure response (including clinically relevant endpoints), or creative approaches in study design targeting small patient populations during clinical development. We believe that quantitative tools encompassed within model-informed drug development (MIDD) framework fits quite well in the development of such novel therapies, as they enable us to benefit from the totality of data approach in order to support dose selection as well as optimize clinical trial designs, end point selection, and patient enrichment. In this thought leadership paper, we provide our collective experiences, identify challenges, and suggest areas of improvement in applications of modeling and innovative trial design in development of AAV-based GT products and reflect on the challenges and opportunities for incorporating MIDD tools and more in rational development of these products.

Assessment of Clinical Response to V937 Oncolytic Virus After Intravenous or Intratumoral Administration Using Physiologically-Based Modeling

Oncolytic viruses (OVs) represent a potential therapeutic strategy in cancer treatment. However, there is currently a lack of comprehensive quantitative models characterizing clinical OV kinetics and distribution to the tumor. In this work, we present a mechanistic modeling framework for V937 OV, after intratumoral (i.t.) or intravascular (i.v.) administration in patients with cancer. A minimal physiologically-based pharmacokinetic model was built to characterize biodistribution of OVs in humans. Viral dynamics was incorporated at the i.t. cellular level and linked to tumor response, enabling the characterization of a direct OV killing triggered by the death of infected tumor cells and an indirect killing induced by the immune response. The model provided an adequate description of changes in V937 mRNA levels and tumor size obtained from phase I/II clinical trials after V937 administration. The model showed prominent role of viral clearance from systemic circulation and infectivity in addition to known tumor aggressiveness on clinical response. After i.v. administration, i.t. exposure of V937 was predicted to be several orders of magnitude lower compared with i.t. administration. These differences could be overcome if there is high virus infectivity and/or replication. Unfortunately, the latter process could not be identified at the current clinical setting. This work provides insights on selecting optimal OV considering replication rate and infectivity.

Translation of Monoclonal Antibodies Pharmacokinetics from Animal to Human Using Physiologically Based Modeling in Open Systems Pharmacology (OSP) Suite: A Retrospective Analysis of Bevacizumab

Interspecies translation of monoclonal antibodies (mAbs) pharmacokinetics (PK) in presence of target-mediated drug disposition (TMDD) is particularly challenging. Incorporation of TMDD in physiologically based PK (PBPK) modeling is recent and needs to be consolidated and generalized to provide better prediction of TMDD regarding inter-species translation during preclinical and clinical development steps of mAbs. The objective of this study was to develop a generic PBPK translational approach for mAbs using the open-source software (PK-Sim® and Mobi®). The translation of bevacizumab based on data in non-human primates (NHP), healthy volunteers (HV), and cancer patients was used as a case example for model demonstration purpose. A PBPK model for bevacizumab concentration-time data was developed using data from literature and the Open Systems Pharmacology (OSP) Suite version 10. PK-sim® was used to build the linear part of bevacizumab PK (mainly FcRn-mediated), whereas MoBi® was used to develop the target-mediated part. The model was first developed for NHP and used for a priori PK prediction in HV. Then, the refined model obtained in HV was used for a priori prediction in cancer patients. A priori predictions were within 2-fold prediction error (predicted/observed) for both area under the concentration-time curve (AUC) and maximum concentration (Cmax) and all the predicted concentrations were within 2-fold average fold error (AFE) and average absolute fold error (AAFE). Sensitivity analysis showed that FcRn-mediated distribution and elimination processes must be accounted for at all mAb concentration levels, whereas the lower the mAb concentration, the more significant the target-mediated elimination. This project is the first step to generalize the full PBPK translational approach in Model-Informed Drug Development (MIDD) of mAbs using OSP Suite.

Translational physiologically-based pharmacokinetic model for ocular disposition of monoclonal antibodies

We have previously published a PBPK model comprising the ocular compartment to characterize the disposition of monoclonal antibodies (mAbs) in rabbits. While rabbits are commonly used preclinical species in ocular research, non-human primates (NHPs) have the most phylogenetic resemblance to humans including the presence of macula in the eyes as well as higher sequence homology. However, their use in ocular research is limited due to the strict ethical guidelines. Similarly, in humans the ocular samples cannot be collected except for the tapping of aqueous humor (AH). Therefore, we have translated this rabbit model to monkeys and human species using literature-reported datasets. Parameters describing the tissue volumes, physiological flows, and FcRn-binding were obtained from the literature, or estimated by fitting the model to the data. In the monkey model, the values for the rate of lysosomal degradation for antibodies (Kdeg), intraocular reflection coefficients (σaq, σret, σcho), bidirectional rate of fluid circulation between the vitreous chamber and the aqueous chamber (QVA), and permeability-surface area product of lens (PSlens) were estimated; and were found to be 31.5 h-1, 0.7629, 0.6982, 0.9999, 1.64 × 10-5 L/h, and 4.62 × 10-7 L/h, respectively. The monkey model could capture the data in plasma, aqueous humor, vitreous humor and retina reasonably well with the predictions being within twofold of the observed values. For the human model, only the value of Kdeg was estimated to fit the model to the plasma pharmacokinetics (PK) of mAbs and was found to be 24.4 h-1 (4.14%). The human model could also capture the ocular PK data reasonably well with the predictions being within two- to threefold of observed values for the plasma, aqueous and vitreous humor. Thus, the proposed framework can be used to characterize and predict the PK of mAbs in the eye of monkey and human species following systemic and intravitreal administration. The model can also facilitate the development of new antibody-based therapeutics for the treatment of ocular diseases as well as predict ocular toxicities of such molecules following systemic administration.

Mechanism-based pharmacokinetic and pharmacodynamic modeling for bispecific antibodies: challenges and opportunities

INTRODUCTION

Unlike conventional antibodies, bispecific antibodies (bsAbs) are engineered antibody- or antibody fragment-based molecules that can simultaneously recognize two different epitopes or antigens. Over the past decade, there has been an explosion of bsAbs being developed across therapeutic areas. Development of bsAbs presents unique challenges and mechanism-based pharmacokinetic/pharmacodynamic (PK/PD) modeling has served as a powerful tool to optimize their development and realize their clinical utility.

AREAS COVERED

In this review, the guiding principles and case examples of how fit-for-purpose, mechanism-based PK/PD models have been applied to answer questions commonly encountered in bsAb development are presented. Such models characterize the key pharmacological elements of bsAbs, and they can be utilized for model-informed drug development. We also include the discussion of challenges, knowledge gaps and future direction for such models.

EXPERT OPINION

Mechanistic PK/PD modeling is a powerful tool to support the development of bsAbs. These models can be extrapolated to predict treatment outcomes based on mechanisms of action (MoA) and clinical observations to form positive learn-and-confirm cycles during drug development, due to their abilities to differentiate system- and drug-specific parameters. Meanwhile, the models should keep being adapted according to novel drug design and MoA, providing continuous opportunities for model-informed drug development.

Modeling and Simulation as a Tool to Assess Voriconazole Exposure in the Central Nervous System

Voriconazole is a triazole antifungal used empirically for the treatment of complicated meningitis associated with Cryptococcus neoformans. Biopsy studies show that the drug exhibits adequate brain penetration although levels of cerebral spinal fluid (CSF) are highly variable. Considering that CSF is one of the main surrogates for CNS exposure, the present work proposed the building of a population pharmacokinetic modeling (popPK) model able to describing the exposure achieved by voriconazole in the plasma, interstitial cerebral fluid and CSF of healthy and infected rats. The developed popPK model was described by four compartments, including total plasma, free brain and total CSF concentrations. The following PK parameters were determined: K_m = 4.76 mg/L, V_max = 1.06 mg/h, Q₁ = 2.69 L, Q_in = 0.81 h^-1 and Q_out = 0.63 h^-1. Infection was a covariate in the Michaelis-Menten constant (K_m) and intercompartmental clearance from the brain tissue compartment to central compartment (Q_out). Simulations performed with the popPK model to determine the probability of reaching the therapeutic target of fAUC > MIC showed that VRC has sufficient tissue exposure in the interstitial fluid and in the CSF for the treatment of fungal infections in these tissues at prevalent minimum inhibitory concentrations.

Development of a Minimalistic Physiologically Based Pharmacokinetic (mPBPK) Model for the Preclinical Development of Spectinamide Antibiotics

Spectinamides 1599 and 1810 are lead spectinamide compounds currently under preclinical development to treat multidrug-resistant (MDR) and extensively drug-resistant (XDR) tuberculosis. These compounds have previously been tested at various combinations of dose level, dosing frequency, and route of administration in mouse models of Mycobacterium tuberculosis (Mtb) infection and in healthy animals. Physiologically based pharmacokinetic (PBPK) modeling allows the prediction of the pharmacokinetics of candidate drugs in organs/tissues of interest and extrapolation of their disposition across different species. Here, we have built, qualified, and refined a minimalistic PBPK model that can describe and predict the pharmacokinetics of spectinamides in various tissues, especially those relevant to Mtb infection. The model was expanded and qualified for multiple dose levels, dosing regimens, routes of administration, and various species. The model predictions in mice (healthy and infected) and rats were in reasonable agreement with experimental data, and all predicted AUCs in plasma and tissues met the two-fold acceptance criteria relative to observations. To further explore the distribution of spectinamide 1599 within granuloma substructures as encountered in tuberculosis, we utilized the Simcyp granuloma model combined with model predictions in our PBPK model. Simulation results suggest substantial exposure in all lesion substructures, with particularly high exposure in the rim area and macrophages. The developed model may be leveraged as an effective tool in identifying optimal dose levels and dosing regimens of spectinamides for further preclinical and clinical development.

Safety, Pharmacokinetics, Serum Neutralizing Titers, and Immunogenicity of Adintrevimab, a Monoclonal Antibody Targeting SARS-CoV-2: A Randomized, Double-Blind, Placebo-Controlled, Phase 1 Dose-escalation Study in Healthy Adults

INTRODUCTION

Adintrevimab is a fully human immunoglobulin G1 extended half-life monoclonal antibody that was developed to have broad neutralization against SARS-CoV, SARS-CoV-2, and other SARS-like CoVs with pandemic potential. Here we report the safety, pharmacokinetics (PK), serum viral neutralizing antibody (sVNA) titers, and immunogenicity results of the first three cohorts evaluated in the first-in-human study of adintrevimab in healthy adults.

METHODS

This is a phase 1, randomized, placebo-controlled, single ascending-dose study of adintrevimab administered intramuscularly (IM) or intravenously (IV) to healthy adults aged ≥ 18-55 years with no current or prior SARS-CoV-2 infection. Participants were randomized 8:2 to adintrevimab or placebo in each of three dose cohorts: adintrevimab 300 mg IM (cohort 1), 500 mg IV (cohort 2), and 600 mg IM (cohort 3). Follow-up was 12 months. Blood samples were taken predose and at multiple time points postdose up to month 12 to assess sVNA, PK, and antidrug antibodies (ADAs).

RESULTS

Thirty participants received a single dose of adintrevimab (n = 24; 8 per cohort) or placebo (n = 6). All except one adintrevimab participant in cohort 1 completed the study. No participants in any treatment arm experienced a study drug-related adverse event. Across adintrevimab-treated participants, 11 (45.8%) experienced at least one TEAE. All but one TEAE were mild in severity, and all were either viral infection or respiratory symptoms. There were no serious adverse events, discontinuations due to adverse events, or deaths. Adintrevimab exhibited a linear and dose-proportional PK profile and extended serum half-life (mean 96, 89, and 100 days in cohorts 1, 2, and 3, respectively). Participants receiving adintrevimab demonstrated dose-dependent increased sVNA titers and breadth across multiple variants.

CONCLUSION

Adintrevimab at doses of 300 mg IM, 500 mg IV, and 600 mg IM was well tolerated in healthy adults. Adintrevimab demonstrated dose-proportional exposure, rapid development of neutralizing antibody titers, and an extended half-life.

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.

Assessment of Epithelial Lining Fluid Partitioning of Systemically Administered Monoclonal Antibodies in Rats

For systemically administered monoclonal antibodies (mAbs) with pharmacological targets in the epithelial lining fluid (ELF), information on the partitioning of mAb between plasma and ELF is instrumental for dose predictions. Bronchoalveolar lavage (BAL) combined with measurements of urea as indicator of sample dilution is often used to estimate ELF concentrations of a drug. However, unbalanced extraction of mAb and urea could potentially lead to a systematic bias in the back-calculated ELF concentration. In the present study 0.5, 1, or 4 mL phosphate-buffered saline was instilled to lungs of rats to obtain lavage samples after systemic dosing of mAb and tool small molecule (n≥4/group). Furthermore, extraction of urea, mAb and the small molecule was assessed by repeatedly lavaging the lung (n = 4). There was no statistically significant difference in the calculated partitioning into ELF between the evaluated instillation volumes. Repeated BAL demonstrated that urea and the small molecule were extracted from other sources than the ELF. In contrast, there was limited to none in-flow of mAb into the lavage fluid. The unbalanced extraction of urea and mAb could theoretically result in underestimated ELF concentrations and the calculated partitioning of 0.17±0.062 might therefore constitute a lower boundary for the true partitioning.

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.

Mechanistic incorporation of FcRn binding in plasma and endosomes in a whole body PBPK model for large molecules

Monoclonal antibodies, endogenous IgG, and serum albumin bind to FcRn in the endosome for salvaging and recycling after pinocytotic uptake, which prolongs their half-life. This mechanism has been broadly recognized and is incorporated in currently available PBPK models. Newer types of large molecules have been designed and developed, which also bind to FcRn in the plasma space for various mechanistic reasons. To incorporate FcRn binding affinity in PBPK models, binding in the plasma space and subsequent internalisation into the endosome needs to be explicitly represented. This study investigates the large molecules model in PK-Sim^® and its applicability to molecules with FcRn binding affinity in plasma. With this purpose, simulations of biologicals with and without plasma binding to FcRn were performed with the large molecule model in PK-Sim^®. Subsequently, this model was extended to ensure a more mechanistic description of the internalisation of FcRn and the FcRn-drug complexes. Finally, the newly developed model was used in simulations to explore the sensitivity for FcRn binding in the plasma space, and it was fitted to an in vivo dataset of wild-type IgG and FcRn inhibitor plasma concentrations in Tg32 mice. The extended model demonstrated a strongly increased sensitivity of the terminal half-life towards the plasma FcRn binding affinity and could successfully fit the in vivo dataset in Tg32 mice with meaningful parameter estimates.

Nanoparticle biodistribution coefficients: A quantitative approach for understanding the tissue distribution of nanoparticles

The objective of this manuscript is to provide quantitative insights into the tissue distribution of nanoparticles. Published pharmacokinetics of nanoparticles in plasma, tumor and 13 different tissues of mice were collected from literature. A total of 2018 datasets were analyzed and biodistribution of graphene oxide, lipid, polymeric, silica, iron oxide and gold nanoparticles in different tissues was quantitatively characterized using Nanoparticle Biodistribution Coefficients (NBC). It was observed that typically after intravenous administration most of the nanoparticles are accumulated in the liver (NBC = 17.56 %ID/g) and spleen (NBC = 12.1 %ID/g), while other tissues received less than 5 %ID/g. NBC values for kidney, lungs, heart, bones, brain, stomach, intestine, pancreas, skin, muscle and tumor were found to be 3.1 %ID/g, 2.8 %ID/g, 1.8 %ID/g, 0.9 %ID/g, 0.3 %ID/g, 1.2 %ID/g, 1.8 %ID/g, 1.2 %ID/g, 1.0 %ID/g, 0.6 %ID/g and 3.4 %ID/g, respectively. Significant variability in nanoparticle distribution was observed in certain organs such as liver, spleen and lungs. A large fraction of this variability could be explained by accounting for the differences in nanoparticle physicochemical properties such as size and material. A critical overview of published nanoparticle physiologically-based pharmacokinetic (PBPK) models is provided, and limitations in our current knowledge about in vitro and in vivo pharmacokinetics of nanoparticles that restrict the development of robust PBPK models is also discussed. It is hypothesized that robust quantitative assessment of whole-body pharmacokinetics of nanoparticles and development of mathematical models that can predict their disposition can improve the probability of successful clinical translation of these modalities.