Whole-body Disposition Research Articles

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 (150kDa), F(ab)2 (100kDa), minibody (80kDa), Fc-containing proteins (205, 200, 110, 105, 92, 84, 81, 65, or 60kDa), albumin conjugate (85.7kDa), albumin (67kDa), Fab (50kDa), diabody (50kDa), scFv (27kDa), dAb2 (23.5kDa), proteins with an albumin-binding domain (26, 23.5, 22, 16, 14, or 13kDa), nanobody (13kDa), and other proteins (110, 65, or 60kDa). 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.

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.

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.

Evaluating Islatravir Administered Via Microneedle Array Patch for Long-Acting HIV Pre-exposure Prophylaxis Using Physiologically Based Pharmacokinetic Modelling.

Technologies for long-acting administration of antiretrovirals (ARVs) for the prevention and treatment of HIV are at the forefront of research initiatives aiming to tackle issues surrounding drug adherence with the current standard of once-daily oral administration. Islatravir (ISL) is an emerging ARV that shows promising characteristics for long-acting prevention and treatment both orally as well as through alternative routes of administration. Microneedle array patches (MAPs) are a pain-free and discreet transdermal delivery technology that offer extended-release administration of nanoparticulate drugs. This study aimed to utilise physiologically based pharmacokinetic (PBPK) modelling to predict the pharmacokinetics resulting from ISL administered via MAP and to identify key MAP characteristics required to sustain effective concentrations over extended dosing intervals. A PBPK model describing the conversion of ISL to ISL-triphosphate (ISL-TP) and its whole-body disposition was developed and verified against observed clinical data for orally administered ISL in healthy adults. An intradermal PBPK model was integrated with the ISL PBPK model to predict the dose and nanoparticle release rate required for MAP administration strategies capable of achieving a minimum ISL-TP target concentration of 0.05 pmol/106 PBMCs over extended dosing intervals. MAP design was limited to a maximum therapeutic area of 20 cm2 with a dose loading of 4.09 mg/cm2 and a minimum duration of 3 months. Due to the lack of available clinical data, a range of nanoparticle release rates and MAP bioavailability scenarios were simulated to provide an overview of potential clinical outcomes. The ISL PBPK model was successfully verified, with predicted vs observed ratios falling within 0.5-2-fold. ISL MAP doses ranging from 15 to 80 mg were predicted to sustain ISL-TP concentrations above the minimum target concentration at 3, 6 and 12 months after administration. Nanoparticle release rate and MAP bioavailability were found to have a major impact on whether dosing strategies achieved the criteria. Minimum doses of 15 mg and 60 mg with a nanoparticle release rate of 0.0005 h-1 and bioavailability ranging from 25 to 100% were predicted to achieve effective ISL-TP concentrations up to 3 and 6 months, respectively. Doses of 15 mg and 30 mg with a nanoparticle release rate of 0.0005 h-1 were also able to attain the target concentration up to 6 months after MAP administration, albeit with a minimum bioavailability of 75% and 50%, respectively. Furthermore, when simulating a bioavailability of 100%, an 80 mg ISL MAP was predicted to sustain ISL-TP concentrations above the minimum target concentration up to 12 months after administration. The ISL PBPK model successfully predicted ISL and ISL-TP pharmacokinetics across a range of orally administered regimens. The integrated intradermal PBPK model outlined optimal MAP dose and nanoparticle release rates for effective ISL-TP concentrations up to 12 months, providing justification for further investigation of ISL as a candidate for MAP administration.

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

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. 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. 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. The ADC PBPK model developed here serves as a platform for translational and clinical investigations of MMAE containing ADCs.

Whole-Body Pharmacokinetics and Physiologically Based Pharmacokinetic Model for Monomethyl Auristatin E (MMAE).

