Trabecular Meshwork Motion Profile from Pulsatile Pressure Transients: A New Platform to Simulate Transitory Responses in Humans and Nonhuman Primates

Volumetric imaging of trabecular meshwork dynamic motion using 600 kHz swept source optical coherence tomography

The motion of the trabecular meshwork (TM) facilitates the aqueous drainage from the anterior chamber to the venous system, thereby maintaining normal intraocular pressure. As such, characterizing the TM motion is valuable for assessing the functionality of the aqueous outflow system, as demonstrated by previous phase-sensitive optical coherence tomography (OCT) studies. Current methods typically acquire motion from a single cross-sectional plane along the circumference of the anterior chamber. While effective, the lateral scan pattern only intersects one spatial location on the TM at a time, significantly limiting examination throughput. In this study, we introduce the first volumetric imaging approach for assessing TM motion. Rather than monitoring a single cross-sectional plane, our method employs repeated volumetric scans, allowing for simultaneous observation of a continuous TM band spanning two millimeters. We also show that the field of view could be further expanded by stitching multiple scans. To ensure robust data processing, we developed a customized volume registration algorithm to correct motion artifacts and an automated segmentation algorithm to identify the TM boundary based on the correlation of OCT phase dynamics with heartbeats. Imaging results from a healthy subject confirmed the feasibility of our approach, revealing considerable variation in TM motions at different spatial locations through the stitching process. This proposed methodology offers unprecedented capabilities and examination throughput in the biomechanical imaging of the TM, providing significant scientific insights and diagnostic value for identifying abnormalities in aqueous outflow.

Comparative analysis of traction forces in normal and glaucomatous trabecular meshwork cells within a 3D, active fluid-structure interaction culture environment

This study presents a 3D in vitro cell culture model, meticulously 3D printed to replicate the conventional aqueous outflow pathway anatomical structure, facilitating the study of trabecular meshwork (TM) cellular responses under glaucomatous conditions. Glaucoma affects TM cell functionality, leading to extracellular matrix (ECM) stiffening, enhanced cell-ECM adhesion, and obstructed aqueous humor outflow. Our model, reconstructed from polyacrylamide gel with elastic moduli of 1.5 and 21.7 kPa, is based on serial block-face scanning electron microscopy images of the outflow pathway. It allows for quantifying 3D, depth-dependent, dynamic traction forces exerted by both normal and glaucomatous TM cells within an active fluid-structure interaction (FSI) environment. In our experimental design, we designed two scenarios: a control group with TM cells observed over 20 hours without flow (static setting), focusing on intrinsic cellular contractile forces, and a second scenario incorporating active FSI to evaluate its impact on traction forces (dynamic setting). Our observations revealed that active FSI results in higher traction forces (normal: 1.83-fold and glaucoma: 2.24-fold) and shear strains (normal: 1.81-fold and glaucoma: 2.41-fold), with stiffer substrates amplifying this effect. Glaucomatous cells consistently exhibited larger forces than normal cells. Increasing gel stiffness led to enhanced stress fiber formation in TM cells, particularly in glaucomatous cells. Exposure to active FSI dramatically altered actin organization in both normal and glaucomatous TM cells, particularly affecting cortical actin stress fiber arrangement. This model while preliminary offers a new method in understanding TM cell biomechanics and ECM stiffening in glaucoma, highlighting the importance of FSI in these processes. STATEMENT OF SIGNIFICANCE: This pioneering project presents an advanced 3D in vitro model, meticulously replicating the human trabecular meshwork's anatomy for glaucoma research. It enables precise quantification of cellular forces in a dynamic fluid-structure interaction, a leap forward from existing 2D models. This advancement promises significant insights into trabecular meshwork cell biomechanics and the stiffening of the extracellular matrix in glaucoma, offering potential pathways for innovative treatments. This research is positioned at the forefront of ocular disease study, with implications that extend to broader biomedical applications.

Segmental biomechanics of the normal and glaucomatous human aqueous outflow pathway

