Extended Pharmacokinetic Model of the Rabbit Eye for Intravitreal and Intracameral Injections of Macromolecules: Quantitative Analysis of Anterior and Posterior Elimination Pathways

Computer Modeling of Bevacizumab Drug Distribution After Intravitreal Injection in Rabbit and Human Eyes

Age-related macular degeneration (AMD) is a condition brought on by macular deterioration caused primarily by inflammation and cell death in the retina. There is no cure for the disease and current treatments for advanced (wet) AMD rely on intravitreal injections of anti-vascular endothelial growth factor (anti-VEGF) therapeutics. One common off-label anti-VEGF drug used in AMD treatment is bevacizumab. There have been experimental efforts to investigate the pharmacokinetic (PK) behavior of bevacizumab in the vitreous and aqueous humor. Still the quantitative effect of elimination routes and drug concentration in the macula are not well understood. In our study, we developed two spatial models representing rabbit and human vitreous humor to better understand the PK behavior of bevacizumab. We explored convective effects on the vitreous while considering the anterior elimination alone or coupled with posterior elimination. We compared our models with available experimental data and calculated an approximate macula concentration. Our results show that both anterior and posterior elimination play a role in bevacizumab clearance from the eye. Furthermore, an effective bevacizumab concentration close to the macula region is maintained for shorter time periods when compared to the whole vitreous region. This model can improve knowledge and understanding of AMD treatment.

Experimental and mathematical approaches for drug delivery for the treatment of wet age-related macular degeneration

Several chronic eye diseases affect the posterior segment of the eye. Among them age-related macular degeneration can cause vision loss if left untreated and is one of the leading causes of visual impairment in the world. Most treatments are based on intravitreally injected therapeutics that inhibit the action of vascular endothelial growth factor. However, due to the need for monthly injections, this method is associated with poor patient compliance. To address this problem, numerous drug delivery systems (DDSs) have been developed. This review covers a selection of particulate systems, non-stimuli responsive hydrogels, implants, and composite systems that have been developed in the last few decades. Depending on the type of DDS, polymer material, and preparation method, different mechanical properties and drug release profiles can be achieved. Furthermore, DDS development can be optimized by implementing mathematical modeling of both drug release and pharmacokinetic aspects. Several existing mathematical models for diffusion-controlled, swelling-controlled, and erosion-controlled drug delivery from polymeric systems are summarized. Compartmental and physiologically based models for ocular drug transport and pharmacokinetics that have studied drug concentration profiles after intravitreal delivery or release from a DDS are also reviewed. The coupling of drug release models with ocular pharmacokinetic models can lead to obtaining much more efficient DDSs for the treatment of age-related macular degeneration and other diseases of the posterior segment of the eye.

Quantitative intravitreal pharmacokinetics in mouse as a step towards inter-species translation

Although mouse models are widely used in retinal drug development, pharmacokinetics in mouse eye is poorly understood. In this study, we applied non-invasive in vivo fluorophotometry to study pharmacokinetics of intravitreal fluorescein sodium (molecular weight 0.38 kDa) and fluorescein isothiocyanate-dextran (FD-150; molecular weight 150 kDa) in mice. Intravitreal half-lives of fluorescein and FD-150 in mouse eyes were 0.53 ± 0.06 h and 2.61 ± 0.86 h, respectively. These values are 8-230 times shorter than the elimination half-lives of similar compounds in the human vitreous. The apparent volumes of distribution in the mouse vitreous were close to the anatomical volume of the mouse vitreous (FD-150, 5.1 μl; fluorescein, 9.6 μl). Dose scaling factors were calculated from mouse to rat, rabbit, monkey and human translation. Based on pharmacokinetic modelling and compound concentrations in the vitreous and anterior chamber, fluorescein is mainly eliminated posteriorly across blood-retina barrier, but FD-150 is cleared via aqueous humour outflow. The results of this study improve the knowledge of intravitreal pharmacokinetics in mouse and facilitate inter-species scaling in ocular drug development.

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.

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.

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 Model of Macromolecular Drug Transport in a Partially Liquefied Vitreous Humor

