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

A Complete Extension of Classical Hepatic Clearance Models Using Fractional Distribution Parameter fd in Physiologically Based Pharmacokinetics

The classical organ clearance models have been proposed to relate the plasma clearance CLp to probable mechanism(s) of hepatic clearance. However, the classical models assume the intrinsic capability of drug elimination (CLu,int) that is physically segregated from the vascular blood but directly acts upon the unbound drug concentration in the blood (fubCavg), and do not handle the transit-time delay between the inlet/outlet concentrations in their closed-form clearance equations. Therefore, we propose unified model structures that can address the internal blood concentration patterns of clearance organs in a more mechanistic/physiological manner, based on the fractional distribution parameter fd operative in PBPK. The basic partial/ordinary differential equations for four classical models are revisited/modified to yield a more complete set of extended clearance models, i.e., the Rattle, Sieve, Tube, and Jar models, which are the counterparts of the dispersion, series-compartment, parallel-tube, and well-stirred models. We demonstrate the feasibility of applying the resulting extended models to isolated perfused rat liver data for 11 compounds and an example dataset for in vitro-in vivo extrapolation of the intrinsic to the systemic clearances. Based on their feasibilities to handle such real data, these models may serve as an improved basis for applying clearance models in the future.

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.

Mini-PBPK-Based Population Model and Covariate Analysis to Assess the Complex Pharmacokinetics and Pharmacodynamics of RO7449135, an Anti-KLK5/KLK7 Bispecific Antibody in Cynomolgus Monkeys

RO7449135, an anti-kallikrein (KLK)5/KLK7 bispecific antibody, is in development as a potential therapy against Netherton's syndrome (NS). In cynomolgus monkey studies, RO7449135 bound to KLK5 and KLK7, causing considerable accumulation of total KLKs, but with non-dose-proportional increase. To understand the complex PKPD, a population model with covariate analysis was developed accounting for target binding in skin and migration of bound targets from skin to blood. The covariate analysis suggested the animal batch as the categorical covariate impacting the different KLK5 synthesis rates between the repeat-dose study and single-dose study, and the dose as continuous covariate impacting the internalization rate of the binary and ternary complexes containing KLK7. To comprehend the mechanism underlying, we hypothesized that inhibition of KLK5 by RO7449135 prevented its cleavage of the pro-enzyme of KLK7 (pro-KLK7) and altered the proportion between pro-KLK7 and KLK7. Besides the pro-KLK7, RO7449135 can interact with other proteins like LEKTI through KLK7 connection in a dose-dependent manner. The different high-order complexes formed by RO7449135 interacting with pro-KLK7 or LEKTI-like proteins can be subject to faster internalization rate. Accounting for the dose and animal batch as covariates, the model-predicted free target suppression is well aligned with the visual target engagement check. The population PKPD model with covariate analysis provides the scientific input for the complex PKPD analysis, successfully predicts the target suppression in cynomolgus monkeys, and thereby can be used for the human dose projection of RO7449135.

Complex PK-PD of an engineered IL-15/IL-15Rα-Fc fusion protein in cynomolgus monkeys: QSP modeling of lymphocyte dynamics

