Full-field deformation of bovine cornea under constrained inflation conditions

Spatial targeted delivery of riboflavin with a controlled corneal iontophoresis delivery system in theranostic-guided UV-A light photo-therapy

Seven human donor eye globes underwent corneal cross-linking using theranostic UV-A device with accessory corneal iontophoresis system for patterned delivery of a 0.22% riboflavin solution. Theranostic-guided UV-A light illumination assessed riboflavin distribution and treated corneas at 10 mW/cm2 for 9 min with a 5.0-mm beam size. Corneal topography maps were taken at baseline and 2-h post-treatment. Analysis utilized corneal topography elevation data, with results showing controlled riboflavin delivery led to a consistent gradient, with 40% higher levels centrally (248 ± 79 μg/cm3) than peripherally (180 ± 72 μg/cm3 at ±2.5 mm from the center). Theranostic-guided UV-A light irradiation resulted in significant changes in corneal topography, with a decrease in best-fit sphere value (-0.7 ± 0.2 D; p < 0.001) and consistent downward shift in corneal elevation map (-11.7 ± 3.7 μm). The coefficient of variation was 2.5%, indicating high procedure performance in achieving significant and reliable corneal flattening.

Structural control of corneal transparency, refractive power and dynamics

The cornea needs to be transparent to visible light and precisely curved to provide the correct refractive power. Both properties are governed by its structure. Corneal transparency arises from constructive interference of visible light due to the relatively ordered arrangement of collagen fibrils in the corneal stroma. The arrangement is controlled by the negatively charged proteoglycans surrounding the fibrils. Small changes in fibril organisation can be tolerated but larger changes cause light scattering. Corneal keratocytes do not scatter light because their refractive index matches that of the surrounding matrix. When activated, however, they become fibroblasts that have a lower refractive index. Modelling shows that this change in refractive index significantly increases light scatter. At the microscopic level, the corneal stroma has a lamellar structure, the parallel collagen fibrils within each lamella making a large angle with those of adjacent lamellae. X-ray scattering has shown that the lamellae have preferred orientations in the human cornea: inferior-superior and nasal-temporal in the central cornea and circumferential at the limbus. The directions at the centre of the cornea may help withstand the pull of the extraocular muscles whereas the pseudo-circular arrangement at the limbus supports the change in curvature between the cornea and sclera. Elastic fibres are also present; in the limbus they contain fibrillin microfibrils surrounding an elastin core, whereas at the centre of the cornea, they exist as thin bundles of fibrillin-rich microfibrils. We present a model based on the structure described above that may explain how the cornea withstands repeated pressure changes due to the ocular pulse.

The ocular mathematical virtual simulator: A validated multiscale model for hemodynamics and biomechanics in the human eye

We present our continuous efforts from a modeling and numerical viewpoint to develop a powerful and flexible mathematical and computational framework called Ocular Mathematical Virtual Simulator (OMVS). The OMVS aims to solve problems arising in biomechanics and hemodynamics within the human eye. We discuss our contribution towards improving the reliability and reproducibility of computational studies by performing a thorough validation of the numerical predictions against experimental data. The OMVS proved capable of simulating complex multiphysics and multiscale scenarios motivated by the study of glaucoma. Furthermore, its modular design allows the continuous integration of new models and methods as the research moves forward, and supports the utilization of the OMVS as a promising non-invasive clinical investigation tool for personalized research in ophthalmology.

Characterization of the Immediate and Delayed Biomechanical Response to UV-A Crosslinking of Human Corneas

PURPOSE

Keratoconus leads to visual deterioration due to irregular astigmatism and corneal thinning. Riboflavin-based corneal UV-A crosslinking (CXL) induces novel intramolecular and intermolecular links resulting in corneal tissue stiffening, thereby halting disease progression. The purpose of this study was to analyze the immediate and delayed biomechanical responses of human donor corneas to CXL.

METHODS

CXL was performed according to the Dresden protocol to corneas not suitable for transplantation. Biomechanical properties were subsequently monitored by measuring the Young modulus using nanoindentation. The immediate tissue response was determined after 0, 1, 15, and 30 minutes of irradiation. Delayed biomechanical effects were investigated with follow-up measurements immediately and 1, 3, and 7 days after CXL.

RESULTS

Young's modulus indicated a linear trend in direct response to increasing irradiation times (mean values: total 61.31 kPa [SD 25.53], 0 minutes 48.82 kPa [SD 19.73], 1 minute 53.44 kPa [SD 25.95], 15 minutes 63.56 kPa [SD 20.99], and 30 minutes 76.76 kPa [SD 24.92]). The linear mixed model for the elastic response of corneal tissue was 49.82 kPa + (0.91 kPa/min × time [minutes]); P < 0.001. The follow-up measurements showed no significant delayed changes in the Young modulus (mean values: total 55,28 kPa [SD 15.95], immediately after CXL 56,83 kPa [SD 18.74], day 1 50.28 kPa [SD 14.15], day 3 57.08 kPa [SD 14.98], and day 7 56.83 kPa [SD 15.07]).

CONCLUSIONS

This study suggests a linear increase of corneal Young modulus as a function of CXL timing. No significant short-term delayed biomechanical changes posttreatment were observed.

Corneal Permeability and Uptake of Twenty-Five Drugs: Species Comparison and Quantitative Structure-Permeability Relationships

The purpose of this study was to determine corneal permeability and uptake in rabbit, porcine, and bovine corneas for twenty-five drugs using an N-in-1 (cassette) approach and relate these parameters to drug physicochemical properties and tissue thickness through quantitative structure permeability relationships (QSPRs). A twenty-five-drug cassette containing β-blockers, NSAIDs, and corticosteroids in solution at a micro-dose was exposed to the epithelial side of rabbit, porcine, or bovine corneas mounted in a diffusion chamber, and the corneal drug permeability and tissue uptake were monitored using an LC-MS/MS method. Data obtained were used to construct and evaluate over 46,000 quantitative structure-permeability (QSPR) models using multiple linear regression, and the best-fit models were cross-validated by Y-randomization. Drug permeability was generally higher in rabbit cornea and comparable between bovine and porcine corneas. Permeability differences between species could be explained in part by differences in corneal thickness. Corneal uptake between species correlated with a slope close to 1, indicating generally similar drug uptake per unit weight of tissue. A high correlation was observed between bovine, porcine, and rabbit corneas for permeability and between bovine and porcine corneas for uptake (R² ≥ 0.94). MLR models indicated that drug characteristics such as lipophilicity (LogD), heteroatom ratio (HR), nitrogen ratio (NR), hydrogen bond acceptors (HBA), rotatable bonds (RB), index of refraction (IR), and tissue thickness (TT) are of great influence on drug permeability and uptake. When data for all species along with thickness as a parameter was used in MLR, the best fit equation for permeability was Log (% transport/cm²·s) = 0.441 LogD - 8.29 IR + 8.357 NR - 0.279 HBA - 3.833 TT + 10.432 (R² = 0.826), and the best-fit equation for uptake was Log (%/g) = 0.387 LogD + 4.442 HR + 0.105 RB - 0.303 HBA - 2.235 TT + 1.422 (R² = 0.750). Thus, it is feasible to explain corneal drug delivery in three species using a single equation.

Typical localised element-specific finite element anterior eye model

Purpose

The study presents an averaged anterior eye geometry model combined with a localised material model that is straightforward, appropriate and amenable for implementation in finite element (FE) modelling.

Methods

Both right and left eye profile data of 118 subjects (63 females and 55 males) aged 22-67 years (38.5 ± 7.6) were used to build an averaged geometry model. Parametric representation of the averaged geometry model was achieved through two polynomials dividing the eye into three smoothly connected volumes. This study utilised the collagen microstructure x-ray data of 6 ex-vivo healthy human eyes, 3 right eyes and 3 left eyes in pairs from 3 donors, 1 male and 2 females aged between 60 and 80 years, to build a localised element-specific material model for the eye.

Results

