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

Finite Element Analysis of Cervical Spine Kinematic Response during Ejection Utilising a Hill-Type Dynamic Muscle Model

To determine the impact of active muscle on the dynamic response of a pilot's neck during simulated emergency ejection, a detailed three-dimensional (3D) cervical spine (C0-T1) finite element (FE) model integrated with active muscles was constructed. Based on the Hill-type model characterising the muscle force activation mechanics, 13 major neck muscles were modelled. The active force generated by each muscle was simulated as functions of (i) active state (Na), (ii) velocity (Fv(v)), and (iii) length (FL(L)). An acceleration-time profile with an initial acceleration rate of 125 G·s-1 in the 0-80 ms period, reaching peak acceleration of 10 G, then kept constant for a further 70 ms, was applied. The rotational angles of each cervical segment under these ejection conditions were compared with those without muscles and with passive muscles derived from the previous study. Similar trends of segmental rotation were observed with S- and C-curvature of the cervical spine in the 150 ms span analysed. With active muscles, the flexion motion of the C0-C2 segments exhibited higher magnitudes of rotation compared to those without muscle and passive muscle models. The flexion motion increased rapidly and peaked at about 95-105 ms, then decreased rapidly to a lower magnitude. Lower C2-T1 segments exhibited less variation in flexion and extension motions. Overall, during emergency ejections, active muscle activities effectively reduce the variability in rotational angles across cervical segments, except C0-C2 segments in the 60-120 ms period. The role of the active state dynamics of the muscles was crucial to the magnitude of the muscle forces demonstrated. This indicates that it is crucial for pilots to consciously contract their muscles before ejection to prevent cervical spine injuries.

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.

Dynamic traction force in trabecular meshwork cells: A 2D culture model for normal and glaucomatous states

Glaucoma, which is associated with intraocular pressure (IOP) elevation, results in trabecular meshwork (TM) cellular dysfunction, leading to increased rigidity of the extracellular matrix (ECM), larger adhesion forces between the TM cells and ECM, and higher resistance to aqueous humor drainage. TM cells sense the mechanical forces due to IOP dynamic and apply multidimensional forces on the ECM. Recognizing the importance of cellular forces in modulating various cellular activities and development, this study is aimed to develop a 2D in vitro cell culture model to calculate the 3D, depth-dependent, dynamic traction forces, tensile/compressive/shear strain of the normal and glaucomatous human TM cells within a deformable polyacrylamide (PAM) gel substrate. Normal and glaucomatous human TM cells were isolated, cultured, and seeded on top of the PAM gel substrate with embedded FluoSpheres, spanning elastic moduli of 1.5 to 80 kPa. Sixteen-hour post-seeding live confocal microscopy in an incubator was conducted to Z-stack image the 3D displacement map of the FluoSpheres within the PAM gels. Combined with the known PAM gel stiffness, we ascertained the 3D traction forces in the gel. Our results revealed meaningfully larger traction forces in the glaucomatous TM cells compared to the normal TM cells, reaching depths greater than 10-µm in the PAM gel substrate. Stress fibers in TM cells increased with gel rigidity, but diminished when stiffness rose from 20 to 80 kPa. The developed 2D cell culture model aids in understanding how altered mechanical properties in glaucoma impact TM cell behavior and aqueous humor outflow resistance. STATEMENT OF SIGNIFICANCE: Glaucoma, a leading cause of irreversible blindness, is intricately linked to elevated intraocular pressures and their subsequent cellular effects. The trabecular meshwork plays a pivotal role in this mechanism, particularly its interaction with the extracellular matrix. This research unveils an advanced 2D in vitro cell culture model that intricately maps the complex 3D forces exerted by trabecular meshwork cells on the extracellular matrix, offering unparalleled insights into the cellular biomechanics at play in both healthy and glaucomatous eyes. By discerning the changes in these forces across varying substrate stiffness levels, we bridge the gap in understanding between cellular mechanobiology and the onset of glaucoma. The findings stand as a beacon for potential therapeutic avenues, emphasizing the gravity of cellular/extracellular matrix interactions in glaucoma's pathogenesis and setting the stage for targeted interventions in its early stages.

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.

Implementing new computational methods for the study of JCT and SC inner wall basement membrane biomechanics and hydrodynamics

PURPOSE

The conventional aqueous outflow pathway, which includes the trabecular meshwork (TM), juxtacanalicular tissue (JCT), and the inner wall endothelium of Schlemm's canal (SC), regulates intraocular pressure (IOP) by controlling the aqueous humor outflow resistance. Despite its importance, our understanding of the biomechanics and hydrodynamics within this region remains limited. Fluid-structure interaction (FSI) offers a way to estimate the biomechanical properties of the JCT and SC under various loading and boundary conditions, providing valuable insights that are beyond the reach of current imaging techniques.

METHODS

In this study, a normal human eye was fixed at a pressure of 7 mm Hg, and two radial wedges of the TM tissues, which included the SC inner wall basement membrane and JCT, were dissected, processed, and imaged using 3D serial block-face scanning electron microscopy (SBF-SEM). Four different sets of images were used to create 3D finite element (FE) models of the JCT and inner wall endothelial cells of SC with their basement membrane. The outer JCT portion was carefully removed as the outflow resistance is not in that region, leaving only the SCE inner wall and a few µm of the tissue, which does contain the resistance. An inverse iterative FE algorithm was then utilized to calculate the unloaded geometry of the JCT/SC complex at an aqueous humor pressure of 0 mm Hg. Then in the model, the intertrabecular spaces, pores, and giant vacuole contents were replaced by aqueous humor, and FSI was employed to pressurize the JCT/SC complex from 0 to 15 mm Hg.

