Changes and quantitative characterization of hyper-viscoelastic biomechanical properties for young corneal stroma after standard corneal cross-linking treatment with different ultraviolet-A energies

Biomechanical and histological changes associated with riboflavin ultraviolet-A-induced CXL with different irradiances in young human corneal stroma

Keratoconus (KC) is a degenerative condition affecting the cornea, characterized by progressive thinning and bulging, which can ultimately result in serious visual impairment. The onset and progression of KC are closely tied to the gradual weakening of the cornea's biomechanical properties. KC progression can be prevented with corneal cross-linking (CXL), but this treatment has shortcomings, and evaluating its tissue stiffening effect is important for determining its efficacy. In this field, the shortage of human corneas has made it necessary for most previous studies to rely on animal corneas, which have different microstructure and may be affected differently from human corneas. In this research, we have used the lenticules obtained through small incision lenticule extraction (SMILE) surgeries as a source of human tissue to assess CXL. And to further improve the results' reliability, we used inflation testing, personalized finite element modeling, numerical optimization and histology microstructure analysis. These methods enabled determining the biomechanical and histological effects of CXL protocols involving different irradiation intensities of 3, 9, 18, and 30 mW/cm2, all delivering the same total energy dose of 5.4 J/cm2. The results showed that the CXL effect did not vary significantly with protocols using 3-18 mW/cm2 irradiance, but there was a significant efficacy drop with 30 mW/cm2 irradiance. This study validated the updated algorithm and provided guidance for corneal lenticule reuse and the effects of different CXL protocols on the biomechanical properties of the human corneal stroma.

Effect of rapid cooling, frozen storage, and thawing on the passive viscoelastic properties and structure of the rat aorta

Biological tissues decay over time after harvesting, which alters their biomechanical properties. This poses logistical challenges for studies investigating passive arterial biomechanics as tissues need to be characterized shortly after excision. Freezing and cryopreservation methods can help alleviate the need for biomechanical testing of fresh tissue in human ex vivo studies. However, these methods tend to eliminate or reduce arterial cell functionality and affect passive biomechanics. Furthermore, their impact on dynamic arterial biomechanics remains unknown despite arterial viscoelastic properties being an integral component contributing to arterial stiffness under in vivo loading conditions. The present study aims to investigate the impact of rapid cooling and subsequent storage at -80 °C on the passive viscoelastic properties of arterial tissue and aid in ascertaining whether this is a suitable method to delay tissue analysis for studies investigating passive arterial biomechanics. Control and frozen abdominal rat aorta segments were quasi-statically and dynamically tested using a biaxial testing set-up. The results were modeled using a constituent-based quasi-linear viscoelastic modeling framework, yielding directional stiffness parameters, individual constituent biomechanical contributions, and a quantification of viscoelastic stiffening under dynamic pressurization conditions. Frozen samples displayed significantly decreased wall thickness, viscoelastic dissipation, viscoelastic stiffening, and significantly decreased circumferential deformation with changes in luminal pressure. Furthermore, frozen samples displayed significantly increased circumferential stiffness, pulse wave velocity, and collagen load bearing. Consequently, these changes should be considered when utilizing this tissue preservation method to delay biomechanical characterization of rat aortic tissue.

Simultaneous tensile and shear measurement of the human cornea in vivo using S0- and A0-wave optical coherence elastography

Understanding corneal stiffness is valuable for improving refractive surgery, detecting corneal abnormalities, and assessing intraocular pressure. However, accurately measuring the elastic properties, specifically the tensile and shear moduli that govern mechanical deformation, has been challenging. To tackle this issue, we have developed guided-wave optical coherence elastography that can simultaneously excite and analyze symmetric (S0) and anti-symmetric (A0) elastic waves in the cornea at around 10 kHz frequencies, enabling us to extract tensile and shear properties from measured wave dispersion curves. We verified the technique using elastomer phantoms and ex vivo porcine corneas and investigated the dependence on intraocular pressure using acoustoelastic theory that incorporates corneal tension and a nonlinear constitutive tissue model. In a pilot study involving six healthy human subjects aged 31 to 62, we measured shear moduli (Gzx) of 94±20 kPa (mean±standard deviation) and tensile moduli (Exx) of 4.0±1.1 MPa at central corneas. Our preliminary analysis of age-dependence revealed contrasting trends: -8.3±4.5 kPa/decade for shear and 0.30±0.21 MPa/decade for tensile modulus. This OCE technique has the potential to become a highly useful clinical tool for the quantitative biomechanical assessment of the cornea. STATEMENT OF SIGNIFICANCE: This article reports an innovative elastography technique using two guided elastic waves, demonstrating the measurement of both tensile and shear moduli in human cornea in vivo with unprecedented precision. This technique paves the way for comprehensive investigations into corneal mechanics and holds clinical significance in various aspects of corneal health and disease management.