Fitting the cornea and the posterior sclera sections to a 5th-order Zernike polynomial resulted in 21 coefficients. The averaged anterior eye geometry model recorded a limbus tangent angle of 37° at a radius of 6.6 mm from the corneal apex. In terms of material models, the difference between the stresses generated in the inflation simulation up to 15 mmHg in the ring-segmented material model and localised element-specific material model were significantly different (p < 0.001) with the ring-segmented material model recording average Von-Mises stress 0.0168 ± 0.0046 MPa and the localised element-specific material model recording average Von-Mises stress 0.0144 ± 0.0025 MPa.

Conclusions

The study illustrates an averaged geometry model of the anterior human eye that is easy to generate through two parametric equations. This model is combined with a localised material model that can be used either parametrically through a Zernike fitted polynomial or non-parametrically as a function of the azimuth angle and the elevation angle of the eye globe. Both averaged geometry and localised material models were built in a way that makes them easy to implement in FE analysis without additional computation cost compared to the limbal discontinuity so-called idealised eye geometry model or ring-segmented material model.

Biaxial hyperelastic and anisotropic behaviors of the corneal anterior central stroma along the preferential fibril orientations. Part II: Quantitative computational analysis of mechanical response of stromal components

To study the hyperelastic and anisotropic behaviors of the central anterior stroma for patients with myopia, 40 corneal stromal specimens extracted after small incision lenticule extraction (SMILE) surgery were used in the biaxial extension test along two preferential fibril orientations. An improved collagen fibril crimping constitutive model with a specific physical meaning was proposed to analyze the hyperelasticity and anisotropy of the stroma. The effective elastic modulus of the two families of preferentially oriented collagen fibrils and the stiffness of the non-collagenous matrix along all three directions were compared according to the specific physical meaning of the parameters. Anisotropic behavior was found in the hyperelastic properties of the corneal anterior central stroma in the preferential fibril orientations. The stiffness of non-collagenous matrix is significantly larger in the optical axis direction than in the nasal-temporal (NT) and superior-inferior (SI) directions. Moreover, individual differences between males and females slightly impact on hyperelastic and anisotropic behaviors. The differences of these behaviors were significant in the comparison of the left and right eyes. These results have a guiding significance for the accurate design of surgical plans for refractive surgery according to a patient's condition and have a driving value for the further exploration of the biomechanical properties of the whole cornea.

Anisotropic Mechanical Response of Bovine Pericardium Membrane Through Bulge Test and In-Situ Confocal-Laser Scanning

In this work, we present a new experimental setup for the assessment of the anisotropic properties of Bovine Pericardium (BP) membranes. The chemically fixed BP samples have been subjected to a bulge test with in situ confocal laser scanning at increasing applied pressure. The high resolution topography provided by the confocal laser scanning has allowed to obtain a quantitative measure of the bulge displacement; after polynomial fitting, principal curvatures have been obtained and a degree of anisotropy (DA) has been defined as the normalized difference between the maximum and minimum principal curvatures. The experiments performed on the BP membranes have allowed us to obtain pressure-displacement data which clearly exhibit distinct principal curvatures indicating an anisotropic response. A comparison with curvatures data obtained on isotropic Nitrile Buthadiene Rubber (NBR) samples has confirmed the effectiveness of the experimental setup for this specific purpose. Numerical simulations of the bulge tests have been performed with the purpose of identifying a range of constitutive parameters which well describes the obtained range of DA on the BP membranes. The DA values have been partially validated with biaxial tests available in literature and with suitably performed uni-axial tensile tests.

In Vivo Biomechanical Measurements of the Cornea

In early corneal examinations, the relationships between the morphological and biomechanical features of the cornea were unclear. Although consistent links have been demonstrated between the two in certain cases, these are not valid in many diseased states. An accurate assessment of the corneal biomechanical properties is essential for understanding the condition of the cornea. Studies on corneal biomechanics in vivo suggest that clinical problems such as refractive surgery and ectatic corneal disease are closely related to changes in biomechanical parameters. Current techniques are available to assess the mechanical characteristics of the cornea in vivo. Accordingly, various attempts have been expended to obtain the relevant mechanical parameters from different perspectives, using the air-puff method, ultrasound, optical techniques, and finite element analyses. However, a measurement technique that can comprehensively reflect the full mechanical characteristics of the cornea (gold standard) has not yet been developed. We review herein the in vivo measurement techniques used to assess corneal biomechanics, and discuss their advantages and limitations to provide a comprehensive introduction to the current state of technical development to support more accurate clinical decisions.

Delayed reorganisation of F-actin cytoskeleton and reversible chromatin condensation in scleral fibroblasts under simulated pathological strain

Mechanical loading regulates the functional capabilities of the ocular system, particularly in the sclera (‘white of the eye’) – the principal load-bearing tissue of the ocular globe. Resident fibroblasts of the scleral eye wall are continuously subjected to fluctuating mechanical strains arising from eye movements, cerebrospinal fluid pressure and, most influentially, intra-ocular pressure (IOP). Whilst fibroblasts are hypothesised to actively participate in scleral biomechanics, to date limited information has been reported on how the macroscopic stresses and strains are transmitted via their cytoskeletal networks. In this study, the effect of applying either a ‘physiological load’ (simulating healthy IOP) or a ‘pathological load’ (simulating an elevated glaucomatous IOP) to bovine scleral fibroblasts, as a model of human glaucoma, was conducted to characterise cytoskeletal organisation, chromatin condensation and cell dimensions using immunofluorescence confocal microscopy. Quantification of cell parameters and cytoskeletal element anisotropy were subsequently performed using FibrilTool, and chromatin condensation parameter assessment through a bespoke MATLAB script. The novel findings suggest that physiological load-induced F-actin rearrangement is transient, whereas pathological load, recapitulating in vivo glaucomatous IOP levels, had a reversible and inhibitory influence on remodelling of the cytoskeletal architecture and, further, induction of chromatin condensation. Ultimately, this could compromise cell behaviour. These findings could provide valuable insight into the mechanism(s) used by scleral fibroblasts to mechanically adapt to support biomechanical tissue integrity, and how it could be potentially modified for therapeutic avenues targeting mechanically mediated ocular pathologies such as glaucoma.

•

Physiological strain induced a transient F-actin rearrangement in scleral fibroblasts.

•

In contrast, pathological strain reversibly delayed F-actin rearrangement.

•

Vimentin and β-tubulin networks were largely unaffected by strain regimens.

•

Pathological strain reversibly increased chromatin condensation parameter.

•

Pathological strain may induce ‘inhibition delay’ to confer cytoprotection.

Biomechanical Evaluation of Decellularized and Crosslinked Corneal Implants Manufactured From Porcine Corneas as a Treatment Option for Advanced Keratoconus

Currently corneal transplantation is the main treatment for late-stage keratoconus; however, transplantation procedures are accompanied by significant risk of post-surgical complications; this in addition to supply limitations imposed by a worldwide shortage of human donor corneas, has driven the development of alternative therapies. One such therapy is the use of corneal implants derived from porcine corneas (Xenia^®, Gebauer Medizintechnik GmbH, Neuhausen, DE). In contrast to human donor tissue, these implants can be produced on demand and due to the processes used pose no risks for host-immune rejection. Their use has already been demonstrated clinically in patients for preventing the progression of topographic changes in keratoconus whilst improving visual acuity. The implants are derived from natural tissue and not standardised synthetic material, whilst this likely reduces the risk of issues with bio-incompatibility, there is inevitably variability in their intrinsic mechanical properties which requires investigation. Here, speckle interferometry is employed to examine the biomechanical properties, in response to physiologically representative forces, of native porcine corneal tissue prior to processing and after a proprietary 4-stage process involving decellularization, washing, compression and crosslinking. The control lenticules had an average Young’s modulus (E) of 11.11 MPa (range 8.39–13.41 MPa), following processing average E of the lenticules increased by 127% over that of the unprocessed tissue to 25.23 MPa (range 18.32–32.9 MPa). The variability in E of the lenticules increased significantly after processing suggesting variability in the propensity of the native tissue to processing. In summary, it is possible to produce thin (<90 µm) lenticules from porcine corneas with enhanced stiffness that are effective for treating late-stage keratoconus. Due to the observed variability in the responses of lenticules to processing, interferometry could be a useful technique for ensuring quality control in commercial production via biomechanical screening.

Tropocollagen springs allow collagen fibrils to stretch elastically