The purpose of this study is to investigate the effect of partial liquefaction (due to ageing) of the vitreous humor on the transport of ocular drugs. In our model, the gel part of the vitreous is treated as a Darcy-type porous medium. A spherical region within the porous part of vitreous is in a liquid state which, for computational purposes, is also treated as a porous medium but with a much higher permeability. Using the finite element method, a time-dependent, three-dimensional model has been developed to computationally simulate (using the Petrov-Galerkin method) the transport of intravitreally injected macromolecules where both convection and diffusion are present. From a fluid physics and transport phenomena perspective, the results show many interesting features. For pressure-driven flow across the vitreous, the flow streamlines converge into the liquefied region as the flow seeks the fastest path of travel. Furthermore, as expected, with increased level of liquefaction, the overall flow rate increases for a given pressure drop. We have quantified this effect for various geometrical considerations. The flow convergence into the liquefied region has important implication for convective transport. One effect is the clear diversion of the drug as it reaches the liquefied region. In some instances, the entry point of the drug in the retinal region gets slightly shifted due to liquefaction. While the model has many approximations and assumptions, the focus is illustrating the effect of liquefaction as one of the building blocks toward a fully comprehensive model.

Ocular Fluid Mechanics and Drug Delivery: A Review of Mathematical and Computational Models

The human eye is a complex biomechanical structure with a range of biomechanical processes involved in various physiological as well as pathological conditions. Fluid flow inside different domains of the eye is one of the most significant biomechanical processes that tend to perform a wide variety of functions and when combined with other biophysical processes play a crucial role in ocular drug delivery. However, it is quite difficult to comprehend the effect of these processes on drug transport and associated treatment experimentally because of ethical constraints and economic feasibility. Computational modeling on the other hand is an excellent means to understand the associated complexity between these aforementioned processes and drug delivery. A wide range of computational models specific to different types of fluids present in different domains of the eye as well as varying drug delivery modes has been established to understand the fluid flow behavior and drug transport phenomenon in an insilico manner. These computational models have been used as a non-invasive tool to aid ophthalmologists in identifying the challenges associated with a particular drug delivery mode while treating particular eye diseases and to advance the understanding of the biomechanical behavior of the eye. In this regard, the author attempts to summarize the existing computational and mathematical approaches proposed in the last two decades for understanding the fluid mechanics and drug transport associated with different domains of the eye, together with their application to modify the existing treatment processes.

Extraocular, periocular, and intraocular routes for sustained drug delivery for glaucoma

Although once daily anti-glaucoma drug therapy is a current clinical reality, most therapies require multiple dosing and there is an unmet need to develop convenient, safe, and effective sustained release drug delivery systems for long-term treatment to improve patient adherence and outcomes. One of the first sustained release drug delivery systems was approved for the reduction of intraocular pressure in glaucoma patients. It is a polymeric reservoir-type insert delivery system, Ocusert™, placed under the eyelid and on the ocular surface for zero-order drug release over one week. The insert, marketed in two strengths, released pilocarpine on the eye surface. While many clinicians appreciated this drug product, it was eventually discontinued. No similar sustained release non-invasive drug delivery system has made it to the market to date for treating glaucoma. Drug delivery systems under development include punctal plugs, ring-type systems, contact lenses, implants, microspheres, nanospheres, gels, and other depot systems placed in the extraocular, periocular, or intraocular regions including intracameral, supraciliary, and intravitreal spaces. This article discusses the advantages and disadvantages of the various routes of administration and delivery systems for sustained glaucoma therapy. It also provides the reader with some examples and discussion of drug delivery systems that could potentially be applied for glaucoma treatment. Interestingly, one intracamerally injected implant, Durysta™, was approved recently for sustained intraocular pressure reduction. However, long-term acceptance of such devices has yet to be established. The ultimate success of the delivery system will depend on efficacy relative to eye drop dosing, safety, reimbursement options, and patient acceptance. Cautious development efforts are warranted considering prior failed approaches for sustained glaucoma drug delivery. Neuroprotective approaches for glaucoma therapy including cell, gene, protein, and drug-combination therapies, mostly administered intravitreally, are also rapidly progressing towards assessment in humans.

