Factors influencing optic nerve head biomechanics

PURPOSE

The biomechanical environment within the optic nerve head (ONH) may play a role in retinal ganglion cell loss in glaucomatous optic neuropathy. This was a systematic analysis in which finite element methods were used to determine which anatomic and biomechanical factors most influenced the biomechanical response of the ONH to acute changes in IOP.

METHODS

Based on a previously described computational model of the eye, each of 21 input factors, representing the biomechanical properties of relevant ocular tissues, the IOP, and 14 geometric factors were independently varied. The biomechanical response of the ONH tissues was quantified through a set of 29 outcome measures, including peak and mean stress and strain within each tissue, and measures of geometric changes in ONH tissues. Input factors were ranked according to their aggregated influence on groups of outcome measures.

RESULTS

The five input factors that had the largest influence across all outcome measures were, in ranked order: stiffness of the sclera, radius of the eye, stiffness of the lamina cribrosa, IOP, and thickness of the scleral shell. The five least influential factors were, in reverse ranked order: retinal thickness, peripapillary rim height, cup depth, cup-to-disc ratio, and pial thickness. Factor ranks were similar for various outcome measure groups and factor ranges.

CONCLUSIONS

The model predicts that ONH biomechanics are strongly dependent on scleral biomechanical properties. Acute deformations of ONH tissues, and the consequent high levels of neural tissue strain, were less strongly dependent on the action of IOP directly on the internal surface of the ONH than on the indirect effects of IOP on the sclera. This suggests that interindividual variations in scleral properties could be a risk factor for the development of glaucoma. Eye size and lamina cribrosa biomechanical properties also have a strong influence on ONH biomechanics.

Consensus recommendations for trabecular meshwork cell isolation, characterization and culture

Cultured trabecular meshwork (TM) cells are a valuable model system to study the cellular mechanisms involved in the regulation of conventional outflow resistance and thus intraocular pressure; and their dysfunction resulting in ocular hypertension. In this review, we describe the standard procedures used for the isolation of TM cells from several animal species including humans, and the methods used to validate their identity. Having a set of standard practices for TM cells will increase the scientific rigor when used as a model, and enable other researchers to replicate and build upon previous findings.

Biomechanics of the optic nerve head

Biomechanical factors acting at the level of the lamina cribrosa (LC) are postulated to play a role in retinal ganglion cell dysfunction and loss in glaucoma. In support of this postulate, we now know that a number of cell types (including lamina cribrosa cells) are mechanosensitive. Here we briefly review data indicating cellular stretching, rate of stretching and substrate stiffness may be important mechanosensitivity factors in glaucoma. We then describe how experiments and finite element modeling can be used to quantify the biomechanical environment within the LC, and how this environment depends on IOP. Generic and individual-specific models both suggest that peripapillary scleral properties have a strong influence on LC biomechanics, which can be explained by the observation that scleral deformation drives much of the IOP-dependent straining of the LC. Elegant reconstructions of the LC in monkey eyes have shown that local strains experienced by LC cells depend strongly on laminar beam microarchitecture, which can lead to large local strain elevations. Further modeling, suitably informed by experiments, is needed to better understand lamina cribrosa biomechanics and how they may be involved in glaucomatous optic neuropathy.

Finite element modeling of optic nerve head biomechanics

PURPOSE

Biomechanical factors have been implicated in the development of glaucomatous optic neuropathy, particularly at the level of the lamina cribrosa. The goal of this study was to characterize the biomechanics of the optic nerve head using computer modeling techniques.

METHODS

Several models of the optic nerve head tissues (pre- and postlaminar neural tissue, lamina cribrosa, central retinal vessel, sclera, and pia mater) were constructed. Stresses, deformations, and strains were computed using finite element modeling for a range of normal and elevated intraocular pressures. Computed retinal surface deformations were compared with measured deformation patterns in enucleated human eyes. A sensitivity analysis was performed in which tissue properties and selected geometric features were varied.

RESULTS