The mechanical properties of connective tissues are tailored to their specific function, and changes can lead to dysfunction and pathology. In most mammalian tissues the mechanical environment is governed by the micro- and nano-scale structure of collagen and its interaction with other tissue components, however these hierarchical properties remain poorly understood. In this study we use the human cornea as a model system to characterise and quantify the dominant deformation mechanisms of connective tissue in response to cyclic loads of physiological magnitude. Synchronised biomechanical testing, x-ray scattering and 3D digital image correlation revealed the presence of two dominant mechanisms: collagen fibril elongation due to a largely elastic, spring-like straightening of tropocollagen supramolecular twist, and a more viscous straightening of fibril crimp that gradually increased over successive loading cycles. The distinct mechanical properties of the two mechanisms suggest they have separate roles in vivo. The elastic, spring-like mechanism is fast-acting and likely responds to stresses associated with the cardiac cycle, while the more viscous crimp mechanism will respond to slower processes, such as postural stresses. It is anticipated that these findings will have broad applicability to understanding the normal and pathological functioning of other connective tissues such as skin and blood vessels that exhibit both helical structures and crimp.

Statement of significance

The tropocollagen spring mechanism allows collagen fibrils from some tissues to elongate significantly under small loads, and its recent discovery has the potential to change our fundamental understanding of how tissue deforms. This time-resolved study quantifies the contribution of the spring mechanism to the local strain in stretched tissue and compares it to the contribution associated with the straightening of fibril waviness, the widely accepted primary low-load strain mechanism. The spring mechanism contributed more to the local tissue strain than fibril straightening, and was found to be elastic while fibril straightening was more viscous. The results suggest that the viscoelastic behaviour of a biomaterial is controlled, at least in part, by the relative amount of fibril-scale crimp and tropocollagen supramolecular twist.

The Effect of Axial Length Elongation on Corneal Biomechanical Property

Background: To investigate the correlation between the corneal biomechanical parameter stress-strain index (SSI) and axial length (AL) in moderately elongated eye (MEE) and severely elongated eye (SEE).

Methods: This study included 117 eyes from 117 participants. Among them, 59 (50.4%) had MEE (AL<26 mm) and 58 (49.6%) had SEE (AL≥26 mm). AL was measured using Lenstar LS-900, and central corneal thickness (CCT) and anterior chamber volume (ACV) were measured using Pentacam. SSI was measured via corneal visualisation Scheimpflug technology (Corvis ST). Kolmogorov-Smirnov test, Student’s t-test, and Pearson and partial correlation analyses were used for statistical analyses.

Results: The mean (±SD) SSI was 1.08 ± 0.15 in the MEE group and 0.92 ± 0.13 in the SEE group (p < 0.01). SSI was positively correlated with age (MEE: r = 0.326, p < 0.05; SEE: r = 0.298, p < 0.05) in both groups; it was negatively correlated with AL (r = −0.476, p < 0.001) in the MEE group but not in the SEE group (p > 0.05). CCT was negatively correlated with AL (r = −0.289, p < 0.05) and ACV positively correlated with AL (r = 0.444, p < 0.001) in the MEE group. Neither CCT nor ACV was correlated with AL (p > 0.05) in the SEE group.

Conclusion: Corneal biomechanical parameter SSI, which represents the stiffness of corneal tissue, was lower in the SEE group than in the MEE group. When analyzed separately, SSI was negatively correlated with AL in the MEE group, but not in the SEE group, which may provide insight into different ocular growth patterns between lower myopia and higher myopia.

Developing a Lung Model in the Age of COVID-19: A Digital Image Correlation and Inverse Finite Element Analysis Framework

Pulmonary diseases, driven by pollution, industrial farming, vaping, and the infamous COVID-19 pandemic, lead morbidity and mortality rates worldwide. Computational biomechanical models can enhance predictive capabilities to understand fundamental lung physiology; however, such investigations are hindered by the lung’s complex and hierarchical structure, and the lack of mechanical experiments linking the load-bearing organ-level response to local behaviors. In this study we address these impedances by introducing a novel reduced-order surface model of the lung, combining the response of the intricate bronchial network, parenchymal tissue, and visceral pleura. The inverse finite element analysis (IFEA) framework is developed using 3-D digital image correlation (DIC) from experimentally measured non-contact strains and displacements from an ex-vivo porcine lung specimen for the first time. A custom-designed inflation device is employed to uniquely correlate the multiscale classical pressure-volume bulk breathing measures to local-level deformation topologies and principal expansion directions. Optimal material parameters are found by minimizing the error between experimental and simulation-based lung surface displacement values, using both classes of gradient-based and gradient-free optimization algorithms and by developing an adjoint formulation for efficiency. The heterogeneous and anisotropic characteristics of pulmonary breathing are represented using various hyperelastic continuum formulations to divulge compound material parameters and evaluate the best performing model. While accounting for tissue anisotropy with fibers assumed along medial-lateral direction did not benefit model calibration, allowing for regional material heterogeneity enabled accurate reconstruction of lung deformations when compared to the homogeneous model. The proof-of-concept framework established here can be readily applied to investigate the impact of assorted organ-level ventilation strategies on local pulmonary force and strain distributions, and to further explore how diseased states may alter the load-bearing material behavior of the lung. In the age of a respiratory pandemic, advancing our understanding of lung biomechanics is more pressing than ever before.

Regional variation of corneal stromal deformation measured by high-frequency ultrasound elastography

The cornea's mechanical response to intraocular pressure elevations may alter in ectatic diseases such as keratoconus. Regional variations of mechanical deformation in normal and keratoconus eyes during intraocular pressure elevation have not been well-characterized. We applied a high-frequency ultrasound elastography technique to characterize the regional deformation of normal and keratoconus human corneas through the full thickness of corneal stroma. A cross-section centered at the corneal apex in 11 normal and 2 keratoconus human donor eyes was imaged with high-frequency ultrasound during whole globe inflation from 5 to 30 mmHg. An ultrasound speckle tracking algorithm was used to compute local tissue displacements. Radial, tangential, and shear strains were mapped across the imaged cross-section. Strains in the central (1 mm surrounding apex) and paracentral (1 to 4 mm from apex) regions were analyzed in both normal and keratoconus eyes. Additional regional analysis was performed in the eye with severe keratoconus presenting significant thinning and scarring. Our results showed that in normal corneas, the central region had significantly smaller tangential stretch than the paracentral region, and that within the central region, the magnitudes of radial and shear strains were significantly larger than that of tangential strain. The eye with mild keratoconus had similar shear strain but substantially larger radial strains than normal corneas, while the eye with severe keratoconus had similar overall strains as in normal eyes but marked regional heterogeneity and large strains in the cone region. These findings suggested regional variation of mechanical responses to intraocular pressure elevation in both normal and keratoconus corneas, and keratoconus appeared to be associated with mechanical weakening in the cone region, especially in resisting radial compression. Comprehensive characterization of radial, tangential, and shear strains through corneal stroma may provide new insights to understand the biomechanical alterations in keratoconus.

Effect of different preconditioning protocols on the viscoelastic inflation response of the posterior sclera

