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

Physiologically based pharmacokinetic models for systemic disposition of protein therapeutics in rabbits

Physiologically based pharmacokinetic (PBPK) modelling is an important tool to predict drug disposition in the body. Rabbits play a pivotal role as a highly valued small animal model, particularly in the field of ocular therapeutics, where they serve as a crucial link between preclinical research and clinical applications. In this context, we have developed PBPK models designed specifically for rabbits, with a focus on accurately predicting the pharmacokinetic profiles of protein therapeutics following intravenous administration. Our goal was to comprehend the influence of key physiological factors on systemic disposition of antibodies and their functional derivatives. For the development of the systemic PBPK models, rabbit physiological factors such as gene expression, body weight, neonatal fragment crystallizable receptor (FcRn) binding, target binding, target concentrations, and target turnover rate were meticulously considered. Additionally, key protein parameters, encompassing hydrodynamic radius, binding kinetic constants (KD, koff), internal degradation of the protein-target complex, and renal clearance, were represented in the models. Our final rabbit models demonstrated a robust correlation between predicted and observed serum concentration-time profiles after single intravenous administration in rabbits, covering IgG, Fab, F(ab)2, Fc, and Fc fusion proteins from various publications. These pharmacokinetic simulations offer a promising platform for translating preclinical findings to clinical settings. The presented rabbit intravenous PBPK models lay an important foundation for more specific applications of protein therapeutics in ocular drug development.

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.

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.

Prospective approaches to gene therapy computational modeling - spotlight on viral gene therapy

Clinical studies have found there still exists a lack of gene therapy dose-toxicity and dose-efficacy data that causes gene therapy dose selection to remain elusive. Model informed drug development (MIDD) has become a standard tool implemented throughout the discovery, development, and approval of pharmaceutical therapies, and has the potential to inform dose-toxicity and dose-efficacy relationships to support gene therapy dose selection. Despite this potential, MIDD approaches for gene therapy remain immature and require standardization to be useful for gene therapy clinical programs. With the goal to advance MIDD approaches for gene therapy, in this review we first provide an overview of gene therapy types and how they differ from a bioanalytical, formulation, route of administration, and regulatory standpoint. With this biological and regulatory background, we propose how MIDD can be advanced for AAV-based gene therapies by utilizing physiological based pharmacokinetic modeling and quantitative systems pharmacology to holistically inform AAV and target protein dynamics following dosing. We discuss how this proposed model, allowing for in-depth exploration of AAV pharmacology, could be the key the field needs to treat these unmet disease populations.

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.

Pharmacokinetics of Monoclonal Antibody and Antibody Fragments in The Mouse Eye Following Intravitreal Administration

Mice are rarely used in pharmacokinetic (PK) studies of ocular therapeutics due to the small size of their eyes and challenges in drug administration, tissue collection, and analysis of drug concentrations. Therefore, ocular PK of protein therapeutics in mouse eye following intravitreal (IVT) administration is not known. Here, we have presented the first of its kind investigation, to study the PK of 4 different size non-binding protein therapeutics in mouse plasma, cornea/ICB, vitreous humor, retina, and posterior cup (including choroid) following IVT administration. Administered proteins include trastuzumab (150 kDa) and F(ab)2 (100 kDa), Fab, and scFv (27 kDa) fragments of trastuzumab. An imaging and injection apparatus suitable for performing small (50 nL) IVT injections in mice was developed, and techniques for enucleation of the eye and dissection of ocular tissues were developed. Furthermore, a sensitive enzyme-linked immunosorbent assay (ELISA) for detection of proteins in very small amounts of ocular tissues were developed. It was observed that elimination from the vitreous chamber was the primary driver of PK in the cornea/ICB, retina, posterior cup, and plasma. Trastuzumab displays first-order kinetics in the vitreous humor with a half-life of 18.8 h. F(ab)2, Fab, and ScFv show biphasic PK profiles with distribution phases becoming more rapid as molecular weight decreases, and terminal elimination becoming longer as molecular weight decreases, with terminal half-lives of 16.3, 20.6, and 48.9 h, respectively. The mean residence times of trastuzumab, F(ab)2, Fab, and scFv in the vitreous humor were 26.0, 12.2, 10.7, and 8.16 h, respectively. It was found that the mean residence time in vitreous humor doubles with an increase in molecular weight of ∼69 kDa. Interestingly, the PK of proteins measured in the un-injected eye suggest the presence of a pathway for drug transfer between the eyes, which needs to be further validated. Overall, the findings presented here pave the way for drug discovery and development studies of protein therapeutics for ophthalmic indications in mice.