RESULTS

In the JCT/SC complex, the shear stress of the aqueous humor is not evenly distributed. Areas proximal to the inner wall of SC experience larger stresses, reaching up to 10 Pa, while those closer to the JCT undergo lower stresses, approximately 4 Pa. Within this complex, giant vacuoles with or without I-pore behave differently. Those without I-pores experience a more significant strain, around 14%, compared to those with I-pores, where the strain is roughly 9%.

CONCLUSIONS

The distribution of aqueous humor wall shear stress is not uniform within the JCT/SC complex, which may contribute to our understanding of the underlying selective mechanisms in the pathway.

Biomechanics of the JCT and SC Inner Wall Endothelial Cells with Their Basement Membrane Using 3D Serial Block-Face Scanning Electron Microscopy

BACKGROUND

More than ~70% of the aqueous humor exits the eye through the conventional aqueous outflow pathway that is comprised of the trabecular meshwork (TM), juxtacanalicular tissue (JCT), the inner wall endothelium of Schlemm's canal (SC). The flow resistance in the JCT and SC inner wall basement membrane is thought to play an important role in the regulation of the intraocular pressure (IOP) in the eye, but current imaging techniques do not provide enough information about the mechanics of these tissues or the aqueous humor in this area.

METHODS

A normal human eye was perfusion-fixed and a radial wedge of the TM tissue from a high-flow region was dissected. The tissues were then sliced and imaged using serial block-face scanning electron microscopy. Slices from these images were selected and segmented to create a 3D finite element model of the JCT and SC cells with an inner wall basement membrane. The aqueous humor was used to replace the intertrabecular spaces, pores, and giant vacuoles, and fluid-structure interaction was employed to couple the motion of the tissues with the aqueous humor.

RESULTS

Higher tensile stresses (0.8-kPa) and strains (25%) were observed in the basement membrane beneath giant vacuoles with open pores. The volumetric average wall shear stress was higher in SC than in JCT/SC. As the aqueous humor approached the inner wall basement membrane of SC, the velocity of the flow decreased, resulting in the formation of small eddies immediately after the flow left the inner wall.

CONCLUSIONS

Improved modeling of SC and JCT can enhance our understanding of outflow resistance and funneling. Serial block-face scanning electron microscopy with fluid-structure interaction can achieve this, and the observed micro-segmental flow patterns in ex vivo perfused human eyes suggest a hypothetical mechanism.

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.

The role of conventional and unconventional adaptive routes in lowering of intraocular pressure: Theoretical model and simulation

In this article, we propose a theoretical model leveraging the analogy between fluid and electric variables to investigate the relation among aqueous humor (AH) circulation and drainage and intraocular pressure (IOP), the principal established risk factor of severe neuropathologies of the optic nerve such as glaucoma. IOP is the steady-state result of the balance among AH secretion (AHs), circulation (AHc), and drainage (AHd). AHs are modeled as a given volumetric flow rate electrically corresponding to an input current source. AHc is modeled by the series of two linear hydraulic conductances (HCs) representing the posterior and anterior chambers. AHd is modeled by the parallel of three HCs: a linear HC for the conventional adaptive route (ConvAR), a nonlinear HC for the hydraulic component of the unconventional adaptive route (UncAR), and a nonlinear HC for the drug-dependent component of the UncAR. The proposed model is implemented in a computational virtual laboratory to study the value attained by the IOP under physiological and pathological conditions. Simulation results (i) confirm the conjecture that the UncAR acts as a relief valve under pathological conditions, (ii) indicate that the drug-dependent AR is the major opponent to IOP increase in the case of elevated trabecular meshwork resistance, and (iii) support the use of the model as a quantitative tool to complement in vivo studies and help design and optimize medications for ocular diseases.

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.

Modeling the Endothelial Glycocalyx Layer in the Human Conventional Aqueous Outflow Pathway

A layer of proteoglycans and glycoproteins known as glycocalyx covers the surface of the trabecular meshwork (TM), juxtacanalicular tissue (JCT), and Schlemm's canal (SC) inner wall of the conventional aqueous outflow pathway in the eye. This has been shown to play a role in the mechanotransduction of fluid shear stress and in the regulation of the outflow resistance. The outflow resistance in the conventional outflow pathway is the main determinant of the intraocular pressure (IOP) through an active, two-way, fluid-structure interaction coupling between the outflow tissues and aqueous humor. A 3D microstructural finite element (FE) model of a healthy human eye TM/JCT/SC complex with interspersed aqueous humor was constructed. A very thin charged double layer that represents the endothelial glycocalyx layer covered the surface of the elastic outflow tissues. The aqueous humor was modeled as electroosmotic flow that is charged when it is in contact with the outflow tissues. The electrical-fluid-structure interaction (EFSI) method was used to couple the charged double layer (glycocalyx), fluid (aqueous humor), and solid (outflow tissues). When the IOP was elevated to 15 mmHg, the maximum aqueous humor velocity in the EFSI model was decreased by 2.35 mm/s (9%) compared to the fluid-structure interaction (FSI) model. The charge or electricity in the living human conventional outflow pathway generated by the charged endothelial glycocalyx layer plays a minor biomechanical role in the resultant stresses and strains as well as the hydrodynamics of the aqueous humor.