Immunology and Pathology in Ocular Drug Development

Clear vision is dependent on features that protect the anatomical integrity of the eye (cornea and sclera) and those that contribute to internal ocular homeostasis by conferring hemangiogenic (avascular tissues and antiangiogenic factors), lymphangiogenic (lack of draining lymphatics), and immunologic (tight junctions that form blood-ocular barriers, immunosuppressive cells, and modulators) privileges. The later examples are necessary components that enable the eye to maintain an immunosuppressive environment that responds to foreign invaders in a deviated manner, minimizing destructive inflammation that would impair vision. These conditions allowed for the observations made by Medawar, in 1948, of delayed rejection of allogenic tissue grafts in the anterior chamber of mouse eye and permit the sequestration of foreign invaders (eg, Toxoplasma gondii) within the retina of healthy individuals. Yet successful development of intraocular drugs (biologics and delivery devices) has been stymied by adverse ocular pathology, much of which is driven by immune pathways. The eye can be intolerant of foreign protein irrespective of delivery route, and endogenous ocular cells have remarkable plasticity when recruited to preserve visual function. This article provides a review of current understanding of ocular immunology and the potential role of immune mechanisms in pathology observed with intraocular drug delivery.

Effect of aging on heat transfer, fluid flow and drug transport in anterior human eye: A computational study

There are a lot of geometrical and morphological changes that happen in the human eye with age. Primary open-angle glaucoma, which is caused by the increase in intraocular pressure inside the anterior chamber of the eye is also associated with the physiological aging of the eye. Therefore, it is crucial to understand the effects of aging on drug delivery in the human eye when applied topically. Consequently, a numerical model of topical drug delivery for an aging human eye has been developed using commercial software COMSOL Multiphysics in the current study. Three different age groups (young, middle and old) have been considered and the changes in geometrical and tissue properties of different domains of the eye with age have been included in the numerical model. The effect of aging on heat transfer, aqueous humor flow, intraocular pressure and drug concentration in different domains and orientations of the eye have been investigated. Additionally, an attempt has been made to predict the best class of anti-glaucomatic treatment in silico that should be preferred to treat primary open-angle glaucoma effectively. Results illustrate that there is a decrease in the average corneal temperature and an increase in the temperature deviation across the cornea with age. Further, there is a decrease in the aqueous humor flow magnitude in the anterior chamber of the eye and an increase in intraocular pressure in the anterior chamber of older age groups, which leads to primary open-angle glaucoma. The reduced aqueous humor flow leads to increased drug concentration in the anterior chamber as well as iris and reduced drug concentration in the trabecular mesh of the older age groups, thereby affecting the treatment efficacy. Additionally, our simulated results demonstrate that anti-glaucomatic treatments should be more focused on treating the trabecular mesh rather than the ciliary body of the eye.

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.

Extended Pharmacokinetic Model of the Intravitreal Injections of Macromolecules in Rabbits. Part 2: Parameter Estimation Based on Concentration Dynamics in the Vitreous, Retina, and Aqueous Humor

PURPOSE

To estimate the diffusion coefficients of an IgG antibody (150 kDa) and its antigen-binding fragment (Fab; 50 kDa) in the neural retina (Dret) and the combined retinal pigment epithelium-choroid (DRPE-cho) with a 3-dimensional (3D) ocular pharmacokinetic (PK) model of the rabbit eye.

METHODS

Vitreous, retina, and aqueous humor concentrations of IgG and Fab after intravitreal injection in rabbits were taken from Gadkar et al. (2015). A least-squares method was used to estimate Dret and DRPE-cho with the 3D finite element model where mass transport was defined with diffusion and convection. Different intraocular pressures (IOP), initial distribution volumes (Vinit), and neural retina/vitreous partition coefficients (Kret/vit) were tested. Sensitivity analysis was performed for the final model.

RESULTS