XmAb24306 is a lymphoproliferative interleukin (IL)-15/IL-15 receptor α (IL-15Rα) Fc-fusion protein currently under clinical investigation as an immunotherapeutic agent for cancer treatment. XmAb24306 contains mutations in IL-15 that attenuate its affinity to the heterodimeric IL-15 receptor βγ (IL-15R). We observe substantially prolonged pharmacokinetics (PK) (half-life ∼ 2.5 to 4.5 days) in single- and repeat-dose cynomolgus monkey (cyno) studies compared to wild-type IL-15 (half-life ∼ 1 hour), leading to increased exposure and enhanced and durable expansion of NK cells, CD8+ T cells and CD4-CD8- (double negative [DN]) T cells. Drug clearance varied with dose level and time post-dose, and PK exposure decreased upon repeated dosing, which we attribute to increased target-mediated drug disposition (TMDD) resulting from drug-induced lymphocyte expansion (i.e., pharmacodynamic (PD)-enhanced TMDD). We developed a quantitative systems pharmacology (QSP) model to quantify the complex PKPD behaviors due to the interactions of XmAb24306 with multiple cell types (CD8+, CD4+, DN T cells, and NK cells) in the peripheral blood (PB) and lymphoid tissues. The model, which includes nonspecific drug clearance, binding to and TMDD by IL15R differentially expressed on lymphocyte subsets, and resultant lymphocyte margination/migration out of PB, expansion in lymphoid tissues, and redistribution to the blood, successfully describes the systemic PK and lymphocyte kinetics observed in the cyno studies. Results suggest that after 3 doses of every-two-week (Q2W) doses up to 70 days, the relative contributions of each elimination pathway to XmAb24306 clearance are: DN T cells > NK cells > CD8+ T cells > nonspecific clearance > CD4+ T cells. Modeling suggests that observed cellular expansion in blood results from the influx of cells expanded by the drug in lymphoid tissues. The model is used to predict lymphoid tissue expansion and to simulate PK-PD for different dose regimens. Thus, the model provides insight into the mechanisms underlying the observed PK-PD behavior of an engineered cytokine and can serve as a framework for the rapid integration and analysis of data that emerges from ongoing clinical studies in cancer patients as single-agent or given in combination.

A Machine Learning Model to Estimate Toxicokinetic Half-Lives of Per- and Polyfluoro-Alkyl Substances (PFAS) in Multiple Species

Per- and polyfluoroalkyl substances (PFAS) are a diverse group of man-made chemicals that are commonly found in body tissues. The toxicokinetics of most PFAS are currently uncharacterized, but long half-lives (t_½) have been observed in some cases. Knowledge of chemical-specific t_½ is necessary for exposure reconstruction and extrapolation from toxicological studies. We used an ensemble machine learning method, random forest, to model the existing in vivo measured t_½ across four species (human, monkey, rat, mouse) and eleven PFAS. Mechanistically motivated descriptors were examined, including two types of surrogates for renal transporters: (1) physiological descriptors, including kidney geometry, for renal transporter expression and (2) structural similarity of defluorinated PFAS to endogenous chemicals for transporter affinity. We developed a classification model for t_½ (Bin 1: <12 h; Bin 2: <1 week; Bin 3: <2 months; Bin 4: >2 months). The model had an accuracy of 86.1% in contrast to 32.2% for a y-randomized null model. A total of 3890 compounds were within domain of the model, and t_½ was predicted using the bin medians: 4.9 h, 2.2 days, 33 days, and 3.3 years. For human t_½, 56% of PFAS were classified in Bin 4, 7% were classified in Bin 3, and 37% were classified in Bin 2. This model synthesizes the limited available data to allow tentative extrapolation and prioritization.

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

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

A Non-radiometric Approach to Determine Tissue Vascular Blood Volume in Biodistribution Studies

The aim of this research was to develop a reliable non-radiometric method to measure the residual blood in tissue without the need for perfusion or radiometric measurements in biodistribution studies. It was found that the perfusion method not only was ineffective in removing blood from tissue, but also introduced additional variability in the determination of tissue drug exposure and was not reproducible across studies. In addition, the use of hemoglobin as an endogenous protein and biomarker for tissue blood content was studied and it was found that hemoglobin measurement in tissue was not a reliable and effective approach for determination of residual blood level in tissue. To evaluate an alternative method for addressing the tissue blood level in biodistribution studies, animals were dosed with a Residual Blood Determinant Reagent (RBDR) 5 min prior to tissue harvesting. The level of RBDR, an exogenous protein, was measured in whole blood homogenate and in tissue lysate. Based on the level of the RBDR, the vascular blood volume (VBV) in tissue was calculated and then the tissue exposures were corrected based on the blood volumes. The tissue VBVs measured by the RBDR method were comparable with the literature values obtained by radiometric measurements.