Acute IOP-induced deformation of the vitreoretinal interface was highly dependent on optic cup shape but showed a characteristic "W-shaped" profile that did not match the deformation of the anterior surface of the lamina cribrosa. The central retinal vasculature had surprisingly little effect on optic nerve head biomechanics. At an IOP of 50 mm Hg, strains (fractional elongation) in the lamina cribrosa averaged 4% to 5.5%, dependent on model geometry, with maximum strains up to 7.7%. Strains in the lamina cribrosa were more dependent on scleral stiffness, scleral thickness, and scleral canal diameter than on lamina cribrosa stiffness and optic cup shape.

CONCLUSIONS

Computed levels of strain in the lamina cribrosa are biologically significant and capable of contributing to the development of glaucomatous optic neuropathy, even without considering the probable accentuating effect of the lamina cribrosa's microarchitecture. Depending on optic cup shape, IOP-induced deformation of the vitreoretinal interface may not match lamina cribrosa deformation. This finding implies that scanning laser tomography has limited ability to estimate lamina cribrosa deformation when imaging the anterior topography of the optic nerve head. Biomechanical effects in the lamina cribrosa depend strongly on scleral properties.

Effects of TGF-beta2 in perfused human eyes

PURPOSE

TGF-beta2 is known to be present at elevated levels in the aqueous humor of patients with primary open-angle glaucoma (POAG). Studies have shown that TGF-beta2 influences cultured trabecular meshwork (TM) cells, but the effects of this cytokine on intact TM and outflow facility have not been studied. The purpose of this study was to investigate whether TGF-beta2 treatment induces changes in outflow facility and morphologic changes in the TM tissue and whether these changes are comparable to those previously recorded in glaucomatous eyes.

METHODS

Baseline facility was measured in paired human eyes (n = 8 pairs), with a constant-flow anterior segment culture system. Medium perfusing experimental eyes was then supplemented with activated human recombinant TGF-beta2 (3.0 ng/mL, comparable to or slightly greater than measured aqueous humor levels in patients with POAG), and facility was measured for at least 8 days. At the conclusion of the perfusion, eyes were fixed and processed for light microscopy, transmission electron microscopy, and immunolabeling studies.

RESULTS

TGF-beta2 perfusion reduced outflow facility by 27% (P = 0.03) and promoted focal accumulation of fine fibrillar extracellular material in multilayered structures under the inner wall of Schlemm's canal. In treated eyes, Schlemm's canal was 27% shorter (P = 0.02), and the length of the inner wall apparently available for fluid flow was 33% less (P = 0.001), both compared with paired control eyes.

CONCLUSIONS

TGF-beta2 reduces outflow facility when perfused into cultured human anterior segments. Furthermore, TGF-beta2 affects the extracellular matrix of the trabecular meshwork in a manner that is consistent with the observed reduction in outflow facility. Although the distribution of accumulated fibrillar material was different in these perfused eyes than that in POAG, the difference could be due to variation in biomechanical environment for TM cells in cultured anterior segments compared with the living eye. Overall, these results support the hypothesis that elevated TGF-beta2 levels in the aqueous humor play a role in the pathogenesis of the ocular hypertension in POAG.

Deformation of the lamina cribrosa by elevated intraocular pressure

The purpose of this study was to determine the mechanical response of the lamina cribrosa (LC) to elevated intraocular pressure (IOP) so as to identify possible mechanisms of optic nerve damage in early glaucoma. Ten pairs of normal human eyes were fixed after 24 hours' exposure to 50 mm Hg pressure (experimental eyes) or 5 mm Hg pressure (contralateral control eyes). Photomicrographs of the central region of the optic nerve head (ONH) were taken to examine the LC morphologically and to measure the dimensions of the LC. It was found that elevated IOP caused the LC to deflect posteriorly without affecting its thickness. The majority of the posterior displacement in the LC occurred near the periphery of the ONH. This shape change is consistent with a model of force distribution within the LC in which shear stresses are dominant; such stresses are maximal at the periphery and minimal at the centre of the ONH. These findings support a model in which mechanical forces, specifically shearing stresses within the peripheral lamina, play a direct role in the pathology of glaucomatous optic neuropathy.