Preconditioning by repeated cyclic loads is routinely used in ex vivo mechanical testing of soft biological tissues. The goal of preconditioning is to achieve a steady and repeatable mechanical response and to measure material properties that are representative of the in vivo condition. Preconditioning protocols vary across studies, and their effect on the viscoelastic response of tested soft tissue is typically not reported or analyzed. We propose a methodology to systematically analyze the preconditioning process with application to inflation testing. We investigated the effect of preconditioning on the viscoelastic inflation response of tree shrew posterior sclera using two preconditioning protocols: (i) continuous cyclic loading-unloading without rest and (ii) cyclic loading-unloading with 15-min rest between cycles. Posterior scleral surface strain was measured using three-dimensional Digital Image Correlation (3D-DIC). We used five variables of characterizing features of the stress-strain loop curve to compare the two preconditioning protocols. Our results showed protocol-dependent differences in the tissue response during preconditioning and at the preconditioned state. Incorporating a resting time between preconditioning cycles significantly decreased the number of cycles (10.5 ± 2.9 cycles vs. 3.1 ± 0.5 cycles, p < 0.001) but increased the total time (15.8 ± 4.4 min vs. 51.2 ± 8.3 min, p < 0.001) needed to reach preconditioned state. At the preconditioned state, 2 of 5 characteristic variables differed significantly between protocols: hysteresis loop area (difference=0.023 kJ/m³, p = 0.0020) and elastic modulus at high IOPs (difference=24.0 MPa, p = 0.0238). Our results suggest that the analysis of the preconditioning process is an essential part of inflation experiments and a prerequisite to properly characterize the tissue viscoelastic response. Furthermore, material properties obtained at the preconditioned state can be impacted by the resting time used during preconditioning and may not be directly compared across studies if the resting time varies by 15 min between studies. STATEMENT OF SIGNIFICANCE: Although applying a preconditioning protocol by repeated cyclic loads is common practice in ex vivo mechanical characterization of soft tissues, the tissue response is typically not reported or analyzed, and the protocol's potential effect on the response remains unclear. This is partially caused by lack of a standardized methodology to precondition soft tissues. We present the first systematic analysis of two representative preconditioning protocols used during inflation testing in ocular biomechanics. Our results show protocol-dependent differences in the viscoelastic response during the preconditioning process and at the preconditioned state. Consequently, the analysis of the preconditioning response represents an essential part of mechanical testing and a prerequisite to properly characterize the tissue viscoelastic response. The effect of preconditioning on the preconditioned state response must be considered when comparing results across studies with different preconditioning protocols.

Interferometric Ex Vivo Evaluation of the Spatial Changes to Corneal Biomechanics Introduced by Topographic CXL: A Pilot Study

PURPOSE

To determine the efficacy of interferometry for examining the spatial changes to the corneal biomechanical response to simulated intraocular pressure (IOP) fluctuations that occur after corneal cross-linking (CXL) applied in different topographic locations.

METHODS

Displacement speckle pattern interferometry (DSPI) was used to measure the total anterior surface displacement of human and porcine corneas in response to pressure variations up to 1 mm Hg from a baseline pressure of 16.5 mm Hg, both before and after CXL treatment, which was applied in isolated topographic locations (10-minute riboflavin soak [VibeX-Xtra; Avedro, Inc], 8-minute ultraviolet-A exposure at 15 mW/cm²). Alterations to biomechanics were evaluated by directly comparing the responses before and after treatment for each cornea.

RESULTS

Before CXL, the corneal response to loading indicated spatial variability in mechanical properties. CXL treatments had a variable effect on the corneal response to loading dependent on the location of treatment, with reductions in regional displacement of up to 80% in response to a given pressure increase.

CONCLUSIONS

Selectively cross-linking in different topographic locations introduces position-specific changes to mechanical properties that could potentially be used to alter the refractive power of the cornea. Changes to the biomechanics of the cornea after CXL are complex and appear to vary significantly depending on treatment location and initial biomechanics. Hence, further investigations are required on a larger number of corneas to allow the development of customized treatment protocols. In this study, laser interferometry was demonstrated to be an effective and valuable tool to achieve this. [J Refract Surg. 2021;37(4):263-273.].

A review of corneal biomechanics: Mechanisms for measurement and the implications for refractive surgery

Detailed clinical assessment of corneal biomechanics has the potential to revolutionize the ophthalmic industry through enabling quicker and more proficient diagnosis of corneal disease, safer and more effective surgical treatments, and the provision of customized and optimized care. Despite these wide-ranging benefits, and an outstanding clinical need, the provision of technology capable of the assessment of corneal biomechanics in the clinic is still in its infancy. While laboratory-based technologies have progressed significantly over the past decade, there remain significant gaps in our knowledge regarding corneal biomechanics and how they relate to shape and function, and how they change in disease and after surgical intervention. Here, we discuss the importance, relevance, and challenges associated with the assessment of corneal biomechanics and review the techniques currently available and underdevelopment in both the laboratory and the clinic.

Experimental evaluation of the viscoelasticity of porcine vitreous

This study aims to estimate the material properties of the porcine vitreous while testing it in close to its natural physiological conditions. Eighteen porcine eyes were tested within 48 h post-mortem. A custom-built computer-controlled test rig was designed to support, load and monitor the behaviour of eye globes while being subjected to dynamic rotation cycles mimicking saccade eye movement. Specimens were glued to the base of a container, surrounded by gelatin, frozen and cut in half to expose the vitreous. After thawing, the container was subjected to concentric dynamic rotations of up to 5°, 10° or 15°, while taking 50 MP photos of the specimen every 2 ms. The images were analysed by a digital image correlation algorithm to trace the movement of marked points on the vitreous surface with different radii from the centre of the posterior chamber. The initial camera image was used in building a finite-element model of the test set-up, which was used in an inverse analysis exercise to estimate the material properties of the vitreous. Angular displacements of the monitored points were up to 3.3°, 4.1° and 3.9° in response to eye rotations of 5°, 10° and 15°, respectively. With the experimental relationships between eye rotation and angular displacements used as target behaviour, the inverse analysis exercise estimated the initial shear modulus, the long-term shear modulus and the viscoelastic decay constant of the porcine vitreous as 2.10 ± 0.15 Pa, 0.50 ± 0.04 Pa and 1.20 ± 0.09 s^-1, respectively. Consideration of the viscoelasticity of the vitreous was essential to represent its experimental behaviour. Testing the vitreous in close to its normal physiological conditions produced estimations of the initial shear modulus and long-term shear modulus that were, respectively, smaller and larger than reported values (Zimberlin et al. 2010 Soft Matter 6, 3632-3635. (doi:10.1039/b925407b), Liu et al. 2013 J. Biomech. 46, 1321-7. (doi:10.1016/j.jbiomech.2013.02.006), Rossi et al. 2011 Invest. Ophthalmol. Vis. Sci. 52, 3994-4002. (doi:10.1167/iovs.10-6477)).

Quasi-Static Optical Coherence Elastography to Characterize Human Corneal Biomechanical Properties

Purpose

Quasi-static optical coherence elastography (OCE) is an emerging technology to investigate corneal biomechanical behavior in situations similar to physiological stress conditions. Herein OCE was applied to evaluate previously inaccessible biomechanical characteristics of human corneal tissue and to study the role of Bowman's layer in corneal biomechanics.

Methods

Human corneal donor buttons (n = 23) were obtained and Descemet's membrane and endothelium were removed. In 11 corneas, Bowman's layer was ablated by a 20 µm stromal excimer laser ablation. Buttons were mounted on an artificial anterior chamber and subjected to a pressure modulation from 10 to 30 mm Hg, and back to 10 mm Hg, in steps of 1 mm Hg. At each step, a spectral-domain optical coherence tomography scan was obtained. Displacements were analyzed by optical flow tracking, and strain over the entire stromal depth was retrieved from the phase gradient of the complex interference signal.

Results

During pressure increase, corneal tissue moved upward (486–585 nm/mm Hg) but did not fully recover (Δ= 2.63 to 8.64 µm) after pressure decrease. Vertical corneal strain distribution was negative in the anterior and positive in the posterior cornea, indicating simultaneous corneal compression and expansion, respectively. Bowman's layer caused minor localized differences in corneal strain distribution.

Conclusions

Corneal strain distribution is more complex than previously assumed, with a fundamental difference in mechanical response between the anterior and posterior stroma. Clinically, OCE technology might be used to monitor the progression of corneal ectatic diseases and to determine the success of corneal cross-linking.

An interferometric ex vivo study of corneal biomechanics under physiologically representative loading, highlighting the role of the limbus in pressure compensation

Background

The mechanical properties of the cornea are complex and regionally variable. This paper uses an original method to investigate the biomechanics of the cornea in response to hydrostatic loading over the typical physiological range of intra-ocular pressure (IOP) fluctuations thereby increasing understanding of clinically relevant corneal biomechanical properties and their contributions to the refractive properties of the cornea.

Methods

Displacement speckle pattern interferometry (DSPI) was used to measure the total surface displacement of 40 porcine and 6 human corneal-scleral specimens in response to pressure variations up to 1 mmHg from a baseline of 16.5 mmHg. All specimens were mounted in a modified artificial anterior chamber (AAC) and loaded hydrostatically. Areas of high strain in response to loading were identified by comparing the displacements across different regions.