Monomethyl auristatin E (MMAE) is one of the most commonly used payloads for developing antibody–drug conjugates (ADC). However, limited studies have comprehensively evaluated the whole-body disposition of MMAE. Consequently, here, we have investigated the whole-body pharmacokinetics (PK) of MMAE in tumor-bearing mice. We show that while MMAE is rapidly eliminated from the plasma, it shows prolonged and extensive distribution in tissues, blood cells, and tumor. Highly perfused tissues (e.g., lung, kidney, heart, liver, and spleen) demonstrated tissue-to-plasma area under the concentration curve (AUC) ratios > 20, and poorly perfused tissues (e.g., fat, pancreas, skin, bone, and muscle) had ratios from 1.3 to 2.4. MMAE distribution was limited in the brain, and tumor had 8-fold higher exposure than plasma. A physiological-based pharmacokinetic (PBPK) model was developed to characterize the whole-body PK of MMAE, which accounted for perfusion/permeability-limited transfer of drug in the tissue, blood cell distribution of the drug, tissue/tumor retention of the drug, and plasma protein binding. The model was able to characterize the PK of MMAE in plasma, tissues, and tumor simultaneously, and model parameters were estimated with good precision. The MMAE PBPK model presented here can facilitate the development of a platform PBPK model for MMAE containing ADCs and help with their preclinical-to-clinical translation and clinical dose optimization.

Current progress and limitations of AAV mediated delivery of protein therapeutic genes and the importance of developing quantitative pharmacokinetic/pharmacodynamic (PK/PD) models.

While protein therapeutics are one of the most successful class of drug molecules, they are expensive and not suited for treating chronic disorders that require long-term dosing. Adeno-associated virus (AAV) mediated in vivo gene therapy represents a viable alternative, which can deliver the genes of protein therapeutics to produce long-term expression of proteins in target tissues. Ongoing clinical trials and recent regulatory approvals demonstrate great interest in these therapeutics, however, there is a lack of understanding regarding their cellular disposition, whole-body disposition, dose-exposure relationship, exposure-response relationship, and how product quality and immunogenicity affects these important properties. In addition, there is a lack of quantitative studies to support the development of pharmacokinetic-pharmacodynamic models, which can support the discovery, development, and clinical translation of this delivery system. In this review, we have provided a state-of-the-art overview of current progress and limitations related to AAV mediated delivery of protein therapeutic genes, along with our perspective on the steps that need to be taken to improve clinical translation of this therapeutic modality.

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 ABCIS 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.

Image-guided mathematical modeling for pharmacological evaluation of nanomaterials and monoclonal antibodies.

While plasma concentration kinetics has traditionally been the predictor of drug pharmacological effects, it can occasionally fail to represent kinetics at the site of action, particularly for solid tumors. This is especially true in the case of delivery of therapeutic macromolecules (drug‐loaded nanomaterials or monoclonal antibodies), which can experience challenges to effective delivery due to particle size‐dependent diffusion barriers at the target site. As a result, disparity between therapeutic plasma kinetics and kinetics at the site of action may exist, highlighting the importance of target site concentration kinetics in determining the pharmacodynamic effects of macromolecular therapeutic agents. Assessment of concentration kinetics at the target site has been facilitated by non‐invasive in vivo imaging modalities. This allows for visualization and quantification of the whole‐body disposition behavior of therapeutics that is essential for a comprehensive understanding of their pharmacokinetics and pharmacodynamics. Quantitative non‐invasive imaging can also help guide the development and parameterization of mathematical models for descriptive and predictive purposes. Here, we present a review of the application of state‐of‐the‐art imaging modalities for quantitative pharmacological evaluation of therapeutic nanoparticles and monoclonal antibodies, with a focus on their integration with mathematical models, and identify challenges and opportunities.This article is categorized under:Therapeutic Approaches and Drug Discovery > Nanomedicine for Oncologic DiseaseDiagnostic Tools > in vivo Nanodiagnostics and ImagingNanotechnology Approaches to Biology > Nanoscale Systems in Biology

Interspecies prediction of pharmacokinetics and tissue distribution of doxorubicin by physiologically-based pharmacokinetic modeling.