Ocular biomechanics and biotransport

The eye transduces light, and we usually do not think of it as a biomechanical structure. Yet it is actually a pressurized, thick-walled shell that has an internal and external musculature, a remarkably complex internal vascular system, dedicated fluid production and drainage tissues, and a variety of specialized fluid and solute transport systems. Biomechanics is particularly involved in accommodation (focusing near and far), as well as in common disorders such as glaucoma, macular degeneration, myopia, and presbyopia. In this review, we give a (necessarily brief) overview of many of the interesting biomechanical aspects of the eye, concluding with a list of open problems.

Factors influencing blood flow patterns in the human right coronary artery

Evidence suggests that atherogenesis is linked to local hemodynamic factors such as wall shear stress. We investigated the velocity and wall shear stress patterns within a human right coronary artery (RCA), an important site of atherosclerotic lesion development. Emphasis was placed on evaluating the effect of flow waveform and inlet flow velocity profile on the hemodynamics in the proximal, medial, and distal arterial regions. Using the finite-element method, velocity and wall shear stress patterns in a rigid, anatomically realistic model of a human RCA were computed. Steady flow simulations (ReD=500) were performed with three different inlet velocity profiles; pulsatile flow simulations utilized two different flow waveforms (both with Womersley parameter=1.82, mean ReD=233), as well as two of the three inlet profiles. Velocity profiles showed Dean-like secondary flow features that were remarkably sensitive to the local curvature of the RCA model. Particularly noteworthy was the "rotation" of these Dean-like profiles, which produced large local variations in wall shear stress along the sidewalls of the RCA model. Changes in the inlet velocity profiles did not produce significant changes in the arterial velocity and wall shear stress patterns. Pulsatile flow simulations exhibited remarkably similar cycle-average wall shear stress distributions regardless of waveform and inlet velocity profile. The oscillatory shear index was very small and was attributed to flow reversal in the waveform, rather than separation. Cumulatively, these results illustrate that geometric effects (particularly local three-dimensional curvature) dominate RCA hemodynamics, implying that studies attempting to link hemodynamics with atherogenesis should replicate the patient-specific RCA geometry.

Simulation of non-Newtonian blood flow in an end-to-side anastomosis

In order to quantify the importance of non-Newtonian blood rheology on anastomotic flow patterns, the characteristics of Newtonian and non-Newtonian blood flows were compared in a 2-D, 45 degree end-to-side anastomosis model under both steady and unsteady flow conditions. All flows were assumed to be two-dimensional, and were simulated numerically using parameters consistent with blood flow in the femoral artery. A novel, purely viscous constitutive relation, based on a generalized form of the power law relation, was developed to model the non-Newtonian rheology of blood. The resulting wall shear stress patterns indicate that for the parameter values under consideration, non-Newtonian blood rheology has a significant effect on steady flow wall shear stresses, but no significant effect on unsteady flow wall shear stresses. Based on these and other simulations, a parameter is formulated that gives an indication of the importance of non-Newtonian blood rheology under a given set of flow conditions. In addition, an argument is presented for allowing the conclusions from this two-dimensional study to be extended to three-dimensional blood flow.

Dexamethasone alters F-actin architecture and promotes cross-linked actin network formation in human trabecular meshwork tissue

Elevated intraocular pressure is an important risk factor for the development of glaucoma, a leading cause of irreversible blindness. This ocular hypertension is due to increased hydrodynamic resistance to the drainage of aqueous humor through specialized outflow tissues, including the trabecular meshwork (TM) and the endothelial lining of Schlemm's canal. We know that glucocorticoid therapy can cause increased outflow resistance and glaucoma in susceptible individuals, that the cytoskeleton helps regulate aqueous outflow resistance, and that glucocorticoid treatment alters the actin cytoskeleton of cultured TM cells. Our purpose was to characterize the actin cytoskeleton of cells in outflow pathway tissues in situ, to characterize changes in the cytoskeleton due to dexamethasone treatment in situ, and to compare these with changes observed in cell culture. Human ocular anterior segments were perfused with or without 10(-7) M dexamethasone, and F-actin architecture was investigated by confocal laser scanning microscopy. We found that outflow pathway cells contained stress fibers, peripheral actin staining, and occasional actin "tangles." Dexamethasone treatment caused elevated IOP in several eyes and increased overall actin staining, with more actin tangles and the formation of cross-linked actin networks (CLANs). The actin architecture in TM tissues was remarkably similar to that seen in cultured TM cells. Although CLANs have been reported previously in cultured cells, this is the first report of CLANs in tissue. These cytoskeletal changes may be associated with increased aqueous humor outflow resistance after ocular glucocorticoid treatment.

