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.

Fluid mixing operations are central to possibly all chemical, petrochemical, and pharmaceutical industries either being related to biphasic blending in polymerisation processes, cell suspension for biopharmaceuticals production, and fractionation of complex oil mixtures. This work aims at providing a fundamental understanding of the mixing and stretching dynamics occurring in a reactor in the presence of a vortical structure, and the vortex ring was selected as a flow paradigm of vortices commonly encountered in stirred and shaken reactors in laminar flow conditions. High resolution laser induced fluorescence and particle imaging velocimetry measurements were carried out to fully resolve the flow dissipative scales and provide a complete data set to fully assess macro- and micro-mixing characteristics. The analysis builds upon the Lamb-Oseen vortex work of Meunier and Villermaux [“How vortices mix,” J. Fluid Mech. 476, 213–222 (2003)] and the engulfment model of Baldyga and Bourne [“Simplification of micromixing calculations. I. Derivation and application of new model,” Chem. Eng. J. 42, 83–92 (1989); “Simplification of micromixing calculations. II. New applications,” ibid. 42, 93–101 (1989)] which are valid for diffusion-free conditions, and a comparison is made between three methodologies to assess mixing characteristics. The first method is commonly used in macro-mixing studies and is based on a control area analysis by estimating the variation in time of the concentration standard deviation, while the other two are formulated to provide an insight into local segregation dynamics, by either using an iso-concentration approach or an iso-concentration gradient approach to take into account diffusion.

New in vitro prototypes (PK-Eye™) were tested with and without eye movement to understand diffusion and convection effects on intraocular clearance. Port placement in front ((i) ciliary inflow model) and behind the model lens ((ii) posterior inflow model) was used to study bevacizumab (1.25 mg/50 µL) and dexamethasone (0.1 mg/100 µL) in phosphate-buffered saline (PBS, pH 7.4) and simulated vitreal fluid (SVF). Dexamethasone was studied in a (iii) retinal-choroid-sclera (RCS) outflow model (with ciliary inflow and two outflow pathways). Ciliary vs. posterior inflow placement did not affect the half-life for dexamethasone at 2.0 µL/min using PBS (4.7 days vs. 4.8 days) and SVF (4.9 days with ciliary inflow), but it did decrease the half-life for bevacizumab in PBS (20.4 days vs. 2.4 days) and SVF (19.2 days vs. 10.8 days). Eye movement only affected the half-life of dexamethasone in both media. Dexamethasone in the RCS model showed approximately 20% and 75% clearance from the RCS and anterior outflows, respectively. The half-life of the protein was comparable to human data in the posterior inflow model. Shorter half-life values for a protein in a ciliary inflow model can be achieved with other eye movements. The RCS flow model with eye movement was comparable to human half-life data for dexamethasone.

Real-time measurement is important in modern dissolution testing to aid in parallel drug characterisation and quality control (QC). The development of a real-time monitoring platform (microfluidic system, a novel eye movement platform with temperature sensors and accelerometers and a concentration probe setup) in conjunction with an in vitro model of the human eye (PK-Eye™) is reported. The importance of surface membrane permeability when modelling the PK-Eye™ was determined with a “pursing model” (a simplified setup of the hyaloid membrane). Parallel microfluidic control of PK-Eye™ models from a single source of pressure was performed with a ratio of 1:6 (pressure source:models) demonstrating scalability and reproducibility of pressure-flow data. Pore size and exposed surface area helped obtain a physiological range of intraocular pressure (IOP) within the models, demonstrating the need to reproduce in vitro dimensions as closely as possible to the real eye. Variation of aqueous humour flow rate throughout the day was demonstrated with a developed circadian rhythm program. Capabilities of different eye movements were programmed and achieved with an in-house eye movement platform. A concentration probe recorded the real-time concentration monitoring of injected albumin-conjugated Alexa Fluor 488 (Alexa albumin), which displayed constant release profiles. These results demonstrate the possibility of real-time monitoring of a pharmaceutical model for preclinical testing of ocular formulations.

Clinical grading systems following glaucoma filtration surgery do not include any effects of the bleb on the intra-ocular pressure and are relatively subjective, therefore carrying the risk of inter and/or intra-observer variability. The main objective of the study is to quantify and model the effect of subconjunctival bleb on flow pressure for assessment of clinical grading following glaucoma surgery. Subconjunctival bleb was created by inserting a tube into ex vivo rabbit eyes via an ab externo approach through the anterior chamber and exiting into the subconjunctival space. Sterile dyed water was injected through the tube into the developing bleb. For the in vitro approach a silicone bleb was created by clamping a circular silicone sheet, injecting dyed water through a fixed resistance outlet tube. Photographic measurements of the bleb height, planform area and pressure were taken as a function of time. Clinical blebs were also collected over a few months. Mathematical algorithm software was used to build the bleb model. Bleb height and volume increase as pressure in the bleb increases. The bleb planform area tended to a constant determined by the section of conjunctiva prior to shunt insertion. These increases were in accordance with the bleb model developed in the Appendix. They show that the pressure in the bleb is related to the resistance of the outflow. The linearity of clinical grading systems is reviewed and a new grading approach is proposed. The pressure in the bleb has a strong dependence on bleb extent, height and a weak dependence on conjunctival thickness. The pressure in a bleb can be estimated from bleb height, radius, and flow rate inlet in agreement with the bleb flow model. These results provide support for an improved bleb categorization system.