The aim of the study was to develop a physiologically-based pharmacokinetic (PBPK) model to describe and predict whole-body disposition of doxorubicin following intravenous administration. The PBPK model was established using previously published data in mice and included 10 tissue compartments: lungs, heart, brain, muscle, kidneys, pancreas, intestine, liver, spleen, adipose tissue, and plasma. Individual tissues were described by either perfusion-limited or permeability-limited models. All parameters were simultaneously estimated and the final model was able to describe murine data with good precision. The model was used for predicting doxorubicin disposition in rats, rabbits, dogs, and humans using interspecies scaling approaches and was qualified using plasma and tissue observed data. Reasonable prediction of the plasma pharmacokinetics and tissue distribution was achieved across all species. In conclusion, the PBPK model developed based on a rich dataset obtained from mice, was able to reasonably predict the disposition of doxorubicin in other preclinical species and humans. Applicability of the model for special populations, such as patients with hepatic impairment, was also demonstrated. The proposed model will be a valuable tool for optimization of exposure profiles of doxorubicin in human patients.

Physiologically-based modeling and interspecies prediction of paclitaxel pharmacokinetics.

The objective was to develop a physiologically-based pharmacokinetic (PBPK) model to characterize the whole-body disposition of paclitaxel (formulated in Cremophor EL and ethanol-Taxol®) in mice and to evaluate the utility of this model for predicting pharmacokinetics in other species. Published studies that reported paclitaxel plasma and tissue concentration-time data following single intravenous bolus administration of Taxol® to mice were used; and the PBPK model included plasma, liver, lungs, kidneys, spleen, heart, gastrointestinal tract, and remainder compartments. The final model resulted in a good description of the experimental plasma and tissues data in mice, where all tissues were represented by a single compartment, except the remainder that included two sub-compartments. The predictive performance of the PBPK model was assessed by evaluating its utility in predicting pharmacokinetics of paclitaxel in rats and humans. The relationship between species body weights (mice, rats, rabbits, and humans) and plasma clearance was determined by power-based regression, and resulting allometric exponent was 0.86. The model demonstrated reasonable predictions of plasma and tissue paclitaxel concentration-time profiles in rats and plasma profiles in humans. The proposed PBPK model represents an important basis that can be further utilized for characterization of novel formulations of paclitaxel.

Prediction of Drug Disposition on the Basis of its Chemical Structure

The chemical structure of any drug determines its pharmacokinetics and pharmacodynamics. Detailed understanding of relationships between the drug chemical structure and individual disposition pathways (i.e., distribution and elimination) is required for efficient use of existing drugs and effective development of new drugs. Different approaches have been developed for this purpose, ranging from statistics-based quantitative structure-property (or structure-pharmacokinetic) relationships (QSPR) analysis to physiologically based pharmacokinetic (PBPK) models. This review critically analyzes currently available approaches for analysis and prediction of drug disposition on the basis of chemical structure. Models that can be used to predict different aspects of disposition are presented, including: (a) value of the individual pharmacokinetic parameter (e.g., clearance or volume of distribution), (b) efficiency of the specific disposition pathway (e.g., biliary drug excretion or cytochrome P450 3A4 metabolism), (c) accumulation in a specific organ or tissue (e.g., permeability of the placenta or accumulation in the brain), and (d) the whole-body disposition in the individual patients. Examples of presented pharmacological agents include "classical" low-molecular-weight compounds, biopharmaceuticals, and drugs encapsulated in specialized drug-delivery systems. The clinical efficiency of agents from all these groups can be suboptimal, because of inefficient permeability of the drug to the site of action and/or excessive accumulation in other organs and tissues. Therefore, robust and reliable approaches for chemical structure-based prediction of drug disposition are required to overcome these limitations. PBPK models are increasingly being used for prediction of drug disposition. These models can reflect the complex interplay of factors that determine drug disposition in a mechanistically correct fashion and can be combined with other approaches, for example QSPR-based prediction of drug permeability and metabolism, pharmacogenomic data and tools, pharmacokinetic-pharmacodynamic modeling approaches, etc. Moreover, the PBPK models enable detailed analysis of clinically relevant scenarios, for example the effect of the specific conditions on the time course of the analyzed drug in the individual organs and tissues, including the site of action. It is expected that further development of such combined approaches will increase their precision, enhance the effectiveness of drugs, and lead to individualized drug therapy for different patient populations (geriatric, pediatric, specific diseases, etc.).