eNOS, a pressure-dependent regulator of intraocular pressure

PURPOSE

Pathology in the primary drainage pathway for aqueous humor in the eye is responsible for ocular hypertension, the only treatable risk factor in patients with glaucoma. Unfortunately, the mechanisms that regulate pressure-dependent drainage of aqueous humor and thus intraocular pressure (IOP) are unknown. To better understand one possible underlying molecular factor that regulates IOP, nitric oxide (NO), pressure-dependent drainage in transgenic mice overexpressing endothelial NO synthase (eNOS) was studied.

METHODS

IOP was measured by rebound tonometry in mice, and pressure versus flow data were measured by ex vivo perfusion at multiple pressures between 8 and 45 mm Hg, using mock AH ±100 μM L-NAME. A subset of eyes was examined histologically using standard techniques or was assayed for fusion protein expression by Western blot analysis.

RESULTS

IOP was lower (9.6 ± 2.7 vs. 11.4 ± 2.5 mm Hg; mean ± SD; P = 0.04) and pressure-dependent drainage was higher (0.0154 ± 0.006 vs. 0.0066 ± 0.0009 μL/min/mm Hg; P = 0.002) in the transgenic mice than in the wild-type animals; however, pressure-independent drainage was unaffected. The NOS inhibitor L-NAME normalized pressure-dependent drainage in transgenic animals. For IOP >35 mm Hg, the slope of the pressure-flow curve in wild-type mice increased to match that seen in transgenic mice. Shear stress in the pressure-dependent pathway at elevated pressures was calculated to be in a range known to affect eNOS expression and activity in vascular endothelia.

CONCLUSIONS

Endothelial NOS overexpression lowers IOP by increasing pressure-dependent drainage in the mouse eye. Data are consistent with NO's having a mechanoregulatory role in aqueous humor dynamics, with eNOS induction at elevated IOPs leading to increased pressure-dependent outflow.

Effects of cardiac motion on right coronary artery hemodynamics

The purpose of this work was to investigate the effects of physiologically realistic cardiac-induced motion on hemodynamics in human right coronary arteries. The blood flow patterns were numerically simulated in a modeled right coronary artery (RCA) having a uniform circular cross section of 2.48 mm diam. Arterial motion was specified based on biplane cineangiograms, and incorporated physiologically realistic bending and torsion. Simulations were carried out with steady and pulsatile inflow conditions (mean ReD=233, alpha=1.82) in both fixed and moving RCA models, to evaluate the relative importance of RCA motion, flow pulsation, and the interaction between motion and flow pulsation. RCA motion with a steady inlet flow rate caused variations in wall shear stress (WSS) magnitude up to 150% of the inlet Poiseuille value. There was significant spatial variability in the magnitude of this motion-induced WSS variation. However, the time-averaged WSS distribution was similar to that predicted in a static model representing the time-averaged geometry. Furthermore, the effects of flow pulsatility dominated RCA motion-induced effects; specifically, there were only modest differences in the WSS history between simulations conducted in fixed and moving RCA models with pulsatile inflow. RCA motion has little effect on time-averaged WSS patterns. It has a larger effect on the temporal variation of WSS, but even this effect is overshadowed by the variations in WSS due to flow pulsation. The hemodynamic effects of RCA motion can, therefore, be ignored as a first approximation in modeling studies.

Finite element modeling of the human sclera: influence on optic nerve head biomechanics and connections with glaucoma