Optomechanical assessment of photorefractive corneal cross-linking via optical coherence elastography

Purpose: Corneal cross-linking (CXL) has recently been used with promising results to positively affect corneal refractive power in the treatment of hyperopia and mild myopia. However, understanding and predicting the optomechanical changes induced by this procedure are challenging. Methods: We applied ambient pressure modulation based optical coherence elastography (OCE) to quantify the refractive and mechanical effects of patterned CXL and their relationship to energy delivered during the treatment on porcine corneas. Three different patterned treatments were performed, designed according to Zernike polynomial functions (circle, astigmatism, coma). In addition, three different irradiation protocols were analyzed: standard Dresden CXL (fluence of 5.4 J/cm2), accelerated CXL (fluence of 5.4 J/cm2), and high-fluence CXL (fluence of 16.2 J/cm2). The axial strain distribution in the stroma induced by ocular inflation (Δp = 30 mmHg) was quantified, maps of the anterior sagittal curvature were constructed and cylindrical refraction was assessed. Results: Thirty minutes after CXL, there was a statistically significant increase in axial strain amplitude (p < 0.050) and a reduction in sagittal curvature (p < 0.050) in the regions treated with all irradiation patterns compared to the non-irradiated ones. Thirty-6 hours later, the non-irradiated regions showed compressive strains, while the axial strain in the CXL-treated regions was close to zero, and the reduction in sagittal curvature observed 30 minutes after the treatment was maintained. The Dresden CXL and accelerated CXL produced comparable amounts of stiffening and refractive changes (p = 0.856), while high-fluence CXL produced the strongest response in terms of axial strain (6.9‰ ± 1.9‰) and refractive correction (3.4 ± 0.9 D). Tripling the energy administered during CXL resulted in a 2.4-fold increase in the resulting refractive correction. Conclusion: OCE showed that refractive changes and alterations in corneal biomechanics are directly related. A patient-specific selection of both, the administered UV fluence and the irradiation pattern during CXL is promising to allow customized photorefractive corrections in the future.

Biaxial hyperelastic and anisotropic behaviors of the corneal anterior central stroma along the preferential fibril orientations. Part I: Measurement and calibration of personalized stress-strain curves

Lacking specimens is the biggest limitation of studying the mechanical behaviors of human corneal. Extracting stress-strain curves is the crucial step in investigating hyperelastic and anisotropic properties of human cornea. 15 human corneal specimens extracted from the small incision lenticule extraction (SMILE) surgery were applied in this study. To accurately measure the personalized true stress-strain curve using corneal lenticules, the digital image correlation (DIC) method and finite element method were used to calibrate the stress and the strain of the biaxial extension test. The hyperelastic load-displacement curves obtained from the biaxial extension test were performed in preferential fibril orientations, which are arranged along the nasal-temporal (NT) and the superior-inferior (SI) directions within the anterior central stroma. The displacement and strain fields were accurately calibrated and calculated using the digital image correlation (DIC) method. A conversion equation was given to convert the effective engineering strain to the true strain. The stress field distribution, which was simulated using the finite element method, was verified. Based on this, the effective nominal stress with personalized characteristics was calibrated. The personalized stress-strain curves containing individual characteristic, like diopter and anterior surface curvature, was accurately measured in this study. These results provide an experimental method using biaxial tensile test with corneal lenticules. It is the foundation for investigating the hyperelasticity and anisotropy of the central anterior stroma of human cornea.

Computational modeling of electromechanical coupling in human cardiomyocyte applied to study hypertrophic cardiomyopathy and its drug response

BACKGROUND AND OBJECTIVE

Knowledge of electromechanical coupling in cardiomyocyte and how it is influenced by various pathophysiological factors is fundamental to understanding the pathogenesis of myocardial disease and its response to medication, which is however hard to be thoroughly addressed by clinical/experimental studies due to technical limitations. At this point, computational modeling offers an alternative approach. The main objective of the study was to develop a computational model capable of simulating the process of electromechanical coupling and quantifying the roles of various factors in play in the human left ventricular cardiomyocyte.

METHODS