The conventional aqueous outflow pathway, encompassing the trabecular meshwork (TM), juxtacanalicular connective tissue (JCT), and inner wall endothelium of Schlemm's canal (SC), governs intraocular pressure (IOP) regulation. This study targets the biomechanics of low-flow (LF) and high-flow (HF) regions within the aqueous humor outflow pathway in normal and glaucomatous human donor eyes, using a combined experimental and computational approach. LF and HF TM/JCT/SC complex tissues from normal and glaucomatous eyes underwent uniaxial tensile testing. Dynamic motion of the TM/JCT/SC complex was recorded using customized green-light optical coherence tomography during SC pressurization in cannulated anterior segment wedges. A hyperviscoelastic model quantified TM/JCT/SC complex properties. A fluid-structure interaction model simulated tissue-aqueous humor interaction. FluoSpheres were introduced into the pathway via negative pressure in the SC, with their motion tracked using two-photon excitation microscopy. Tensile test results revealed that the elastic moduli of the LF and HF regions in glaucomatous eyes are 3.5- and 1.5-fold stiffer than the normal eyes, respectively. The FE results also showed significantly larger shear moduli in the TM, JCT, and SC of the glaucomatous eyes compared to the normal subjects. The LF regions in normal eyes demonstrated larger elastic moduli compared to the HF regions in glaucomatous eyes. The resultant strain in the outflow tissues and velocity of the aqueous humor in the FSI models were in good agreement with the digital volume correlation and 3D particle image velocimetry data, respectively. This study uncovers stiffer biomechanical responses in glaucomatous eyes, with LF regions stiffer than HF regions in both normal and glaucomatous eyes. STATEMENT OF SIGNIFICANCE: This study delves into the biomechanics of the conventional aqueous outflow pathway, a crucial regulator of intraocular pressure and ocular health. By analyzing mechanical differences in low-flow and high-flow regions of normal and glaucomatous eyes, this research unveils the stiffer response in glaucomatous eyes. The distinction between regions' properties offers insights into aqueous humor outflow regulation, while the integration of experimental and computational methods enhances credibility. These findings have potential implications for disease management and present a vital step toward innovative ophthalmic interventions. This study advances our understanding of glaucoma's biomechanical basis and its broader impact on ocular health.

Developing an experimental-computational workflow to study the biomechanics of the human conventional aqueous outflow pathway

The aqueous humor actively interacts with the trabecular meshwork (TM), juxtacanalicular tissue (JCT), and Schlemm's canal (SC) through a dynamic fluid-structure interaction (FSI) coupling. Despite the fact that intraocular pressure (IOP) undergoes significant fluctuations, our understanding of the hyperviscoelastic biomechanical properties of the aqueous outflow tissues is limited. In this study, a quadrant of the anterior segment from a normal human donor eye was dynamically pressurized in the SC lumen, and imaged using a customized optical coherence tomography (OCT). The TM/JCT/SC complex finite element (FE) with embedded collagen fibrils was reconstructed based on the segmented boundary nodes in the OCT images. The hyperviscoelastic mechanical properties of the outflow tissues' extracellular matrix with embedded viscoelastic collagen fibrils were calculated using an inverse FE-optimization method. Thereafter, the 3D microstructural FE model of the TM, with adjacent JCT and SC inner wall, from the same donor eye was constructed using optical coherence microscopy and subjected to a flow load-boundary from the SC lumen. The resultant deformation/strain in the outflow tissues was calculated using the FSI method, and compared to the digital volume correlation (DVC) data. TM showed larger shear modulus (0.92 MPa) compared to the JCT (0.47 MPa) and SC inner wall (0.85 MPa). Shear modulus (viscoelastic) was larger in the SC inner wall (97.65 MPa) compared to the TM (84.38 MPa) and JCT (56.30 MPa). The conventional aqueous outflow pathway is subjected to a rate-dependent IOP load-boundary with large fluctuations. This necessitates addressing the biomechanics of the outflow tissues using hyperviscoelastic material-model. STATEMENT OF SIGNIFICANCE: While the human conventional aqueous outflow pathway is subjected to a large-deformation and time-dependent IOP load-boundary, we are not aware of any studies that have calculated the hyperviscoelastic mechanical properties of the outflow tissues with embedded viscoelastic collagen fibrils. A quadrant of the anterior segment of a normal humor donor eye was dynamically pressurized from the SC lumen with relatively large fluctuations. The TM/JCT/SC complex were OCT imaged and the mechanical properties of the tissues with embedded collagen fibrils were calculated using the inverse FE-optimization algorithm. The resultant displacement/strain in the FSI outflow model was validated versus the DVC data. The proposed experimental-computational workflow may significantly contribute to understanding of the effects of different drugs on the biomechanics of the conventional aqueous outflow pathway.

Morphological and biomechanical analyses of the human healthy and glaucomatous aqueous outflow pathway: Imaging-to-modeling

BACKGROUND AND OBJECTIVE