Scleral thickness, especially near the optic nerve head (ONH), is a potential factor of interest in the development of glaucomatous optic neuropathy. Large differences in the dimensions of the sclera, the principal load-bearing tissue of the eye, have been observed between individuals. This study aimed to characterize the effects of these differences on ONH biomechanics. Eleven enucleated human globes (7 normal and 4 ostensibly glaucomatous) were imaged using high-field microMRI and segmented to produce 3-D individual-specific corneoscleral shells. An identical, idealized ONH geometry was inserted into each shell. Finite element modeling predicted the effects of pressurizing the eyes to an IOP of 30 mmHg, with the results used to characterize the effect of inter-individual differences in scleral dimensions on the biomechanics of the ONH. Measurements of the individual-specific corneoscleral shells were used to construct a 2-D axisymmetric idealized model of the corneoscleral shell and ONH. A sensitivity analysis based on this model quantified the relative importance of different geometrical characteristics of the scleral shell on the biomechanics of the ONH. Significant variations were observed in various measures of strain in the idealized lamina cribrosa (LC) across the seven normal corneoscleral shells, implying large differences in individual biomechanics due to scleral anatomy variations alone. The sensitivity analysis revealed that scleral thickness adjacent to the ONH was responsible for the vast majority of variation. Remarkably, varying peripapilary scleral thickness over the physiologically measured range changed the peak (95th percentile) first principal strain in the LC and radial displacement of the ONH canal by an amount that was equivalent to a change in IOP of 15 mmHg. Inter-individual variations in scleral thickness, particularly peripapillary scleral thickness, can result in vastly different biomechanical responses to IOP. These differences may be significant for understanding the interactions between IOP and scleral biomechanics in the pathogenesis of glaucomatous optic neuropathy. The relationship between scleral thickness and material properties needs to be studied in human eyes.

Predicted extension, compression and shearing of optic nerve head tissues

Glaucomatous optic neuropathy may be in part due to an altered biomechanical environment within the optic nerve head (ONH) produced by an elevated intraocular pressure (IOP). Previous work has characterized the magnitude of the IOP-induced deformation of ONH tissues but has not focused specifically on the mode of deformation (strain), i.e. whether the ONH tissues and cells are stretched, compressed or sheared. Circumstantial evidence indicates that the mode of deformation has biological consequences. Here we use computational models to study the different modes of deformation that occur in an ONH as a result of an increase in IOP. One generic and three individual-specific 3D models of the human ONH were reconstructed as previously described. Each model consisted of five tissue regions: pre and post-laminar neural tissue, lamina cribrosa, sclera and pia mater. Finite element methods were then used to predict the biomechanical response to changes in IOP. For each model we computed six local measures of strain, including the magnitude and direction of maximum stretching, maximum compression and maximum shearing strain. We compared the spatial and population distributions of the various measures of strain by using semi-quantitative (contour plots) and quantitative (histograms) methods. For all models, as IOP increased, the tissues of the ONH were subjected simultaneously to various modes of strain, including compression, extension and shearing. The highest magnitudes of all modes of strain occurred within the neural tissue regions. There were substantial differences in the magnitudes of the various modes of strain, with the largest strains being in compression, followed by shearing and finally by extension. The biomechanical response of an individual-specific ONH to changes in IOP is complex and cannot be fully captured by one measure of deformation. We predict that cells within the ONH are subjected to very different modes of deformation as IOP increases. The largest deformations are compressive, followed by shearing and stretching. Models of IOP-induced RGC damage need to be further refined by characterizing the cellular response to these different modes of strain.

Schlemm's canal and primary open angle glaucoma: correlation between Schlemm's canal dimensions and outflow facility