Results

The nature of the response of the corneal surface to loading demonstrated high regional topographic variation. Mechanical properties were shown to be asymmetrical, and deformation of the limbal and pre-limbal regions dominated these responses respectively with over 90% (N-T) and 60% (S-I) of the total maximum displacement occurring in these regions indicating high-strain. In contrast, the curvature of the central cornea remained relatively unchanged merely translating in position.

Conclusions

The limbal and pre-limbal regions of the cornea appear to be fundamental to the absorption of small pressure fluctuations facilitating the curvature of the central cornea to remain relatively unchanged. The differential mechanical properties of this region could have important implications for the application of corneal surgery and corneal crosslinking, warranting further investigation.

Topical Corneal Cross-Linking Solution Delivered Via Corneal Reservoir in Dutch-Belted Rabbits

Purpose

A topical corneal cross-linking solution that can be used as an adjunct or replacement to standard photochemical cross-linking (UV-riboflavin) methods remain an attractive possibility. Optimal concentration and delivery method for such topical corneal stabilization in the living rabbit eye were developed.

Methods

A series of experiments were carried out using Dutch-belted rabbits (3 months old, weighing 1.0–1.5 kg) and topical cross-linking solutions (sodium hydroxymethylglycinate) (10–250 mM) delivered via corneal reservoir. The application regimen included a one-time 30-minute application (10–40 mM sodium hydroxymethylglycinate) as well as a once per week 5-minute application (250 mM sodium hydroxymethylglycinate) for 7 weeks. Animals were evaluated serially for changes in IOP, pachymetry, epithelial integrity, and endothelial cell counts. Keratocyte changes were identified using intravital laser scanning confocal microscopy. Post mortem efficacy was evaluated by mechanical inflation testing.

Results

Overall, there were very few differences observed in right eye treated versus left eye controls with respect to intraocular pressure, pachymetry, and endothelial cell counts, although 30-minute cross-linking techniques did cause transient increases in thickness resolving within 7 days. Epithelial damage was noted in all of the 30-minute applications and fully resolved within 72 hours. Keratocyte changes were significant, showing a wound healing pattern similar to that after riboflavin UVA photochemical cross-linking in rabbits and humans. Surprisingly, post mortem inflation testing showed that the lower concentration of 20 mM delivered over 30 minutes showed the most profound stiffening/strengthening effect.

Conclusions

Topical cross-linking conditions that are safe and can increase corneal stiffness/strength in the living rabbit eye have been identified.

Translational Relevance

A topical corneal cross-linking solution delivered via corneal reservoir is shown to be both safe and effective at increasing tissue strength in living rabbit eyes and could now be tested in patients suffering from keratoconus and other conditions marked by corneal tissue weakness.

Full-field displacement measurement of corneoscleral shells by combining multi-camera speckle interferometry with 3D shape reconstruction

Changes in the biomechanical properties of the connective tissue of the eye occur with age and underlie the development of several ocular diseases, such as glaucoma, myopia, and keratoconus. The biomechanical dynamics of ocular connective tissue are measured by ex vivo inflation testing, in which intraocular pressure (IOP) is varied and optical methods are used to produce maps of corneal and scleral displacement. Current optical methods are limited by acquisition rate, occlusions, poor spatial resolution, and insufficient 3D mapping. We developed an interferometric optical method integrates four-camera electronic speckle pattern interferometry (ESPI) and a novel three-dimensional (3D) shape reconstruction process to measure shape and full-field mechanical deformations of corneal and scleral shells during ex vivo inflation testing. Each camera provides accurate measurements of the laser beam phase related to deformations of the specimen surface; a multi-view stereovision method generates the shape of the specimen and a functional form that links every pixel of a given camera to 3D points on the specimen's visible surface. In this way, dynamic deformations of the specimen are localized, with quantification of the time-dependent 3D displacements of the specimen at nanometric accuracy. The ESPI-3D system is suitable for analyzing scleral deformation and morphological changes caused by time-varying IOP.

Inflation experiments and inverse finite element modelling of posterior human sclera

The complexity of inverse finite element modelling methods used in ocular biomechanics research has significantly increased in recent years in order to produce material parameters that capture microscale tissue behaviour. This study presents a more accessible method for researchers to optimise sclera material parameters for use in finite element studies where macroscale sclera displacements are required. Five human donor sclerae aged between 36 and 72 years were subjected to cycles of internal pressure up to 61 mmHg using a custom-built inflation rig. Displacements were measured using a laser beam and two cameras through a digital image correlation algorithm. Specimen-specific finite element models incorporating regional thickness variation and sclera surface topography were divided into six circumferential regions. An inverse finite element procedure was used to optimise Ogden material parameters for each region. The maximum root mean squared (RMS) error between the numerical and experimental displacements within individual specimens was 17.5 µm. The optimised material parameters indicate a gradual reduction in material stiffness (as measured by the tangent modulus) from the equator to the posterior region at low-stress levels up to 0.005 MPa. The variation in stiffness between adjacent regions became gradually less apparent and statistically insignificant at higher stresses. The study demonstrated how inflation testing combined with inverse modelling could be used to effectively characterise regional material properties capable of reproducing global sclera displacements. The material properties were found to vary between specimens, and it is expected that age could be a contributing factor behind this variation.

Cornea Full-field Displacement and Strain Measurement in Vivo Using Three-dimensional Digital Image Correlation

SIGNIFICANCE

Corneal biomechanics are becoming increasingly important for clinical evaluation and diagnosis, such as in refractive surgery, glaucoma, and keratoconus. We developed and used a new technique to measure the full-field displacement and strain of the cornea under elevated intraocular pressure in vivo.

PURPOSE

In this study, we propose a three-dimensional digital image correlation method for the determination of corneal biomechanical properties.

METHODS

Corneal deformation change on the rabbit eyeball was investigated under different inflation conditions with intraocular pressure levels of 2.3, 3.3, 5.3, and 6.9 kPa in vivo.

RESULTS

The proposed method was able to measure the displacement of the corneal surface and to construct real-time full-field three-dimensional deformation vector fields and strain mapping. The results show that the strain distribution is not uniform on the corneal surface at each intraocular pressure level. Compression strains, rather than extensile strains, are dominant at the corneal apex region.

CONCLUSIONS

This technique has the potential to be used as an assistive tool for the determination of corneal biomechanical properties in ophthalmologic investigations.

A Chest-Laminated Ultrathin and Stretchable E-Tattoo for the Measurement of Electrocardiogram, Seismocardiogram, and Cardiac Time Intervals

Seismocardiography (SCG) is a measure of chest vibration associated with heartbeats. While skin soft electronic tattoos (e-tattoos) have been widely reported for electrocardiogram (ECG) sensing, wearable SCG sensors are still based on either rigid accelerometers or non-stretchable piezoelectric membranes. This work reports an ultrathin and stretchable SCG sensing e-tattoo based on the filamentary serpentine mesh of 28-µm-thick piezoelectric polymer, polyvinylidene fluoride (PVDF). 3D digital image correlation (DIC) is used to map chest vibration to identify the best location to mount the e-tattoo and to investigate the effects of substrate stiffness. As piezoelectric sensors easily suffer from motion artifacts, motion artifacts are effectively reduced by performing subtraction between a pair of identical SCG tattoos placed adjacent to each other. Integrating the soft SCG sensor with a pair of soft gold electrodes on a single e-tattoo platform forms a soft electro-mechano-acoustic cardiovascular (EMAC) sensing tattoo, which can perform synchronous ECG and SCG measurements and extract various cardiac time intervals including systolic time interval (STI). Using the EMAC tattoo, strong correlations between STI and the systolic/diastolic blood pressures, are found, which may provide a simple way to estimate blood pressure continuously and noninvasively using one chest-mounted e-tattoo.

Line-Field Optical Coherence Tomography as a tool for In vitro characterization of corneal biomechanics under physiological pressures