With the final IgG model (IOP 10.1 Torr, Vinit 400 μl, Kret/vit 0.5), the estimated Dret and DRPE-cho were 36.8 × 10-9 cm2s-1 and 4.11 × 10-9 cm2s-1, respectively, and 76% of the dose was eliminated via the anterior chamber. Modeling of Fab revealed that a physiological model parameter "aqueous humor formation rate" sets constraints that need to be considered in the parameter estimation.

CONCLUSIONS

This study extends the use of 3D ocular PK models for parameter estimation using simultaneously macromolecule concentrations in three ocular tissues.

Intravitreal hydrogels for sustained release of therapeutic proteins

This review highlights how hydrogel formulations can improve intravitreal protein delivery to the posterior segment of the eye in order to increase therapeutic outcome and patient compliance. Several therapeutic proteins have shown excellent clinical successes for the treatment of various intraocular diseases. However, drug delivery to the posterior segment of the eye faces significant challenges due to multiple physiological barriers preventing drugs from reaching the retina, among which intravitreal protein instability and rapid clearance from the site of injection. Hence, frequent injections are required to maintain therapeutic levels. Moreover, because the world population ages, the number of patients suffering from ocular diseases, such as age-related macular degeneration (AMD) and diabetic retinopathy (DR) is increasing and causing increased health care costs. Therefore, there is a growing need for suitable delivery systems able to tackle the current limitations in retinal protein delivery, which also may reduce costs. Hydrogels have shown to be promising delivery systems capable of sustaining release of therapeutic proteins and thus extending their local presence. Here, an extensive overview of preclinically developed intravitreal hydrogels is provided with attention to the rational design of clinically useful intravitreal systems. The currently used polymers, crosslinking mechanisms, in vitro/in vivo models and advancements are discussed together with the limitations and future perspective of these biomaterials.

Shape Memory Tube Plug for Fine-control of Intraocular Pressure by Glaucoma Devices

As a standard approach to treating glaucoma, a silicone tube device is inserted into the eye to drain water and thereby lower the intraocular pressure (IOP). However, the fixed diameter of the tube often results in too much initial water drain. Suture-sheaving of the outer tube wall also leads to random IOP drops over time. In this study, self-expansion of the inner tube diameter was achieved by inserting a shape memory tube into the silicone tube. The difference in controlling small to large IOP drops before and after tube diameter expansion was demonstrated via computational modeling, a flow pumping system, and rabbit experiments.

Prodrug Approach for Posterior Eye Drug Delivery: Synthesis of Novel Ganciclovir Prodrugs and in Vitro Screening with Cassette Dosing

Because of poor ocular drug bioavailability, intravitreal injections have become the gold standard for drug delivery to the posterior eye. The prodrug approach can be used for optimizing the biopharmaceutical properties of intravitreal drugs. The preclinical screening of prodrugs' properties, such as hydrolysis and bioconversion, should be conducted in a resource-efficient way for an extensive set of synthesized compounds with validated methods. Our objective was to explore cassette dosing in in vitro prodrug hydrolysis and bioconversion studies in buffer, vitreous, and retinal pigment epithelium (RPE) homogenate for rapid medium-throughput screening. Moreover, our aim was to correlate the prodrug structure with hydrolytic behavior. We synthesized 18 novel ganciclovir prodrugs and first studied their hydrolysis in aqueous buffer and porcine vitreous in vitro with cassette dosing for 35 h. A method for vitreous homogenate pH equilibration to a physiological level by using buffer and incubation under 5% carbon dioxide was validated. The hydrolysis of the prodrugs was evaluated in porcine RPE homogenate in vitro with cassette dosing, and five prodrugs were assayed individually to examine their bioconversion into ganciclovir in RPE after 2 h. Lastly, the prodrugs' binding to melanin was studied in vitro. The prodrugs showed a wide spectrum of hydrolysis rates, ranging from a few percentages to 100% in the vitreous and RPE; in general, hydrolysis in RPE was faster than in vitreous. Prodrugs with long carbon chains and disubstitution showed lability in the tissue homogenates, whereas prodrugs with branched carbon chains and aromatic groups were stable. All five prodrugs chosen for the bioconversion study in RPE were hydrolyzed into ganciclovir, and their hydrolytic behavior matched results from the cassette mix experiment, supporting the cassette mix approach for hydrolysis and bioconversion studies. None of the prodrugs bound highly to melanin (<50% bound). In conclusion, cassette dosing proved useful for the rapid screening of prodrug hydrolysis and bioconversion properties. Analyzing several compounds simultaneously can complicate the analytics, and thus, choosing the compounds of the cassette mix should be done carefully to avoid mutual interference of the compounds with the results. The methodology and results of the work are applicable in ocular drug research and prodrug design.