The purpose of the study was to examine the correlation between outflow facility and morphometric measurements of Schlemm's canal (SC) in normal and glaucomatous human eyes. Outflow facility was measured in ten normal and five glaucomatous (POAG) human eyes prior to perfusion-fixation at a pressure of 15 mmHg. Two separate sections per quadrant (eight sections per eye) were analysed using a morphometric analysis system. SC cross sectional area, perimeter, and inner wall length were measured. Inner wall endothelial nuclei were counted. When pooled data were examined, a statistically significant correlation was found between facility and SC area (P < or = 0.01). Most importantly, POAG eyes had a significantly smaller SC cross-sectional area, SC perimeter and SC inner wall length compared to normal eyes. If examined as a reduction in SC filtering area, this decrease in SC inner wall length alone could account for approximately 41-55% of the difference in outflow facility observed between normal and POAG eyes. The dimensions of Schlemm's canal in glaucomatous human eyes were significantly smaller than those in normal eyes. This reduction in SC dimensions may account for approximately half of the decrease in outflow facility observed in POAG eyes.

Modeling individual-specific human optic nerve head biomechanics. Part I: IOP-induced deformations and influence of geometry

Glaucoma, the second most common cause of blindness worldwide, is an ocular disease characterized by progressive loss of retinal ganglion cell (RGC) axons. Biomechanical factors are thought to play a central role in RGC loss, but the specific mechanism underlying this disease remains unknown. Our goal was to characterize the biomechanical environment in the optic nerve head (ONH)--the region where RGC damage occurs--in human eyes. Post mortem human eyes were imaged, fixed at either 5 or 50 mmHg pressure and processed histologically to acquire serial sections through the ONH. Three-dimensional models of the ONH region were reconstructed from these sections and embedded in a generic scleral shell to create a model of an entire eye. We used finite element simulations to quantify the effects of an acute change in intraocular pressure from 5 to 50 mmHg on the ONH biomechanical environment. Computed strains varied substantially within the ONH, with the pre-laminar neural tissue and the lamina cribrosa showing the greatest strains. The mode of strain having the largest magnitude was third principal strain (compression), reaching 12-15% in both the lamina cribrosa and the pre-laminar neural tissue. Shear strains were also substantial. The distribution of strains in all ONH tissues was remarkably similar between eyes. Inter-individual variations in ONH geometry (anatomy) have only modest effects on ONH biomechanics, and may not explain inter-individual susceptibility to elevated intraocular pressure. Consistent with previous results using generic ONH models, the displacements of the vitreo-retinal interface and the anterior surface of the lamina cribrosa can differ substantially, suggesting that currently available optical imaging methods do not provide information of the acute deformations within ONH tissues. Predicted strains within ONH tissues are potentially biologically significant and support the hypothesis that biomechanical factors contribute to the initial insult that leads to RGC loss in glaucoma.

Biomechanics of Schlemm's canal endothelium and intraocular pressure reduction

Ocular hypertension in glaucoma develops due to age-related cellular dysfunction in the conventional outflow tract, resulting in increased resistance to aqueous humor outflow. Two cell types, trabecular meshwork (TM) and Schlemm's canal (SC) endothelia, interact in the juxtacanalicular tissue (JCT) region of the conventional outflow tract to regulate outflow resistance. Unlike endothelial cells lining the systemic vasculature, endothelial cells lining the inner wall of SC support a transcellular pressure gradient in the basal to apical direction, thus acting to push the cells off their basal lamina. The resulting biomechanical strain in SC cells is quite large and is likely to be an important determinant of endothelial barrier function, outflow resistance and intraocular pressure. This review summarizes recent work demonstrating how biomechanical properties of SC cells impact glaucoma. SC cells are highly contractile, and such contraction greatly increases cell stiffness. Elevated cell stiffness in glaucoma may reduce the strain experienced by SC cells, decrease the propensity of SC cells to form pores, and thus impair the egress of aqueous humor from the eye. Furthermore, SC cells are sensitive to the stiffness of their local mechanical microenvironment, altering their own cell stiffness and modulating gene expression in response. Significantly, glaucomatous SC cells appear to be hyper-responsive to substrate stiffness. Thus, evidence suggests that targeting the material properties of SC cells will have therapeutic benefits for lowering intraocular pressure in glaucoma.

Modeling individual-specific human optic nerve head biomechanics. Part II: influence of material properties