Investigating brain uptake of a non-targeting monoclonal antibody after intravenous and intracerebroventricular administration

Monoclonal antibodies play an important role in the treatment of various diseases. However, the development of these drugs against neurological disorders where the drug target is located in the brain is challenging and requires a good understanding of the local drug concentration in the brain. In this original research, we investigated the systemic and local pharmacokinetics in the brain of healthy rats after either intravenous (IV) or intracerebroventricular (ICV) administration of EGFRvIII-T-Cell bispecific (TCB), a bispecific monoclonal antibody. We established an experimental protocol that allows serial sampling in serum, cerebrospinal fluid (CSF) and interstitial fluid (ISF) of the prefrontal cortex in freely moving rats. For detection of drug concentration in ISF, a push-pull microdialysis technique with large pore membranes was applied. Brain uptake into CSF and ISF was characterized and quantified with a reduced brain physiologically-based pharmacokinetic model. The model allowed us to interpret the pharmacokinetic processes of brain uptake after different routes of administration. The proposed model capturing the pharmacokinetics in serum, CSF and ISF of the prefrontal cortex suggests a barrier function between the CSF and ISF that impedes free antibody transfer. This finding suggests that ICV administration may not be better suited to reach higher local drug exposure as compared to IV administration. The model enabled us to quantify the relative contribution of the blood-brain barrier (BBB) and Blood-CSF-Barrier to the uptake into the interstitial fluid of the brain. In addition, we compared the brain uptake of three monoclonal antibodies after IV dosing. In summary, the presented approach can be applied to profile compounds based on their relative uptake in the brain and provides quantitative insights into which pathways are contributing to the net exposure in the brain.

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.

Development of a Physiologically-Based Pharmacokinetic Model for Whole-Body Disposition of MMAE Containing Antibody-Drug Conjugate in Mice

PURPOSE

To quantitate and mathematically characterize the whole-body pharmacokinetics (PK) of different ADC analytes following administration of an MMAE-conjugated ADC in tumor-bearing mice.

METHODS

The PK of different ADC analytes (total antibody, total drug, unconjugated drug) was measured following administration of an MMAE-conjugated ADC in tumor-bearing mice. The PK of ADC analytes was compared with the whole-body PK of the antibody and drug obtained following administration of these molecules alone. An ADC PBPK model was developed by linking antibody PBPK model with small-molecule PBPK model, where the drug was assumed to deconjugate in DAR-dependent manner.

RESULTS

Comparison of antibody biodistribution coefficient (ABC) values for total antibody suggests that conjugation of drug did not significantly affect the PK of antibody. Comparison of tissue:plasma AUC ratio (T/P) for the conjugated drug and total antibody suggests that in certain tissues (e.g., spleen) ADC may demonstrate higher deconjugation. It was observed that the tissue distribution profile of the drug can be altered following its conjugation to antibody. For example, MMAE distribution to the liver was found to increase while its distribution to the heart was found to decrease upon conjugation to antibody. MMAE exposure in the tumor was found to increase by ~20-fold following administration as conjugate (i.e., ADC). The PBPK model was able to a priori predict the PK of all three ADC analytes in plasma, tissues, and tumor reasonably well.

CONCLUSIONS

The ADC PBPK model developed here serves as a platform for translational and clinical investigations of MMAE containing ADCs.

Effect of molecular size on interstitial pharmacokinetics and tissue catabolism of antibodies