How can machine learning and multiscale modeling benefit ocular drug development?

The eyes possess sophisticated physiological structures, diverse disease targets, limited drug delivery space, distinctive barriers, and complicated biomechanical processes, requiring a more in-depth understanding of the interactions between drug delivery systems and biological systems for ocular formulation development. However, the tiny size of the eyes makes sampling difficult and invasive studies costly and ethically constrained. Developing ocular formulations following conventional trial-and-error formulation and manufacturing process screening procedures is inefficient. Along with the popularity of computational pharmaceutics, non-invasive in silico modeling & simulation offer new opportunities for the paradigm shift of ocular formulation development. The current work first systematically reviews the theoretical underpinnings, advanced applications, and unique advantages of data-driven machine learning and multiscale simulation approaches represented by molecular simulation, mathematical modeling, and pharmacokinetic (PK)/pharmacodynamic (PD) modeling for ocular drug development. Following this, a new computer-driven framework for rational pharmaceutical formulation design is proposed, inspired by the potential of in silico explorations in understanding drug delivery details and facilitating drug formulation design. Lastly, to promote the paradigm shift, integrated in silico methodologies were highlighted, and discussions on data challenges, model practicality, personalized modeling, regulatory science, interdisciplinary collaboration, and talent training were conducted in detail with a view to achieving more efficient objective-oriented pharmaceutical formulation design.

Quantitative systems pharmacology of the eye: Tools and data for ocular QSP

Good eyesight belongs to the most-valued attributes of health, and diseases of the eye are a significant healthcare burden. Case numbers are expected to further increase in the next decades due to an aging society. The development of drugs in ophthalmology, however, is difficult due to limited accessibility of the eye, in terms of drug administration and in terms of sampling of tissues for drug pharmacokinetics (PKs) and pharmacodynamics (PDs). Ocular quantitative systems pharmacology models provide the opportunity to describe the distribution of drugs in the eye as well as the resulting drug-response in specific segments of the eye. In particular, ocular physiologically-based PK (PBPK) models are necessary to describe drug concentration levels in different regions of the eye. Further, ocular effect models using molecular data from specific cellular systems are needed to develop dose-response correlations. We here describe the current status of PK/PBPK as well as PD models for the eyes and discuss cellular systems, data repositories, as well as animal models in ophthalmology. The application of the various concepts is highlighted for the development of new treatments for postoperative fibrosis after glaucoma surgery.

Smart design of patient centric long-acting products: from preclinical to marketed pipeline trends and opportunities

INTRODUCTION

We see a development in the field of long-acting products to serve patients with chronic diseases by providing benefits in adherence, efficacy and safety of the treatment. This review investigates features of long-acting products on the market/pipeline to understand which drug substance (DS) and drug product (DP) characteristics likely enable a successful patient-centric, low-dosing frequency product.

AREAS COVERED

This review evaluates marketed/pipeline long-acting products with greater than one week release of small molecules and peptides by oral and injectable route of administration (RoA), with particular focus on patient centricity, adherence impact, health outcomes, market trends, and the match of DS/DP technologies which lead to market success.

EXPERT OPINION

Emerging trends are expected to change the field of long-acting products in the upcoming years by increasing capability in engineered molecules (low solubility, long half-life, high potency, etc.), directly developing DP as long-acting oral/injectable, increasing the proportion of products for local drug delivery, and a direction towards more subcutaneous, self-administered products. Among long-acting injectable products, nanosuspensions show a superiority in dose per administration and dosing interval, overwhelming the field of infectious diseases with the recently marketed products.

Quantitative pharmacokinetic analyses of anterior and posterior elimination routes of intravitreal anti-VEGF macromolecules using published human and rabbit data