Biomechanical factors acting within the optic nerve head (ONH) likely play a role in the loss of vision that occurs in glaucoma. In a companion paper (Sigal et al. 2008), we quantified the biomechanical environment within individual-specific ONH models reconstructed from human post mortem eyes. Our goal in this manuscript was to use finite element modeling to investigate the influence of tissue material properties on ONH biomechanics in these same individual-specific models. A sensitivity analysis was carried out by simulating the effects of changing intraocular pressure on ONH biomechanics as tissue mechanical properties were systematically varied over ranges reported in the literature. This procedure was repeated for each individual-specific model described in the companion paper (Sigal et al. 2008). The outcome measures of the analysis were first and third principal strains, as well as the derived quantity of maximum shear strain, in ONH tissues. Scleral stiffness had by far the largest influence in ONH biomechanics, and this result was remarkably consistent across ONH models. The stiffnesses of the lamina cribrosa and pia mater were also influential. These results are consistent with those obtained using generic ONH models. The compressibility of the pre-laminar neural tissue influenced compressive and shearing strains. Overall, tissue material properties had a much greater influence on ONH biomechanics than did tissue geometry, as assessed by comparing results between our individual-specific models. Material properties of ONH tissues, particularly of the peripapillary sclera, play a dominant role in the mechanical response of an ONH to acute changes in IOP and may be important in the pathogenesis of glaucoma. We need to better understand inter-individual differences in scleral biomechanical properties and whether they are clinically important.

Requirements for mesh resolution in 3D computational hemodynamics

Computational techniques are widely used for studying large artery hemodynamics. Current trends favor analyzing flow in more anatomically realistic arteries. A significant obstacle to such analyses is generation of computational meshes that accurately resolve both the complex geometry and the physiologically relevant flow features. Here we examine, for a single arterial geometry, how velocity and wall shear stress patterns depend on mesh characteristics. A well-validated Navier-Stokes solver was used to simulate flow in an anatomically realistic human right coronary artery (RCA) using unstructured high-order tetrahedral finite element meshes. Velocities, wall shear stresses (WSS), and wall shear stress gradients were computed on a conventional "high-resolution" mesh series (60,000 to 160,000 velocity nodes) generated with a commercial meshing package. Similar calculations were then performed in a series of meshes generated through an adaptive mesh refinement (AMR) methodology. Mesh-independent velocity fields were not very difficult to obtain for both the conventional and adaptive mesh series. However, wall shear stress fields, and, in particular, wall shear stress gradient fields, were much more difficult to accurately resolve. The conventional (nonadaptive) mesh series did not show a consistent trend towards mesh-independence of WSS results. For the adaptive series, it required approximately 190,000 velocity nodes to reach an r.m.s. error in normalized WSS of less than 10 percent. Achieving mesh-independence in computed WSS fields requires a surprisingly large number of nodes, and is best approached through a systematic solution-adaptive mesh refinement technique. Calculations of WSS, and particularly WSS gradients, show appreciable errors even on meshes that appear to produce mesh-independent velocity fields.

Ultrastructural changes associated with dexamethasone-induced ocular hypertension in mice

PURPOSE

To determine whether dexamethasone (DEX)-induced ocular hypertension (OHT) in mice mimics the hallmarks of steroid-induced glaucoma (SIG) in humans, including reduced conventional outflow facility (C), increased extracellular matrix (ECM), and myofibroblasts within the outflow pathway.

METHODS

Osmotic mini-pumps were implanted subcutaneously into C57BL/6J mice for systemic delivery of DEX (3-4 mg/kg/d, n = 31 mice) or vehicle (n = 28). IOP was measured weekly by rebound tonometry. After 3 to 4 weeks, mice were euthanized and eyes enucleated for ex vivo perfusion to measure C, for electron microscopy to examine the trabecular meshwork (TM) and Schlemm's canal (SC), or for immunohistochemistry to examine type IV collagen and α-smooth muscle actin. The length of basement membrane material (BMM) was measured along the anterior-posterior extent of SC by electron microscopy. Ultrastructural changes in BMM of DEX-treated mice were compared against archived human SIG specimens.

RESULTS