A new electrophysiological model was firstly built by combining several existing electrophysiological models and incorporating the mechanism of electrophysiological homeostasis, which was subsequently coupled to models representing the cross-bridge dynamics and active force generation during excitation-contraction coupling and the passive mechanical properties of cardiomyocyte to yield an integrative electromechanical model. Model parameters were calibrated or optimized based on a large amount of experimental data. The resulting model was applied to delineate the characteristics of electromechanical coupling and explore underlying determinant factors in hypertrophic cardiomyopathy (HCM) cardiomyocyte, as well as quantify their changes in response to different medications.

RESULTS

Model predictions captured the major electromechanical characteristics of cardiomyocyte under both normal physiological and HCM conditions. In comparison with normal cardiomyocyte, HCM cardiomyocyte suffered from systemic changes in both electrophysiological and mechanical variables. Numerical simulations of drug response revealed that Mavacamten and Metoprolol could both reduce the active contractility and alleviate calcium overload but had marked differential influences on many other electromechanical variables, which theoretically explained why the two drugs have differential therapeutic effects. In addition, our numerical experiments demonstrated the important role of compensatory ion transport in maintaining electrophysiological homeostasis and regulating cytoplasmic volume.

CONCLUSIONS

A sophisticated computational model has the advantage of providing quantitative and integrative insights for understanding the pathogenesis and drug responses of HCM or other myocardial diseases at the level of cardiomyocyte, and hence may contribute as a useful complement to clinical/experimental studies. The model may also be coupled to tissue- or organ-level models to strengthen the physiological implications of macro-scale numerical simulations.

A Critical Assessment of Friedenwald’s Technique for Estimating the Coefficient of Rigidity of the Cornea

Purpose

To determine if Friedenwald's technique for estimating the coefficient of corneal rigidity (Ko, units mmHg/μL), could differentiate between the cornea in keratoconus, normal eyes, and after crosslinking (CXL).

Methods

Two operators (1 and 2) independently measured Ko in three groups (keratoconus, normal, and post-CXL corneas), and repeated the procedure in some where their care remained unchanged and others after routine CXL (>28 days postop, epi-off treatment, 3.0 mW/cm², 30 min). The data were subsequently used to quantify interoperator error, test-retest/intersessional reliability for estimation of Ko, the significance of intergroup differences, and the effect of CXL on Ko.

Results

The major findings were: (i) Ko values were not normally distributed; (ii) mean (±sd, 95% CI) interoperator error was -0.002 (±0.019, −0.006 to 0.003, n = 95) and the limit of agreement between the operators was ±0.039; (iii) RMS differences in the intersessional estimation of Ko values were 0.011 (operator 1) and 0.012 (operator 2); (iv) intergroup differences in Ko were not significant (p > 0.05); (v) intersessional change in Ko (y) was linearly related to Ko estimated (x) at 1^st session (for operator 2 y = 1.187x−0.021, r = 0.755, n = 16, p < 0.01); and (vi) change in Ko (y₁) after CXL was linearly related to Ko (x₁) at preop (for operator 2 y₁ = 0.880x₁−0.016, r = 0.935, n = 20, p < 0.01).

Conclusion

Friedenwald's technique for estimating the Ko is prone to substantial interoperator error and intersessional differences. According to the technique, the change in Ko following CXL is on par with the expected intersessional change observed in controls.

Exploring the Biomechanical Properties of the Human Cornea In Vivo Based on Corvis ST

Purpose: The aim of this study was to provide a method to determine corneal nonlinear viscoelastic properties based on the output data of corneal visualization Scheimpflug technology (Corvis ST).

Methods: The Corvis ST data from 18 eyes of 12 healthy humans were collected. Based on the air-puff pressure and the corneal displacement from the Corvis ST test of normal human eyes, the work done by the air-puff attaining the whole corneal displacement was obtained. By applying a visco-hyperelastic strain energy density function of the cornea, in which the first-order Prony relaxation function and the first-order Ogden strain energy were employed, the corneal strain energy during the Corvis ST test was calculated. Then the work done by the air-puff attaining the whole corneal displacement was completely regarded as the strain energy of the cornea. The identification of the nonlinear viscoelastic parameters was carried out by optimizing the sum of difference squares of the work and the strain energy using the genetic algorithm.

Results: The visco-hyperelastic model gave a good fit to the data of corneal strain energy with time during the Corvis ST test (R² > 0.95). The determined Ogden model parameter μ ranged from 0.42 to 0.74 MPa, and α ranged from 32.76 to 55.63. The parameters A and τ in the first-order Prony function were 0.09–0.36 and 1.21–1.95 ms, respectively.