The purpose of this study was to evaluate the contribution of the anterior elimination route for four anti-vascular endothelial growth factor (anti-VEGF) macromolecules (aflibercept, bevacizumab, pegaptanib and ranibizumab) after intravitreal injection using published human and rabbit data and three previously described pharmacokinetic (PK) modeling methods. A PubMed search was used to identify published studies with concentration-time data. The data were utilized only if the intravitreally injected drugs were used as plain solutions and several criteria for a well-performed PK study were fulfilled. The three methods to analyze rabbit data were (1) the equation for vitreal elimination half-life based molecular size assuming anterior elimination, (2) Maurice equation and plot for the ratio of aqueous humor (AH) to vitreal concentration assuming anterior elimination, and (3) the equation for amount of macromolecule eliminated anteriorly based on the area under the curve in AH. The first and third methods were used for human data. In the second and third methods, AH flow rate is a key model parameter, and it was varied between 2 and 3 μl/min. The methods were applied to data from 9 rabbit studies (1 for aflibercept, 5 for bevacizumab, and 3 for ranibizumab) and 5 human studies (1 for aflibercept, 3 for bevacizumab, and 1 for ranibizumab). Experimental half-lives of anti-VEGF macromolecules in both vitreous and aqueous humor were close to those calculated with the equations for vitreal elimination half-life in humans and rabbits. Rabbit data analyzed with Maurice plot indicated that the contribution of anterior elimination was usually at least 75%. In most human and rabbit studies, the calculated percentage of anterior elimination was at least 51%. Variability between studies was extensive for bevacizumab and ranibizumab. The results suggest that the anterior elimination route dominates after intravitreal injection of anti-VEGF macromolecules. However, the clinical data are sparse and variability is extensive, the latter emphasizing the need of proper experimental design.

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.

Blood retinal barrier and ocular pharmacokinetics: Considerations for the development of oncology drugs

Tyrosine kinase inhibitors (TKIs) are an example of targeted drug therapy to treat cancer while minimizing damage to healthy tissue. In contrast to traditional oncology drugs, the toxicity profile of targeted therapies is less well understood and can include severe ocular adverse events, which are among the most common toxicity reported by these therapeutics. Inhibition of Mer receptor tyrosine kinase (MERTK) promotes innate tumor immunity by decreasing M2-macrophage polarization and efferocytosis. This mechanism offers the opportunity for targeted immunotherapy to treat cancer; however, the ocular expression of MERTK increases the difficulty for developing a targeted drug due to toxicity concerns. In this article we review the pharmacokinetic (PK) parameters and in vitro absorption, distribution, metabolism, and excretion (ADME) assays available to evaluate ocular disposition and assess the relationship between clinical PK and reported ocular events for TKIs to allow backtranslation to preclinical models. Understanding the ocular disposition in the context of PK and safety remains an evolving area and is likely to be a key aspect of developing safe and efficacious oncology drugs, devoid of ocular toxicity.

Ocular Physiologically Based Pharmacokinetic Modeling for Ointment Formulations

Purpose

The purpose of this study is to show how the Ocular Compartmental Absorption & Transit (OCAT™) model in GastroPlus^® can be used to characterize ocular drug pharmacokinetic performance in rabbits for ointment formulations.

Methods

A newly OCAT™ model developed for fluorometholone, as well as a previously verified model for dexamethasone, were used to characterize the aqueous humor (AH) concentration following the administration of multiple ointment formulations to rabbit. The model uses the following parameters: application surface area (SA), a fitted application time, and the fitted Higuchi release constant to characterize the rate of passage of the active pharmaceutical ingredient from the ointment formulations into the tears in vivo.

Results

Parameter sensitivity analysis was performed to understand the impact of ointment formulation changes on ocular exposure. While application time was found to have a significant impact on the time of maximal concentration in AH, both the application SA and the Higuchi release constant significantly influenced both the maximum concentration and the ocular exposure.

Conclusions

This initial model for ointment ophthalmic formulations is a first step to better understand the interplay between physiological factors and ophthalmic formulation physicochemical properties and their impact on in vivo ocular drug pharmacokinetic performance in rabbits.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11095-020-02965-y.