Simulation of the pharmacokinetics of bisoprolol in healthy adults and patients with impaired renal function using whole-body physiologically based pharmacokinetic modeling

To develop and evaluate a whole-body physiologically based pharmacokinetic (WB-PBPK) model of bisoprolol and to simulate its exposure and disposition in healthy adults and patients with renal function impairment. Bisoprolol dispositions in 14 tissue compartments were described by perfusion-limited compartments. Based the tissue composition equations and drug-specific properties such as log P, permeability, and plasma protein binding published in literatures, the absorption and whole-body distribution of bisoprolol was predicted using the 'Advanced Compartmental Absorption Transit' (ACAT) model and the whole-body disposition model, respectively. Renal and hepatic clearances were simulated using empirical scaling methods followed by incorporation into the WB-PBPK model. Model refinements were conducted after a comparison of the simulated concentration-time profiles and pharmacokinetic parameters with the observed data in healthy adults following intravenous and oral administration. Finally, the WB-PBPK model coupled with a Monte Carlo simulation was employed to predict the mean and variability of bisoprolol pharmacokinetics in virtual healthy subjects and patients. The simulated and observed data after both intravenous and oral dosing showed good agreement for all of the dose levels in the reported normal adult population groups. The predicted pharmacokinetic parameters (AUC, C(max), and T(max)) were reasonably consistent (<1.3-fold error) with the observed values after single oral administration of doses ranging from of 5 to 20 mg using the refined WB-PBPK model. The simulated plasma profiles after multiple oral administration of bisoprolol in healthy adults and patient with renal impairment matched well with the observed profiles. The WB-PBPK model successfully predicts the intravenous and oral pharmacokinetics of bisoprolol across multiple dose levels in diverse normal adult human populations and patients with renal insufficiency.

Novel approach to bioequivalence assessment based on physiologically motivated model

The study was conducted to exemplify an approach capable of obtaining a new insight into bioequivalence (BE) assessment, by the use of a physiologically motivated model. Data from an oral BE study of two piroxicam (PXM) products was used as an example. The BE study was carried out with 24 healthy European subjects according to a two-sequence crossover-randomized design. The test and reference formulations were a PXM generic formulation (LaborMed Pharma, Romania) and Feldene ® (Pfizer, USA), respectively. Plasma concentrations of PXM were monitored by a validated high-performance liquid chromatography over a period of 144 h after administration. After the structure of the optimal model was selected, parameters that characterized the whole-body disposition behavior of PXM in the subjects were derived. The paired Student's t-test and Wilkoxon's test were performed on the derived parameters. The null hypothesis of no differences in the parameters of the whole-body disposition behavior of PXM related to the test and reference product was not rejected at 5% level of significance. This result suggested that the compared products were bioequivalent and could be used interchangeably in clinical setting. The presented approach might show a new way, worth incorporating in future BE guidelines.

Recirculatory model of fentanyl disposition with the brain as the target organ

The factors affecting the concentrations of fentanyl in the brain after intravenous administration have not been completely quantified. A model integrating the role of brain, lung and systemic kinetics was developed based on data from conscious instrumented sheep. Brain kinetics were inferred from arterio-sagittal sinus concentration gradients and cerebral blood flow, and lung kinetics from the pulmonary artery-aortic gradient and cardiac output. The best models of the lung and brain were incorporated into a recirculatory model of the whole-body disposition of fentanyl. The validity of the model structure was tested by its ability to describe published data on the effect of hypo-, normo- and hypercarbia on the blood and brain concentrations of fentanyl in anaesthetized dogs. The cerebral kinetics of fentanyl were consistent with partial membrane limitation: the time to 50% equilibration with arterial blood was 10.0 min. Lung kinetics had two distinct components: a shallow compartment that was 50% equilibrated with blood in 0.72 min, and a loss term probably representing sequestration. Despite its simplicity, the recirculatory model was an adequate description of the sheep data. The dog data could be described if cerebral blood flow and cardiac output in the model were allowed to differ between hypo-, normo- and hypercarbic states. The required flow changes were in good agreement with the known effect of these states in the dog. A recirculatory model with the brain as a target organ defined the quantitative relationship between the brain concentrations of fentanyl and the circulatory state.