Intraocular pressure (IOP) is maintained via a dynamic balance between the production of aqueous humor and its drainage through the trabecular meshwork (TM), juxtacanalicular connective tissue (JCT), and Schlemm's canal (SC) endothelium of the conventional outflow pathway. Primary open angle glaucoma (POAG) is often associated with IOP elevation that occurs due to an abnormally high outflow resistance across the outflow pathway. Outflow tissues are viscoelastic and actively interact with aqueous humor dynamics through a two-way fluid-structure interaction coupling. While glaucoma affects the morphology and stiffness of the outflow tissues, their biomechanics and hydrodynamics in glaucoma eyes remain largely unknown. This research aims to develop an image-to-model method allowing the biomechanics and hydrodynamics of the conventional aqueous outflow pathway to be studied.

METHODS

We used a combination of X-ray computed tomography and scanning electron microscopy to reconstruct high-fidelity, eye-specific, 3D microstructural finite element models of the healthy and glaucoma outflow tissues in cellularized and decellularized conditions. The viscoelastic TM/JCT/SC complex finite element models with embedded viscoelastic beam elements were subjected to a physiological IOP load boundary; the stresses/strains and the flow state were calculated using fluid-structure interaction and computational fluid dynamics.

RESULTS

Based on the resultant hydrodynamics parameters across the outflow pathway, the primary site of outflow resistance in healthy eyes was in the JCT and immediate vicinity of the SC inner wall, while the majority of the outflow resistance in the glaucoma eyes occurred in the TM. The TM and JCT in the glaucoma eyes showed 1.32-fold and 1.13-fold larger beam thickness and smaller trabecular space size (2.24-fold and 1.50-fold) compared to the healthy eyes.

CONCLUSIONS

Characterizing the accurate morphology of the outflow tissues may significantly contribute to constructing more accurate, robust, and reliable models, that can eventually help to better understand the dynamic IOP regulation, hydrodynamics of the aqueous humor, and outflow resistance dynamic in the human eyes. This model demonstrates proof of concept for determining changes to outflow resistance in healthy and glaucomatous tissues and thus may be utilized in larger cohorts of donor tissues where disease specificity, race, age, and gender of the eye donors may be accounted for.

The Effect of Intraocular Pressure Load Boundary on the Biomechanics of the Human Conventional Aqueous Outflow Pathway

BACKGROUND

Aqueous humor outflow resistance in the trabecular meshwork (TM), juxtacanalicular connective tissue (JCT), and Schlemm's canal (SC) endothelium of the conventional outflow pathway actively contribute to intraocular pressure (IOP) regulation. Outflow resistance is actively affected by the dynamic outflow pressure gradient across the TM, JCT, and SC inner wall tissues. The resistance effect implies the presence of a fluid-structure interaction (FSI) coupling between the outflow tissues and the aqueous humor. However, the biomechanical interactions between viscoelastic outflow tissues and aqueous humor dynamics are largely unknown.

METHODS

A 3D microstructural finite element (FE) model of a healthy human eye TM/JCT/SC complex was constructed with elastic and viscoelastic material properties for the bulk extracellular matrix and embedded elastic cable elements. The FE models were subjected to both idealized and a physiologic IOP load boundary using the FSI method.

RESULTS

The elastic material model for both the idealized and physiologic IOP load boundary at equal IOPs showed similar stresses and strains in the outflow tissues as well as pressure in the aqueous humor. However, outflow tissues with viscoelastic material properties were sensitive to the IOP load rate, resulting in different mechanical and hydrodynamic responses in the tissues and aqueous humor.

CONCLUSIONS

Transient IOP fluctuations may cause a relatively large IOP difference of ~20 mmHg in a very short time frame of ~0.1 s, resulting in a rate stiffening in the outflow tissues. Rate stiffening reduces strains and causes a rate-dependent pressure gradient across the outflow tissues. Thus, the results suggest it is necessary to use a viscoelastic material model in outflow tissues that includes the important role of IOP load rate.

Viscoelastic Biomechanical Properties of the Conventional Aqueous Outflow Pathway Tissues in Healthy and Glaucoma Human Eyes

BACKGROUND

Although the tissues comprising the ocular conventional outflow pathway have shown strong viscoelastic mechanical response to aqueous humor pressure dynamics, the viscoelastic mechanical properties of the trabecular meshwork (TM), juxtacanalicular connective tissue (JCT), and Schlemm's canal (SC) inner wall are largely unknown.

METHODS

A quadrant of the anterior segment from two human donor eyes at low- and high-flow (LF and HF) outflow regions was pressurized and imaged using optical coherence tomography (OCT). A finite element (FE) model of the TM, the adjacent JCT, and the SC inner wall was constructed and viscoelastic beam elements were distributed in the extracellular matrix (ECM) of the TM and JCT to represent anisotropic collagen. An inverse FE-optimization algorithm was used to calculate the viscoelastic properties of the ECM/beam elements such that the TM/JCT/SC model and OCT imaging data best matched over time.

RESULTS