There has been a lot of interest in accurately characterising corneal biomechanical properties under intraocular pressure (IOP) to help better understand ocular pathologies that are associated with elevated IOP. This study investigates the novel use of Line-Field Optical Coherence Tomography (LF-OCT) as an elastographic tool for accurately measuring mechanical properties of porcine corneas based on volumetric deformation following varying IOPs. A custom-built LF-OCT was used to measure geometrical and corneal surface displacement changes in porcine corneas under a range of IOPs, from 0-60 mmHg. Corneal thickness, elastic properties and hysteresis were calculated as a function of pressure. In addition, the effects of hydration were explored. We found that the elastic modulus increased in a linear fashion with IOP. Corneal thickness was found to reduce with IOP, decreasing 14% from 0 to 60 mmHg. Prolonged hydration in phosphate buffered saline (PBS) was found to significantly increase the elastic modulus and corneal hysteresis. Our study demonstrates that LF-OCT can be used to accurately measure the elastic properties based on volumetric deformation following physiological pressures. Furthermore, we show that prolonged hydration in PBS has a significant effect on the measured corneal properties.

Corneal elevation changes after forced eyelid closure in healthy participants and in patients with keratoconus

BACKGROUND

To assess pointwise corneal elevation changes after forced eyelid closure test (FECT) in the eyes of healthy subjects and in eyes with keratoconus.

METHODS

Twenty-nine subjects with keratoconus and 31 healthy volunteers were evaluated. Patients with keratoconus who had corneal hydrops, apical scarring, corneal thickness ≤ 400 μm, ocular surface disease, contact lens wear on the examination day and a history of corneal cross-linking were excluded. Exclusion criteria for healthy participants were spherical error > +3.00 D and < -3.00 D, corneal astigmatism > 1.50 D, corneal curvature > 47 D, ocular allergy, clinical findings and family history of keratoconus. Pentacam was performed before and after 20 seconds of FECT and raw data were extracted from the built-in software. Pointwise anterior and posterior elevation changes in the central 8 mm cornea were assessed using paired samples t-test and heat maps were constructed to reflect mean changes and statistically significant data points. Statistical significance was assumed at p < 0.01.

RESULTS

Age and gender were similar between healthy subjects (24.5 ± 1.6 years, 46.4 per cent female) and subjects with keratoconus (28.6 ± 9.2 years, 46.4 per cent female, p = 0.19, 0.61, respectively). Healthy eyes displayed posterior depression clustering in the inferotemporal and inferonasal areas (mean change: -4.5 ± 7.8 μm and -5.2 ± 9.8 μm, respectively, all p < 0.01). In contrast, keratoconus eyes exhibited a wider area of posterior elevation clustering in the inferior cornea (mean change: 8.1 ± 14.5 μm, all p < 0.01) with a small extension in the inferotemporal cornea (mean change: 12.1 ± 22.3 μm, all p < 0.01).

CONCLUSION

FECT elicits corneal elevation changes mainly in the inferior cornea with the change being more pronounced and wider in eyes with keratoconus.

Optical elastography: tracking surface waves with digital image correlation

Elastography consists in evaluating the propagation speed of waves into a tissue to estimate its stiffness. Usually this method is based on Ultrasounds, magnetic resonance imaging or optical coherent tomography. This paper proposes a simple optic method using ultrafast cameras. Based on digital image correlation (DIC), the tracking of elastic surface wave from white light intensity pattern, allows estimating the propagation speed as an indicator of the tissue local stiffness. Two configurations are presented: (1) 2D imaging of a flat phantom surface with a single camera and (2) 3D imaging of a curved phantom surface with two cameras. As a feasibility study of the first configuration, surface wave speed was measured on isotropic and anisotropic phantoms. Comparisons with ultrasound methods fully validate this approach. Although more sophisticated, the second configuration account for propagation distortions caused by locally curved topology. Triangulation techniques used to retrieve local topology are named stereo-correlation in the field of biomechanics. Stereo-elastography is thus proposed to determine tissue local elasticity from any soft tissue surface wave.

STUDY OF THE TRANSVERSAL DEFORMATION OF CORNEAL STRIP UNDER UNIAXIAL LOADING

Uniaxial test is easy to access and to obtain accuracy data, but it is difficult to acquire two-dimensional deformation information. We investigated the relationship between the two strain components of corneal strip in uniaxial tests, which is the basis for determining of anisotropic strain energy function of cornea via uniaxial tests. Nine rabbits were taken. The left and right corneas were cut along superior-inferior (SI) and nasal-temporal (NT) direction, respectively. For each strip the uniaxial test was carried out, and the tensile displacements, strip images and loads were recorded. Then the stretching strain, the transversal strains and stress were obtained. Optimization based inverse analysis was utilized to find the best among six fitting models that characterizes the relationship between two strain components in uniaxial tests. All models fitted well the experimental data gathered for corneal strips ([Formula: see text]). According to the model selection index, the power model achieved the best performance index: 0.1268 for SI strips and 0.1063 for NT strips versus 0.151 (SI strips) and 0.107 (NT strips) found at most by other models. Thus, it is the most suitable one for describing the relationship between the two strain components of corneal strip during uniaxial stretching.

A comparative study of the elastic fibre system within the mouse and human cornea

The cornea relies on its organised extracellular matrix for maintaining transparency and biomechanical strength. Studies have identified an elastic fibre system within the human posterior cornea, thought to allow for slight deformations in response to internal pressure fluctuations within the eye. However, the type of elastic fibres that exist within the cornea and their roles remain elusive. The aim of this study was to compare the distribution and organisation of the elastic fibres within the posterior peripheral mouse and human cornea, and elucidate how these fibres integrate with the trabecular meshwork, whilst characterising the distribution of their main likely components (fibrillin-1, elastin and type VI collagen) in different parts of the cornea and adjacent sclera.

We identified key differences in the elastic fibre system between the human and mouse cornea. True elastic fibres (containing elastin) were identified within the human posterior peripheral cornea. Elastic fibres appeared to present as an extensive network throughout the mouse corneal stroma, but as fibrillin-rich microfibril bundles rather than true elastic fibres. However, tropoelastin staining indicated the possibility that true elastic fibres had yet to develop in the young mice studied. Differences were also apparent within the anatomy of the trabecular meshwork. The human trabecular meshwork appeared to insert between the corneal stroma and Descemet's membrane, with elastic fibres continuing into the stroma from the trabecular meshwork anterior to Descemet's membrane. Within the mouse cornea, no clear insertion point of the trabecular meshwork was seen, instead the elastic fibres within the trabecular meshwork continued into Descemet's membrane, with the trabecular meshwork joining posterior to Descemet's membrane.

•

True Elastic fibres (containing elastin) occur within the human posterior peripheral cornea.

•

A fibrillin-rich microfibril system was found throughout the mouse corneal stroma, with no apparent true elastic fibres.

•

Human trabecular meshwork inserts between the posterior corneal stroma and Descemet's membrane.

•

A continuity of elastic fibres was identified between cornea and trabecular meshwork in both species.

The influence of hydration on different mechanical moduli of the cornea

PURPOSE

To determine the interrelation of different elastic moduli of the cornea and to investigate their dependency on corneal hydration.

METHODS

Rabbit eyes were divided into four groups. Corneas were excised and mounted into a Barron artificial anterior chamber. Various corneal hydration steady states were achieved with different dextran T-500 concentrations in the anterior chamber, as well as on the corneal anterior surface. The treatment-solutions of each group contained either 5, 10, 15, or 20% w/w dextran. Ultrasound pachymetry was used to measure central corneal thickness. Brillouin microscopy of the central cornea determined the longitudinal bulk modulus by means of Brillouin frequency shift. Subsequently, a 5-mm-wide central strip was taken for extensiometry to measure the tangential elastic modulus.

RESULTS

The longitudinal bulk modulus was 1.2-times higher in corneas dehydrated with 20% dextran compared to those hydrated with 5% dextran. In contrast, the tangential elastic modulus increased by 4.4 times. The obtained longitudinal bulk moduli were two orders of magnitude bigger than the tangential elastic moduli. Regression analysis of longitudinal bulk modulus and tangential elastic modulus revealed a quadratic relation. The bulk modulus seemed to be independent of tension, whereas the elastic modulus was tension-dependent. Greater corneal hydration led to significantly thicker pachymetry.

CONCLUSION