Dexamethasone increased IOP by 2.6 ± 1.6 mm Hg (mean ± SD) over 3 to 4 weeks and decreased C by 52% ± 17% versus controls. Intraocular pressure elevation correlated with decreased C. Dexamethasone treatment led to increased fibrillar material in the TM, plaque-like sheath material surrounding elastic fibers, and myofibroblasts along SC outer wall. The length of BMM underlying SC was significantly increased in mice with DEX and in humans with SIG, and in mice decreased C correlated with increased BMM.

CONCLUSIONS

Dexamethasone-induced OHT in mice mimics hallmarks of human SIG within 4 weeks of DEX treatment. The correlation between reduced C and newly formed ECM motivates further study using DEX-treated mice to investigate the pathogenesis of conventional outflow obstruction in glaucoma.

Computational modeling of mass transfer and links to atherosclerosis

In the context of atherogenesis, mass transport refers to the movement of atherogenic molecules from flowing blood into the artery wall, or vice versa. Although LDL transport clearly plays a role in atherosclerotic plaque development, it is much less clear whether abnormalities in mass transfer patterns are in themselves atherogenic. A powerful way of addressing this question is through computational modeling, which provides detailed descriptions of local mass transport features. Here we briefly review the strategy and some of the pros and cons of such a modeling approach, and then focus on results gained from studies in a variety of arterial geometries. The general picture is that zones of hypoxia (low oxygen transport from blood to wall) and elevated LDL tend to colocalize with each other, and with areas of atherosclerotic lesion development and/or intimal thickening. The picture is complicated by the fact that such zones also tend to have "abnormal" wall shear stress patterns, which are also believed to be atherogenic. Taken together, these results suggest, but do not prove, a role for mass transport in atherogenesis.

Accuracy of computational hemodynamics in complex arterial geometries reconstructed from magnetic resonance imaging

PURPOSE

Combining computational blood flow modeling with three-dimensional medical imaging provides a new approach for studying links between hemodynamic factors and arterial disease. Although this provides patient-specific hemodynamic information, it is subject to several potential errors. This study quantifies some of these errors and identifies optimal reconstruction methodologies.

METHODS

A carotid artery bifurcation phantom of known geometry was imaged using a commercial magnetic resonance (MR) imager. Three-dimensional models were reconstructed from the images using several reconstruction techniques, and steady and unsteady blood flow simulations were performed. The carotid bifurcation from a healthy, human volunteer was then imaged in vivo, and geometric models were reconstructed.

RESULTS

Reconstructed models of the phantom showed good agreement with the gold standard geometry, with a mean error of approximately 15% between the computed wall shear stress fields. Reconstructed models of the in vivo carotid bifurcation were unacceptably noisy, unless lumenal profile smoothing and approximating surface splines were used.

CONCLUSIONS

All reconstruction methods gave acceptable results for the phantom model, but in vivo models appear to require smoothing. If proper attention is paid to smoothing and geometric fidelity issues, models reconstructed from MR images appear to be suitable for use in computational studies of in vivo hemodynamics.

Steady and pulsatile flow fields in an end-to-side arterial anastomosis model

We investigated the flow field within a rigid-walled in vitro model of an end-to-side 45 degree anastomosis in an attempt to identify possible hemodynamic factors that may contribute to the pathogenesis of distal anastomotic intimal hyperplasia. A high-resolution photochromic tracer technique was used to visualize the flow in orthogonal planes and to determine the axial wall shear stress profiles for both steady and pulsatile flows over a range of physiologically relevant conditions. The flow field showed qualitative similarities to those seen in curved vessel: rapidly moving fluid from the graft section affects the bed of the host vessel, that is, the wall opposite the anastomosis, eventually advancing down the host vessel in a spiraling motion. A small mobile separation zone was noted at the toe of the anastomosis. Comparison of wall shear stress profiles with previously reported preferential sites for the development of intimal hyperplasia supported a low wall shear stress and/or flow separation pathogenesis hypothesis. One notable exception was the bed of the host artery that appeared to be subjected to a complex hemodynamic environment.