Conclusion: It is feasible to determine the corneal nonlinear viscoelastic properties based on the corneal contour information and air-puff pressure of the Corvis ST test.

Effective elastic modulus of an intact cornea related to indentation behavior: A comparison between the Hertz model and Johnson-Kendall-Roberts model

In this study, a macro-indentation test on the submillimeter scale was performed to analyze the indentation behavior of an intact cornea under physiological pressures. The Hertz and Johnson-Kendall-Roberts (JKR) models were employed to solve the elastic modulus (E) of the intact cornea. The relevant detailed analysis showed that the JKR model, which accounted for the contribution from the adhesion energy, could be used to obtain the E values that were more than two-folds of those obtained from the Hertz model, which only considered the external force. Compared with the uniaxial tension test in vitro, unlike the elastic Hertz-model, the E values under physiological pressures that were obtained with the JKR model were between the lower and upper limits of corneal material. This phenomenon indicated that the JKR model could be used to obtain reasonably effective E values of an intact cornea under physiological pressures.

Biomechanics of Ophthalmic Crosslinking

Crosslinking involves the formation of bonds between polymer chains, such as proteins. In biological tissues, these bonds tend to stiffen the tissue, making it more resistant to mechanical degradation and deformation. In ophthalmology, the crosslinking phenomenon is being increasingly harnessed and explored as a treatment strategy for treating corneal ectasias, keratitis, degenerative myopia, and glaucoma. This review surveys the multitude of exogenous crosslinking strategies reported in the literature, both "light" (involving light energy) and "dark" (involving non-photic chemical processes), and explores their mechanisms, cytotoxicity, and stage of translational development. The spectrum of ophthalmic applications described in the literature is then discussed, with particular attention to proposed therapeutic mechanisms in the cornea and sclera. The mechanical effects of crosslinking are then discussed in the context of their proposed site and scale of action. Biomechanical characterization of the crosslinking effect is needed to more thoroughly address knowledge gaps in this area, and a review of reported methods for biomechanical characterization is presented with an attempt to assess the sensitivity of each method to crosslinking-mediated changes using data from the experimental and clinical literature. Biomechanical measurement methods differ in spatial resolution, mechanical sensitivity, suitability for detecting crosslinking subtypes, and translational readiness and are central to the effort to understand the mechanistic link between crosslinking methods and clinical outcomes of candidate therapies. Data on differences in the biomechanical effect of different crosslinking protocols and their correspondence to clinical outcomes are reviewed, and strategies for leveraging measurement advances predicting clinical outcomes of crosslinking procedures are discussed. Advancing the understanding of ophthalmic crosslinking, its biomechanical underpinnings, and its applications supports the development of next-generation crosslinking procedures that optimize therapeutic effect while reducing complications.

Experimental in-vitro investigation on Epi-Off-Crosslinking on porcine corneas

Aim

To evaluate quantitatively the effects of the Epi-Off-CXL irradiance dose on the stromal stiffening of pig corneas.

Setting

Laboratory of Biological structures (LaBS), Politecnico di Milano, Milano, Italy.

Methods

Inflation tests have been carried on 90 excised and de-epithelized pig corneas, monitoring the change of configuration of the corneal dome at specific pressures. Test have been carried out twice on each cornea, once before and once after Epi-Off-CXL performed at a constant irradiance of 9 mW/cm² and variable UV-A exposure times. Corneas were grouped according to the exposure time (2.5, 5, 10, 15 and 20 min), proportional to the irradiation dose (1.35, 2.7, 5.4, 8.1, and 10.8 J/cm²). A theoretical model based on linearized shell theory has been used to estimate the increment of the corneal stiffness.

Results

The linearized shell theory allowed to establish a quantitative relation between the increment of the stiffness parameters and the irradiation dose. Relative to the pre-treatment values, in all experiments the post-treatment corneal stiffness revealed a pronounced increase. In general, the stiffness gain increased with the exposure time. No significant differences in stiffening was observed between tests conducted at 2.5, 5, and 10 min exposure.

Conclusions

Qualitatively, the effectiveness of accelerated CXL treatments observed in pig corneas complies very well with in-vivo clinical results in humans, suggesting that experimental data in pigs can be very useful for the design of the procedure in humans. A larger irradiation dose provides a larger increment of the corneal stiffness. Due to the biological variability of the tissues, however, it is difficult to distinguish quantitatively the level of the reinforcement induced by accelerated protocols (low doses with < = 10 min exposure), less prone to induce damage in the corneal tissue. Therefore, the definition of personalized treatments must be related to the actual biomechanics of the cornea.