No flow through the vitreous humor: How strong is the evidence?

When analyzing vitreal drug delivery, or the pharmacological effects of drugs on intraocular pressure, or when interpreting outflow facility measurements, it is generally accepted that the fluid in the vitreous humor is stagnant. It is accepted that for all practical purposes, the aqueous fluid exits the eye via anterior pathways only, and so there is negligible if any posteriorly directed flow of aqueous through the vitreous humor. This assumption is largely based on the interpretation of experimental data from key sources including Maurice (1957), Moseley (1984), Gaul and Brubaker (1986), Maurice (1987) and Araie et al. (1991). However, there is strong independent evidence suggesting there is a substantial fluid flow across the retinal pigment epithelium from key sources including Cantrill and Pederson (1984), Chihara and Nao-i, Tsuboi (1985), Dahrouj et al. (2014), Smith and Gardiner (2017) and Smith et al. (2019). The conflicting evidence creates a conundrum-how can both interpretations be true? This leads us to re-evaluate the evidence. We demonstrate that the data believed to be supporting no aqueous flow through the vitreous are in fact compatible with a significant normal aqueous flow. We identify strong and independent lines of evidence supporting fluid flow across the RPE, including our new outflow model for the eye. On balance it appears the current evidence favors the view that there is normally a significant aqueous flow across the RPE in vivo. This finding suggests that past and future analyses of outflow facility, interpretations of some drug distributions and the interpretation of some drug effects on eye tissues, may need to be revised.

Remodeling of Cellular Surfaces via Fast Disulfide-Thiol Exchange To Regulate Cell Behaviors

Remodeling of cellular surfaces is shown highly effective in the manipulation and control of cell behaviors via nonbiological means. By 5-thio-2-nitrobenzoate-mediated, fast, and reversible disulfide-thiol exchange, a sequential layer by layer assembly process was developed to grow albumin protein shells on cellular surfaces fixed by a disulfide-linked network, in a cytocompatible manner. The artificial shells, accomplished by a double-assembly process, were sustainable up to >1 day, and thereafter gradually bioabsorbed with unaffected cell viability. The surface engineering process enabled dynamic remodeling of cellular surfaces that effectively controlled cell behaviors including regulated cell proliferation, enhanced uptake efficiency of dextran-fluorescein isothiocyanate that is known for cell-impermeability, and targeted imaging. This unique approach was well-validated on tumor cells (B16), immune cells (DC2.4), and neutrophils, showing its potential universality for most of the cells that are rich in thiols. The new strategy will show promise in cell manipulation and targeted imaging.

Hydrodynamics of Intravitreal Injections into Liquid Vitreous Substitutes

Intravitreal injections have become the cornerstone of retinal care and one of the most commonly performed procedures across all medical specialties. The impact of hydrodynamic forces of intravitreal solutions when injected into vitreous or vitreous substitutes has not been well described. While computational models do exist, they tend to underestimate the starting surface area of an injected bolus of a drug. Here, we report the dispersion profile of a dye bolus (50 µL) injected into different vitreous substitutes of varying viscosities, surface tensions, and volumetric densities. A novel 3D printed in vitro model of the vitreous cavity of the eye was designed to visualize the dispersion profile of solutions when injected into the following vitreous substitutes-balanced salt solution (BSS), sodium hyaluronate (HA), and silicone oils (SO)-using a 30G needle with a Reynolds number (Re) for injection ranging from approximately 189 to 677. Larger bolus surface areas were associated with faster injection speeds, lower viscosity of vitreous substitutes, and smaller difference in interfacial surface tensions. Boluses exhibited buoyancy when injected into standard S1000. The hydrodynamic properties of liquid vitreous substitutes influence the initial injected bolus dispersion profile and should be taken into account when simulating drug dispersion following intravitreal injection at a preclinical stage of development, to better inform formulations and performance.