Corneal biomechanics are highly dependent on the level of corneal hydration. Surprisingly, tangential elastic moduli were more sensitive to hydration changes than longitudinal bulk moduli. A quadratic relation was found between both moduli.

Corneal Biomechanical Properties after FS-LASIK with Residual Bed Thickness Less Than 50% of the Original Corneal Thickness

Background

The changes in corneal biomechanical properties after LASIK remain an unknown but important topic for surgical design and prognostic evaluation. This study aims to observe the postoperative corneal biomechanical properties one month after LASIK with amount of corneal cutting (ACC) greater than 50% of the central corneal thickness (CCT).

Methods

FS-LASIK was performed in 10 left rabbit eyes with ACC being 60% (L60) and 65% (L65) of the CCT, while the right eyes (R) were the control. After 4 weeks, rabbits were executed and corneal strip samples were prepared for uniaxial tensile tests.

Results

At the same strain, the stresses of L65 and L60 were larger than those of R. The elastic moduli of L60 and L65 were larger than those of R when the stress was 0.02 MPa, while they began to be less than those of R when stress exceeds the low-stress region. After 10 s relaxation, the stress of specimens L65, L60, and R increased in turn.

Conclusion

The elastic moduli of the cornea after FS-LASIK with ACC greater than 50% of the CCT do not become less under normal rabbit IOP. The limit stress grows with the rise of ACC when relaxation becomes stable.

Imaging Corneal Biomechanical Responses to Ocular Pulse Using High-Frequency Ultrasound

Imaging corneal biomechanical changes or abnormalities is important for better clinical diagnosis and treatment of corneal diseases. We propose a novel ultrasound-based method, called ocular pulse elastography (OPE), to image corneal deformation during the naturally occurring ocular pulse. Experiments on animal and human donor eyes, as well as synthetic radiofrequency (RF) data, were used to evaluate the efficacy of the OPE method. Using very high-frequency ultrasound (center frequency = 55 MHz), correlation-based speckle tracking yielded an accuracy of less than 10% error for axial tissue displacements of or above. Satisfactory speckle tracking was achieved for out-of-plane displacements up to . Using synthetic RF data with or without a pre-defined uniform strain, the OPE method detected strains down to 0.0001 axially and 0.00025 laterally with an error less than 10%. Experiments in human donor eyes showed excellent repeatability with an intraclass correlation of 0.98. The measurement outcome from OPE was also shown to be highly correlated with that of standard inflation. These results suggest the feasibility of OPE as a potential clinical tool for evaluating corneal biomechanics in vivo.

Assessment of Corneal Biomechanical Properties with Inflation Test Using Optical Coherence Tomography

Biomechanical properties are important for the cornea to maintain its normal shape and function. There is still a need to develop better methods for accurate measurement of corneal mechanical properties. In this study, we propose to introduce the optical coherence tomography (OCT) in inflation test for the imaging of corneal deformation in order to measure its biomechanical properties. Ten cornea-mimicking silicone phantoms with different stiffness and five fresh porcine corneas were tested using the proposed method. Intra-ocular pressure was changed from 10 to 90 mmHg using two different loading rates to observe the pressure-apex displacement relationship and calculate the apparent stiffness of the corneas. Stiffness of the corneal phantoms obtained by the inflation test ranged from 0.2 to 1 MPa, which was highly consistent with the results from the mechanical tensile test (y = 0.70x, p < 0.001). The porcine corneas showed highly viscoelastic behavior with obvious hysteresis in inflation. The apparent stiffness of the porcine corneas was 0.63 ± 0.07 and 1.05 ± 0.08 MPa with loading rates of 3.3 and 33 mmHg/min, respectively. Mapping of corneal surface displacement was also generated for both the phantom and porcine corneas. This study showed that it is feasible to incorporate the high resolution OCT imaging in inflation test to measure the biomechanical properties of the cornea.

The hierarchical response of human corneal collagen to load

Fibrillar collagen in the human cornea is integral to its function as a transparent lens of precise curvature, and its arrangement is now well-characterised in the literature. While there has been considerable effort to incorporate fibrillar architecture into mechanical models of the cornea, the mechanical response of corneal collagen to small applied loads is not well understood. In this study the fibrillar and molecular response to tensile load was quantified using small and wide angle X-ray scattering (SAXS/WAXS), and digital image correlation (DIC) photography was used to calculate the local strain field that gave rise to the hierarchical changes. A molecular scattering model was used to calculate the tropocollagen tilt relative to the fibril axis and changes associated with applied strain. Changes were measured in the D-period, molecular tilt and the orientation and spacing of the fibrillar and molecular networks. These measurements were summarised into hierarchical deformation mechanisms, which were found to contribute at varying strains. The change in molecular tilt is indicative of a sub-fibrillar “spring-like” deformation mechanism, which was found to account for most of the applied strain under physiological and near-physiological loads. This deformation mechanism may play an important functional role in tissues rich in fibrils of high helical tilt, such as skin and cartilage.

Statement of Significance

Collagen is the primary mediator of soft tissue biomechanics, and variations in its hierarchical structure convey the varying amounts of structural support necessary for organs to function normally. Here we have examined the structural response of corneal collagen to tensile load using X-rays to probe hierarchies ranging from molecular to fibrillar. We found a previously unreported deformation mechanism whereby molecules, which are helically arranged relative to the axis of their fibril, change in tilt akin to the manner in which a spring stretches. This “spring-like” mechanism accounts for a significant portion of the applied deformation at low strains (<3%). These findings will inform the future design of collagen-based artificial corneas being developed to address world-wide shortages of corneal donor tissue.

In vivo measurement of regional corneal tangent modulus

Currently available clinical devices are unable to measure corneal biomechanics other than at the central region. Corneal stiffness (S), thickness, and radius of curvature was measured at the central cornea (primary fixation) and 3 mm from the temporal limbus (primary and nasal fixations). The corneal tangent modulus (E) of 25 healthy subjects was calculated from these data. After confirming normality, repeated measures analysis of variance (RMANOVA) revealed significant difference in S (F(2, 48) = 21.36, p < 0.001) at different corneal regions and direction of fixations. E also varied significantly at different corneal regions and direction of fixations (RMANOVA: F(2, 48) = 23.06, p < 0.001). A higher S and a lower E were observed at the temporal region compared with the corneal centre. Nasal fixation further increased S and E values compared with primary fixation. Due to the specific arrangement of corneal collagen fibrils, heterogeneity of corneal biomechanical properties is expected. In future clinical practice, localized corneal biomechanical alternation and measurement might assist corneal disease detection and post-surgery management. In addition, practitioners should be aware of the fixation effect on corneal biomechanical measurement.

3D Characterization of corneal deformation using ultrasound speckle tracking

The three-dimensional (3D) mechanical response of the cornea to intraocular pressure (IOP) elevation has not been previously reported. In this study, we use an ultrasound speckle tracking technique to measure the 3D displacements and strains within the central 5.5 mm of porcine corneas during the whole globe inflation. Inflation tests were performed on dextran-treated corneas (treated with a 10% dextran solution) and untreated corneas. The dextran-treated corneas showed an inflation response expected of a thin spherical shell, with through-thickness thinning and in-plane stretch, although the strain magnitudes exhibited a heterogeneous spatial distribution from the central to more peripheral cornea. The untreated eyes demonstrated a response consistent with swelling during experimentation, with through-thickness expansion overriding the inflation response. The average volume ratios obtained in both groups was near 1 confirming general incompressibility, but local regions of volume loss or expansion were observed. These results suggest that biomechanical measurements in 3D provide important new insight to understand the mechanical response of ocular tissues such as the cornea.

Racioethnic Differences in Human Posterior Scleral and Optic Nerve Stump Deformation

Purpose

The purpose of this study was to quantify the biomechanical response of human posterior ocular tissues from donors of various racioethnic groups to better understand how differences in these properties may play a role in the racioethnic health disparities known to exist in glaucoma.

Methods

Sequential digital image correlation (S-DIC) was used to measure the pressure-induced surface deformations of 23 normal human posterior poles from three racioethnic groups: African descent (AD), European descent (ED), and Hispanic ethnicity (HIS). Regional in-plane principal strains were compared across three zones: the optic nerve stump (ONS), the peripapillary (PP) sclera, and non-PP sclera.