Advances in antibody engineering have enabled the construction of novel molecular formats in diverse shapes and sizes, providing new opportunities for biologic therapies and expanding the need to understand how various structural aspects affect their distribution properties. To assess the effect of antibody size on systemic pharmacokinetics (PK) and tissue distribution with or without neonatal Fc receptor (FcRn) binding, we evaluated a series of non-mouse-binding anti-glycoprotein D monoclonal antibody formats, including IgG [~150 kDa], one-armed IgG [~100 kDa], IgG-HAHQ (attenuated FcRn binding) [~150 kDa], F(ab')₂ [~100 kDa], and F(ab) [~50 kDa]. Tissue-specific concentration-time profiles were corrected for blood content based on vascular volumes and normalized based on interstitial volumes to allow estimation of interstitial concentrations and interstitial:serum concentration ratios. Blood correction demonstrated that the contribution of circulating antibody on total uptake was greatest at early time points and for highly vascularized tissues. Tissue interstitial PK largely mirrored serum exposure profiles. Similar interstitial:serum ratios were obtained for the two FcRn-binding molecules, IgG and one-armed IgG, which reached pseudo-steady-state kinetics in most tissues. For non-FcRn-binding molecules, interstitial:serum ratios changed over time, suggesting that these molecules did not reach steady-state kinetics during the study. Furthermore, concentration-time profiles of both intact and catabolized molecule were measured by a dual tracer approach, enabling quantification of tissue catabolism and demonstrating that catabolism levels were highest for IgG-HAHQ. Overall, these data sets provide insight into factors affecting preclinical distribution and may be useful in estimating interstitial concentrations and/or catabolism in human tissues.

Kinetics of Blood-Brain Barrier Transport of Monoclonal Antibodies Targeting the Insulin Receptor and the Transferrin Receptor

Biologic drugs are large molecule pharmaceuticals that do not cross the blood-brain barrier (BBB), which is formed by the brain capillary endothelium. Biologics can be re-engineered for BBB transport as IgG fusion proteins, where the IgG domain is a monoclonal antibody (MAb) that targets an endogenous BBB transporter, such as the insulin receptor (IR) or transferrin receptor (TfR). The IR and TfR at the BBB transport the receptor-specific MAb in parallel with the transport of the endogenous ligand, insulin or transferrin. The kinetics of BBB transport of insulin or transferrin, or an IRMAb or TfRMAb, can be quantified with separate mathematical models. Mathematical models to estimate the half-time of receptor endocytosis, MAb or ligand exocytosis into brain extracellular space, or receptor recycling back to the endothelial luminal membrane were fit to the brain uptake of a TfRMAb or a IRMAb fusion protein in the Rhesus monkey. Model fits to the data also allow for estimates of the rates of association of the MAb in plasma with the IR or TfR that is embedded within the endothelial luminal membrane in vivo. The parameters generated from the model fits can be used to estimate the brain concentration profile of the MAb over time, and this brain exposure is shown to be a function of the rate of clearance of the antibody fusion protein from the plasma compartment.

Mathematical Models of Blood-Brain Barrier Transport of Monoclonal Antibodies Targeting the Transferrin Receptor and the Insulin Receptor

We develop and analyze mathematical models for receptor-mediated transcytosis of monoclonal antibodies (MAb) targeting the transferrin receptor (TfR) or the insulin receptor (IR), which are expressed at the blood-brain barrier (BBB). The mass-action kinetic model for both the TfR and IR antibodies were solved numerically to generate predictions for the concentrations of all species in all compartments considered. Using these models, we estimated the rates of MAb endocytosis into brain capillary endothelium, which forms the BBB in vivo, the rates of MAb exocytosis from the intra-endothelial compartment into brain extracellular space, and the rates of receptor recycling from the endothelial space back to the luminal endothelial plasma membrane. Our analysis highlights the optimal rates of MAb association with the targeted receptor. An important role of the endogenous ligand, transferrin (Tf) or insulin, in receptor-mediated-transport (RMT) of the associated MAb was found and was attributed to the five order magnitude difference between plasma concentrations of Tf (25,000 nM) and insulin (0.3 nM). Our modeling shows that the very high plasma concentration of Tf leads to only 5% of the endothelial TfR expressed on the luminal endothelial membrane.

Monitoring tumor cell death in murine tumor models using deuterium magnetic resonance spectroscopy and spectroscopic imaging

2H magnetic resonance spectroscopic imaging has been shown recently to be a viable technique for metabolic imaging in the clinic. We show here that 2H MR spectroscopy and spectroscopic imaging measurements of [2,3-2H2]malate production from [2,3-2H2]fumarate can be used to detect tumor cell death in vivo via the production of labeled malate. Production of [2,3-2H2]malate, following injection of [2,3-2H2]fumarate (1 g/kg) into tumor-bearing mice, was measured in a murine lymphoma (EL4) treated with etoposide, and in human breast (MDA-MB-231) and colorectal (Colo205) xenografts treated with a TRAILR2 agonist, using surface-coil localized 2H MR spectroscopy at 7 T. Malate production was also imaged in EL4 tumors using a fast 2H chemical shift imaging sequence. The malate/fumarate ratio increased from 0.016 ± 0.02 to 0.16 ± 0.14 in EL4 tumors 48 h after drug treatment (P = 0.0024, n = 3), and from 0.019 ± 0.03 to 0.25 ± 0.23 in MDA-MB-231 tumors (P = 0.0001, n = 5) and from 0.016 ± 0.04 to 0.28 ± 0.26 in Colo205 tumors (P = 0.0002, n = 5) 24 h after drug treatment. These increases were correlated with increased levels of cell death measured in excised tumor sections obtained immediately after imaging. 2H MR measurements of [2,3-2H2]malate production from [2,3-2H2]fumarate provide a potentially less expensive and more sensitive method for detecting cell death in vivo than 13C MR measurements of hyperpolarized [1,4-13C2]fumarate metabolism, which have been used previously for this purpose.

AT-527, a Double Prodrug of a Guanosine Nucleotide Analog, Is a Potent Inhibitor of SARS-CoV-2 In Vitro and a Promising Oral Antiviral for Treatment of COVID-19

The impact of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), the causative agent of COVID-19, is global and unprecedented. Although remdesivir has recently been approved by the FDA to treat SARS-CoV-2 infection, no oral antiviral is available for outpatient treatment. AT-527, an orally administered double prodrug of a guanosine nucleotide analog, was previously shown to be highly efficacious and well tolerated in hepatitis C virus (HCV)-infected subjects. Here, we report the potent in vitro activity of AT-511, the free base of AT-527, against several coronaviruses, including SARS-CoV-2. In normal human airway epithelial cells, the concentration of AT-511 required to inhibit replication of SARS-CoV-2 by 90% (EC₉₀) was 0.47 μM, very similar to its EC₉₀ against human coronavirus (HCoV)-229E, HCoV-OC43, and SARS-CoV in Huh-7 cells. Little to no cytotoxicity was observed for AT-511 at concentrations up to 100 μM. Substantial levels of the active triphosphate metabolite AT-9010 were formed in normal human bronchial and nasal epithelial cells incubated with 10 μM AT-511 (698 ± 15 and 236 ± 14 μM, respectively), with a half-life of at least 38 h. Results from steady-state pharmacokinetic and tissue distribution studies of nonhuman primates administered oral doses of AT-527, as well as pharmacokinetic data from subjects given daily oral doses of AT-527, predict that twice daily oral doses of 550 mg AT-527 will produce AT-9010 trough concentrations in human lung that exceed the EC₉₀ observed for the prodrug against SARS-CoV-2 replication. This suggests that AT-527 may be an effective treatment option for COVID-19.

Whole-Body Pharmacokinetics of Antibody in Mice Determined using Enzyme-Linked Immunosorbent Assay and Derivation of Tissue Interstitial Concentrations

Here we have reported whole-body disposition of wild-type IgG and FcRn non-binding IgG in mice, determined using Enzyme-Linked Immunosorbent Assay (ELISA). The disposition data generated using ELISA are compared with previously published biodistribution data generated using radiolabelled IgG. In addition, we introduce a novel concept of ABC_IS values, which are defined as percentage ratios of tissue interstitial and plasma AUC values. These values can help in predicting tissue interstitial concentrations of monoclonal antibodies (mAbs) based on the plasma concentrations. Tissue interstitial concentrations derived from our study are also compared with previously reported values measured using microdialysis or centrifugation method. Lastly, the new set of biodistribution data generated using ELISA are used to refine the PBPK model for mAbs.