The ECM of the glaucoma tissues showed significantly larger time-dependent shear moduli compared to the heathy tissues. Significantly larger shear moduli were also observed in the LF regions of both the healthy and glaucoma eyes compared to the HF regions.

CONCLUSIONS

The outflow tissues in both glaucoma eyes and HF regions are stiffer and less able to respond to dynamic IOP.

Biomechanics of human trabecular meshwork in healthy and glaucoma eyes via dynamic Schlemm's canal pressurization

BACKGROUND AND OBJECTIVE

The trabecular meshwork (TM) consists of extracellular matrix (ECM) with embedded collagen and elastin fibers providing its mechanical support. TM stiffness is considerably higher in glaucoma eyes. Emerging data indicates that the TM moves dynamically with transient intraocular pressure (IOP) fluctuations, implying the viscoelastic mechanical behavior of the TM. However, little is known about TM viscoelastic behavior. We calculated the viscoelastic mechanical properties of the TM in n = 2 healthy and n = 2 glaucoma eyes.

METHODS

A quadrant of the anterior segment was submerged in a saline bath, and a cannula connected to an adjustable saline reservoir was inserted into Schlemm's canal (SC). A spectral domain-OCT (SD-OCT) provided continuous cross-sectional B-scans of the TM/JCT/SC complex during pressure oscillation from 0 to 30 mmHg at two locations. The TM/JCT/SC complex boundaries were delineated to construct a 20-µm-thick volume finite element (FE) mesh. Pre-tensioned collagen and elastin fibrils were embedded in the model using a mesh-free penalty-based cable-in-solid algorithm. SC pressure was represented by a position- and time-dependent pressure boundary; floating boundary conditions were applied to the other cut edges of the model. An FE-optimization algorithm was used to adjust the ECM/fiber mechanical properties such that the TM/JCT/SC model and SD-OCT imaging data best matched over time.

RESULTS

Significantly larger short- and long-time ECM shear moduli (p = 0.0032), and collagen (1.82x) and elastin (2.72x) fibril elastic moduli (p = 0.0001), were found in the TM of glaucoma eyes compared to healthy controls.

CONCLUSIONS

These findings provide additional clarity on the mechanical property differences in healthy and glaucomatous outflow pathway under dynamic loading. Understanding the viscoelastic properties of the TM may serve as a new biomarker in early diagnosis of glaucoma.

Valve-Like Outflow System Behavior With Motion Slowing in Glaucoma Eyes: Findings Using a Minimally Invasive Glaucoma Surgery–MIGS-Like Platform and Optical Coherence Tomography Imaging

Purpose

This study aimed to investigate anatomic relationships and biomechanics of pressure-dependent trabecular meshwork and distal valve-like structure deformation in normal and glaucoma eyes using high-resolution optical coherence tomography (HR-OCT).

Methods

We controlled Schlemm’s canal (SC) pressure during imaging with HR-OCT in segments of three normal (NL) and five glaucomatous (GL) ex vivo eyes. The dissected limbal wedges were studied from 15 locations (5 NL and 10 GL). A minimally invasive glaucoma surgery (MIGS)-like cannula was inserted into the SC lumen, whereas the other end was attached to a switch between two reservoirs, one at 0, the other at 30 mm Hg. A steady-state pressure of 30 mm Hg was maintained to dilate SC and collector channels (CC) during 3D volume imaging. The resulting 3D lumen surface relationships were correlated with internal structural features using an image mask that excluded tissues surrounding SC and CC. While imaging with HR-OCT, real-time motion responses in SC and CC areas were captured by switching pressure from 0 to 30 or 30 to 0 mm Hg. NL vs. GL motion differences were compared.

Results

Lumen surface and internal relationships were successfully imaged. We identified SC inlet and outlet valve-like structures. In NL and GL, the mean SC areas measured at the steady-state of 0 and 30 mm Hg were each significantly different (p < 0.0001). Synchronous changes in SC and CC lumen areas occurred in <200 ms. Measured SC area differences at the steady-state 0 and 30 mmHg, respectively, were larger in NL than GL eyes (p < 0.0001). The SC motion curves rose significantly more slowly in GL than NL (p < 0.001). Pressure waves traveled from the cannula end along the SC lumen to CC and deep intrascleral channels.

Conclusion

HR-OCT provided simultaneous measurements of outflow pathway lumen surfaces, internal structures, and biomechanics of real-time pressure-dependent dimension changes. We identified SC inlet and outlet valve-like structures. GL tissues underwent less motion and responded more slowly than NL, consistent with increased tissue stiffness. A MIGS-like shunt to SC permitted pulse waves to travel distally along SC lumen and into CC.