Physiologically based ocular pharmacokinetic modeling using computational methods

By explicitly representing ocular anatomy, computational fluid dynamic simulation methods model drug mass transport both within and between ocular tissue regions, providing reliable animal-to-human translation of bioavailability. Here, we apply physiologically based models to simulate ocular drug administration. A non-anatomical model is used that applies a simple theorem for calculating ocular bioavailability from a topical dose. A computational fluid dynamic model is also described that incorporates ocular physiology in anatomical models for rabbit, monkey and man. This second method applies material properties and boundary conditions for various tissues enabling simulation of fluid flows, pressures, temperatures, convection, and drug advection following various modes of administration. The method provides a regional distribution with a given tissue not available using standard compartmental models, and enables translation of results from animal experiments into predictions for human ocular pharmacokinetics (PK).

Application of Mechanistic Ocular Absorption Modeling and Simulation to Understand the Impact of Formulation Properties on Ophthalmic Bioavailability in Rabbits: a Case Study Using Dexamethasone Suspension

Developing mathematical models to predict changes in ocular bioavailability and pharmacokinetics due to differences in the physicochemical properties of complex topical ophthalmic suspension formulations is important in drug product development and regulatory assessment. Herein, we used published FDA clinical pharmacology review data, in-house, and literature rabbit pharmacokinetic data generated for dexamethasone ophthalmic suspensions to demonstrate how the mechanistic Ocular Compartmental Absorption and Transit model by GastroPlus™ can be used to characterize ocular drug pharmacokinetic performance in rabbits for suspension formulations. This model was used to describe the dose-dependent (0.01 to 0.1%) non-linear pharmacokinetic in ocular tissues and characterize the impact of viscosity (1.67 to 72.9 cP) and particle size (5.5 to 22 μm) on in vivo ocular drug absorption and disposition. Parameter sensitivity analysis (hypothetical suspension particle size: 1 to 10 μm, viscosity: 1 to 100 cP) demonstrated that the interplay between formulation properties and physiological clearance through drainage and tear turnover rates in the pre-corneal compartment drives the ocular drug bioavailability. The quick removal of drug suspended particles from the pre-corneal compartment renders the impact of particle size inconsequential relative to viscosity modification. The in vivo ocular absorption is (1) viscosity non-sensitive when the viscosity is high and the impact of viscosity on the pre-corneal residence time reaches the maximum physiological system capacity or (2) viscosity sensitive when the viscosity is below a certain limit. This study reinforces our understanding of the interplay between physiological factors and ophthalmic formulation physicochemical properties and their impact on in vivo ocular drug PK performance in rabbits.

Pharmacokinetic Simulations of Intravitreal Biologicals: Aspects of Drug Delivery to the Posterior and Anterior Segments

Biologicals are important ocular drugs that are be delivered using monthly and bimonthly intravitreal injections to treat retinal diseases, such as age-related macular degeneration. Long acting delivery systems are needed for prolongation of their dosing interval. Intravitreal biologicals are eliminated from the eye via the aqueous humor outflow. Thus, the anterior and posterior segments are exposed to the drug. We utilized a kinetic simulation model to estimate protein drug concentrations in the vitreous and aqueous humor after bolus injection and controlled release administration to the vitreous. The simulations predicted accurately the experimental levels of 5 biologicals in the vitreous and aqueous humor. The good match between the simulations and experimental data demonstrated almost complete anterior segment bioavailability, and major dose sparing with ocular controlled release systems. Overall, the model is a useful tool in the design of intraocular delivery of biologicals.