Results

The PP scleral tensile strains were found to be lower for ED eyes compared with AD and HIS eyes at 15 mm Hg (P = 0.024 and 0.039, respectively). The mean compressive strains were significantly higher for AD eyes compared with ED eyes at 15 mm Hg (P = 0.018). We also found that the relationship between tensile strain and pressure was significant for those of ED and HIS eyes (P < 0.001 and P = 0.004, respectively), whereas it was not significant for those of AD (P = 0.392).

Conclusions

Our results suggest that, assuming glaucomatous nerve loss is caused by mechanical strains in the vicinity of the optic nerve head, the mechanism of increased glaucoma prevalence may be different in those of AD versus HIS. Our ONS strain analysis also suggested that it may be important to account for ONS geometry and material properties in future scleral biomechanical analysis.

Corneoscleral stiffening increases IOP spike magnitudes during rapid microvolumetric change in the eye

Factors governing the steady-state IOP have been extensively studied; however, the dynamic aspects of IOP are less understood. Clinical studies have suggested that intraocular pressure (IOP) fluctuation may be associated with glaucoma risk. This study aims to investigate how stiffening of corneoscleral biomechanical properties affects IOP spikes induced by rapid microvolumetric change. Porcine eyes (n = 25 in total) were subjected to volumetric infusions before and after external treatment of a circular area (11 mm diameter) in either the central cornea or posterior sclera. The treated area in the control group was immersed in phosphate-buffered saline (PBS) for 40 min, while the treated area of the chemical crosslinking group was immersed in 4% glutaraldehyde/PBS for 40 min. A subset of the sham-treated eyes was also subjected to volumetric infusions at a raised steady-state IOP. The magnitude of IOP spikes increased after localized chemical crosslinking of either the cornea (27.5% increase, p < 0.001) or the sclera (14.3% increase, p < 0.001) with corneal crosslinking having a stronger effect than scleral crosslinking (p = 0.018). We also observed that raising the steady-state IOP from 15 to 25 mmHg resulted in marked increase in IOP spike magnitudes by 63.9% (p < 0.001). These results suggested that an increased corneoscleral stiffness could significantly increase IOP spike magnitudes at the same volumetric change. Corneal stiffness appeared to have a strong impact on the IOP spike magnitude and may play a major role in regulating rapid volume-pressure dynamics. An increase in steady-state IOP also resulted in larger IOP fluctuations due to the increased "apparent" stiffness of the ocular shell, suggesting a potential interaction between the magnitude of IOP and its fluctuations. Corneoscleral properties may represent additional pathways for understanding and managing glaucoma risk and warrant future investigation.

In-vivo corneal pulsation in relation to in-vivo intraocular pressure and corneal biomechanics assessed in-vitro. An animal pilot study

The aim was to ascertain whether the characteristics of the corneal pulse (CP) measured in-vivo in a rabbit eye change after short-term artificial increase of intraocular pressure (IOP) and whether they correlate with corneal biomechanics assessed in-vitro. Eight New Zealand white rabbits were included in this study and were anesthetized. In-vivo experiments included simultaneous measurements of the CP signal, registered with a non-contact method, IOP, intra-arterial blood pressure, and blood pulse (BPL), at the baseline and short-term elevated IOP. Afterwards, thickness of post-mortem corneas was determined and then uniaxial tensile tests were conducted leading to estimates of their Young's modulus (E). At the baseline IOP, backward stepwise regression analyses were performed in which successively the ocular biomechanical, biometric and cardiovascular predictors were separately taken into account. Results of the analysis revealed that the 3rd CP harmonic can be statistically significantly predicted by E and central corneal thickness (Models: R² = 0.662, p < 0.005 and R² = 0.832, p < 0.001 for the signal amplitude and power, respectively). The 1st CP harmonic can be statistically significantly predicted by the amplitude and power of the 1st BPL harmonic (Models: R² = 0.534, p = 0.015 and R² = 0.509, p < 0.018, respectively). For elevated IOP, non-parametric analysis indicated significant differences for the power of the 1st CP harmonic (Kruskal-Wallis test; p = 0.031) and for the mean, systolic and diastolic blood pressures (p = 0.025, p = 0.019, p = 0.033, respectively). In conclusion, for the first time, the association between parameters of the CP signal in-vivo and corneal biomechanics in-vitro was confirmed. In particular, spectral analysis revealed that higher amplitude and power of the 3rd CP harmonic indicates higher corneal stiffness, while the 1st CP harmonic correlates positively with the corresponding harmonic of the BPL signal.

Plastic modification of the cornea by pneumatic force corrects myopia: Pneumatic keratology

PurposeThe eyeball and, in particular, the cornea deform in vitro by the application of a distending force like other visco-elastic tissue. If the force is large enough, the cornea strains beyond the elastic range, and a permanent deformation occurs. Such permanent strain is referred to as 'plastic' strain. The phenomenon, however, has never been observed or produced on living tissue. This report seeks to demonstrate that the central radius of a patient's cornea can be altered in a controlled manner designed to correct refractive errors.MethodsTo plastically deform the living cornea, we applied a vacuum to the cornea of eight rabbits and five human eyes with a novel device. This device consists of a chamber of 11 mm in diameter. The chamber is radially divided into four interconnected sub-chambers.ResultsHere we show that a strain can be achieved in vivo with a force produced by the application of the specially designed chamber where air is evacuated. An anatomical modification of the cornea of humans and rabbits was achieved. The deformation of the cornea was plastic, and therefore permanent.ConclusionsThe method described here-Pneumatic Keratology- can be used to alter the cornea by non-invasive means. A vacuum chamber with radial openings alters the collagen fibers in the stroma and flattens the cornea. A flatter cornea corrects or reduces myopia.

Elastic microfibril distribution in the cornea: Differences between normal and keratoconic stroma

The optical and biomechanical properties of the cornea are largely governed by the collagen-rich stroma, a layer that represents approximately 90% of the total thickness. Within the stroma, the specific arrangement of superimposed lamellae provides the tissue with tensile strength, whilst the spatial arrangement of individual collagen fibrils within the lamellae confers transparency. In keratoconus, this precise stromal arrangement is lost, resulting in ectasia and visual impairment. In the normal cornea, we previously characterised the three-dimensional arrangement of an elastic fiber network spanning the posterior stroma from limbus-to-limbus. In the peripheral cornea/limbus there are elastin-containing sheets or broad fibers, most of which become microfibril bundles (MBs) with little or no elastin component when reaching the central cornea. The purpose of the current study was to compare this network with the elastic fiber distribution in post-surgical keratoconic corneal buttons, using serial block face scanning electron microscopy and transmission electron microscopy. We have demonstrated that the MB distribution is very different in keratoconus. MBs are absent from a region of stroma anterior to Descemet's membrane, an area that is densely populated in normal cornea, whilst being concentrated below the epithelium, an area in which they are absent in normal cornea. We contend that these latter microfibrils are produced as a biomechanical response to provide additional strength to the anterior stroma in order to prevent tissue rupture at the apex of the cone. A lack of MBs anterior to Descemet's membrane in keratoconus would alter the biomechanical properties of the tissue, potentially contributing to the pathogenesis of the disease.

Effect of topical antiinflammatory drugs on mechanical behavior of rabbit cornea

Background

A variety of antiinflammatory therapies are employed to promote corneal wound healing. The effects of steroidal and nonsteroidal antiinflammatory drugs on the biomechanical properties of rabbit cornea were investigated over time using tensile tests.

Methods

Full-thickness incisions were made and used to analyze the effects of dexamethasone sodium phosphate 0.1% and diclofenac sodium 0.1% on corneal biomechanical properties during wound healing at 7, 14 and 21 days after surgery.

Results

The full-thickness incision deeply modified all of the mechanical properties. At 3 weeks after incision, regardless of the drug therapy, the tensile modulus was about 70% of the value for the intact cornea.

Conclusions

Topical treatment with dexamethasone was particularly effective during the first week after surgery; the second week after surgery, a similar result was observed in the corneas treated with diclofenac. Low doses of steroidal and nonsteroidal antiinflammatory drugs would seem to have the potential to improve biomechanical properties only during the early stage of the healing process of the cornea.